A refresher in many-body theory 

Andre-Marie Tremblay 
May 2008 



2 



Contents 



1 Quantum mechanics for interacting systems in condensed matter 7 



1.1 Introduction 7 

1.2 Statistical Physics and Density matrix 7 

1.2.1 Density matrix in ordinary quantum mechanics 7 

1.2.2 Density Matrix in Statistical Physics 8 

1.2.3 Legendre transforms 8 

1.2.4 Legendre transform from the statistical mechanics point of 
view 9 

1.3 Second quantization 10 

1.3.1 Describing symmetrized or antisymmetrized states 10 

1.3.2 Change of basis 11 

1.3.3 Second quantized version of operators 12 

1.4 Hartree-Fock approximation 14 

1.4.1 The theory of everything 14 

1.4.2 Variational theorem 14 

1.4.3 Wick's theorem 15 

1.4.4 Minimization and Hartree-Fock equations 16 

1.5 Model Hamiltonians 17 

1.5.1 The Hubbard model 18 

1.5.2 Heisenberg and t-J model 18 

1.5.3 Anderson lattice model 20 

1.6 Broken symmetry and canonical transformations 21 

1.6.1 The BCS Hamiltonian 21 

1.7 Elementary quantum mechanics and path integrals 25 

1.7.1 Physical interpretation of the propagator 25 

1.7.2 Computing the propagator with the path integral 26 

1.7.3 Coherent-state path integrals 28 

2 Response functions and their general properties 31 

2.1 Relation between correlation functions and experiments 32 

2.2 Linear-response theory 32 

2.2.1 Schrodinger and Heisenberg pictures 33 

2.2.2 Interaction picture and perturbation theory 34 

2.2.3 Linear response 36 

2.3 General properties of correlation functions 37 

2.3.1 Notations and definitions 38 

2.3.2 Symmetry properties of the response functions 39 

2.3.3 Kramers-Kronig relations and causality 44 

2.3.4 Positivity of wx"(w) and dissipation 47 

2.3.5 Fluctuation-dissipation theorem 48 

2.3.6 Sum rules 49 

CONTENTS 3 



3 Green functions 55 

3.0.7 Second quantized operators in the Heisenberg picture ... 55 

3.1 Motivation of the definition of the second quantized Green's func- 
tion G R 56 

3.1.1 Examples with quadratic Hamiltonians: 58 

3.2 Interaction representation and time-ordered product 60 

3.3 Kadanoff-Baym and Keldysh-Schwinger contours 63 

3.4 Matsubara Green's function and its relation to usual Green's func- 
tions. (The case of fermions) 66 

3.4.1 Definition 66 

3.4.2 Antiperiodicity and Fourier expansion (Matsubara frequen- 
cies) 68 

3.4.3 Spectral representation, relation between G R and Q and an- 
alytic continuation 70 

3.4.4 Spectral weight and rules for analytical continuation .... 72 

3.4.5 Matsubara Green's function in momentum space and non- 
interacting case 74 

3.4.6 Sums over Matsubara frequencies 78 

3.5 Physical meaning of the spectral weight: Quasiparticles, effective 
mass, wave function renormalization, momentum distribution. ... 80 

3.5.1 Spectral weight for non-interacting particles 80 

3.5.2 Lehman representation 80 

3.5.3 Probabilistic interpretation of the spectral weight 82 

3.5.4 Angle-resolved photoemission spectroscopy (ARPES) on a 
Fermi liquid compound 83 

3.5.5 Quasiparticles [17] 86 

3.5.6 Fermi liquid interpretation of ARPES 88 

3.5.7 Momentum distribution in an interacting system 91 

3.6 A few more formal matters : asymptotic behavior and causality . . 92 

3.6.1 Asymptotic behavior of Q (k;ik n ) and £ (k;ife„) 93 

3.6.2 Implications of causality for G R and T, R 94 

3.7 Three general theorems 95 

3.7.1 Wick's theorem 96 

3.7.2 Linked cluster theorems 100 

3.7.3 Variational principle and application to Hartree-Fock theory 102 

3.8 Quantum impurities 104 



4 



CONTENTS 



List of Figures 



3-1 Kadanoff-Baym contour to compute G > (t — t') 65 

3-2 Keldysh-Schwinger contour 65 

3-3 Contour for time ordering in imaginary time 68 

3-4 Deformed contour used to relate the Matsubara and the retarded 

Green's functions 71 

3-5 Analytical structure of G(z) in the complex frequency plane. G(z) 
reduces to either G R (uS) , G A (oS) or Q (ioj n ) depending on the value 
of the complex frequency z. There is a branch cut along the real axis. 73 

3-6 Go (p, t) for a value of momentum above the Fermi surface 76 

3-7 Go (p, t) for a value of momentum at the Fermi surface 76 

3-8 Go (p, t) for a value of momentum below the Fermi surface 76 

3-9 Evaluation of fermionic Matsubara frequency sums in the complex 

plane 79 

3-10 Schematic representation of an angle-resolved photoemission exper- 
iment 83 

3-11 ARPES spectrum of 1 - T - TiTe 2 85 

3-12 Figure 1 from Ref.[27] for the ARPES spectrum of 1 T-TiTe 2 mea- 
sured near the Fermi surface crossing along the high-symmetry TM 
direction (6 — is normal emission) . The lines are results of Fermi 
liquid fits and the inset shows a portion of the Brillouin zone with 

the relevant ellipsoidal electron pocket 90 

3-13 Qualitative sketch of the zero-temperature momentum distribution 

in an interacting system 92 



LIST OF FIGURES 5 



LIST OF FIGURES 



1. QUANTUM MECHANICS FOR 
INTERACTING SYSTEMS IN CON- 
DENSED MATTER 



1.1 Introduction 



The purpose of these notes is to point out a few key results in quantum mechanics, 
statistical physics and many-body theory that should be useful to understand the 
application of advanced numerical methods to condensed matter. The degree of 
sophistication necessary to understand these methods is very high. The impor- 
tant developments in the field have been achieved by people familiar with all the 
anlytical tricks of the trade. 

These notes will not be very useful to learn the topics if you have never heard 
of them before. They should be seen as a summary of important concepts and 
methods to know in the context of this School. 



1.2 Statistical Physics and Density matrix 



1.2.1 Density matrix in ordinary quantum mechanics 

Quantum mechanics tells us that the expectation value of an observable O in a 
normalized state \ip) is given by (ip\ O \ip) . Expanding over complete sets of states, 
we obtain 

(vioiv) = ^2(ip\i)(i\o\j)(j\i>) (i.i) 
= ^2(jmmi)(i\o\j) (i.2) 

= ^2(j\p\i)(i\0\j) (1.3) 

i,i 

= Tr[pO] (1.4) 

where the Density Matrix p is defined, as an operator, by 

p = |V> <Vl - (1-5) 

This is when we have a pure state. If the state is prepared in a statistical 
superposition, in other words, if we have a certain probability p n that the state 
that is prepared is \ijj n ) , then the expectation value of an observable will be given 
by the weighted sum of the results in each state, in other words, in the above 
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formula for the average we should use 

n 

This is the density matrix for a mixed state. Note that 

P 2 =^2PnPm\lp n )('4>n\ V'mXV'ml- ( L7 ) 

We have the property p 2 = p only for a pure state. 

When a system of interest is in contact with an environment, it is very useful 
to work with an effective density matrix obtained by taking the trace first over the 
degrees of freedom of the environment. This idea is common in particular in the 
field of quantum information. In this school, we will see that by considering part 
of a large system as the environment, we can greatly reduce the size of the Hilbert 
space that needs to be considered to diagonalize a Hamiltonian, especially in one 
dimension. The optimal way of doing this was found by Steve White and will be 
discussed in the context of the "Density Matrix Renormalization Group". Not so 
surprisingly, quantum information theory has helped to improve even further this 
approach. Uli Schollwock will explain this. 



1.2.2 Density Matrix in Statistical Physics 



Statistical Physics tells us that conserved quantities play a special role. Indeed, 
at equilibrium, the density matrix cannot depend on time, so it depends only on 
conserved quantities. This means that generally, the density matrix is diagonal in 
the energy and number basis for example. All that is left to do is to specify p n . The 
basic postulate of statistical physics is that in an isolated system, all mcroscopic 
states consistent with the value of the conserved quantities are equiprobable. This 
is the microcanonical ensemble where p n is identical for all energy eigenstates 
\ip n ) ■ The other ensembles are derived in the usual way by considering the micro- 
canonical system as including the system of interest and various reservoirs. In the 
canonical ensemble for example, p n = e~@ En jZ where Z is the partition function 
Y. n er PEn ^df3 = (k B T)- 1 . 

Alternatively, the various ensembles are obtained by maximizing the entropy 

S = -k B Tr[plnp] (1.8) 

subject to constraints such as fixed average energy and normalization in the case 
of the canonical ensemble. Important properties of the entropy include extensivity 
and concavity. The entropy also plays a major role in quantum information. 



1.2.3 Legendre transforms 



Legendre transforms are encountered in mechanics when going from a Lagrangian 
to a Hamiltonian formulation. That transformation is extremely useful in sta- 
tistical physics as well and it will be used for example by Gabi Kotliar at this 
School. 

The important idea of statistical physics that we start with is that of potentials. 
If you know the entropy as a function of mechanical quantities, like energy volume 
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and number of particles for example, then you know all the thermodynamics. 
Indeed, 



dE = TdS - P dV + iidN (1.9) 
dS = Us+^dV-^dN (1.10) 

so you can obtain temperature T, pressure p and chemical potential fi simply by 
taking partial derivatives of the entropy. (1/T, S) (p/T, V) (—{i/T, N) are pairs of 
conjugate variables. Instead of using E, V, N as independent variables, given the 
concavity of the entropy and the uniqueness of the equilibrium state, you can write 
S as a function of any three other variables. Nevertheless, the purely mechanical 
variables E, V, N are the most natural ones for the entropy. The entropy plays the 
role of a thermodynamic potential. As a function of all microscopic variables not 
fixed by E, V, N, it is maximum at equilibrium. 

Remark 1 When there are broken symmetries, additional variables must be added. 
For example, for a ferromagnet with magnetization M in a magnetic field H, 

dE = TdS - pdV + fidN + MdH (1.11) 

There are other potentials. For example, if a system is in contact with a heat 
reservoir, the work that will be done at constant temperature will be modified by 
the presence of the reservoir. It is thus physically motivated to define for example 
the Helmholtz free energy 

F = E-TS (1.12) 




(1.13) 



In this case 

dF = dE- SdT - TdS = -SdT - pdV + fidN. (1.14) 

The Helmholtz free energy F can be written in terms of any three thermodynamical 
variables, but T, V, N are the most natural ones. At fixed T, V, N it is the free 
energy that is a minimum instead of the energy because we have to take into 
account the reservoir. The change from S to T as a natural variable has been 
done through the pair of equations (1.12,1.13). This is the general structure of a 
Legendre transform. F and E are potentials, and the subtraction of the product of 
the conjugate variables - (ff ) VJV S does the trick of relating the two potentials 

Remark 2 Note that (d 2 E/dS 2 ) = (dT/dS) = 1/ (dS/dT) = -1/ (d 2 F/0T 2 ) . 



1.2.4 Legendre transform from the statistical mechanics point of view 



Note that since 

T " T 



= S-% (1.15) 



= L ,1.10, 



dE 



V,N 



the quantity —FJT can be seen as the Legendre transform of the microcanonical 
entropy. From the point of view of statistical mechanics, if we define ft (E) as the 
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number of microstate n corresponding to a given energy, then p n = l/Vl{E) for 
every microstate and 

S(E) = -^Tr[pln^- fcs ^^ln^ (1.17) 

= k B lnn{E) (1.18) 

So, from the point of view of statistical mechanics, the Legendre transform of the 
entropy is obtained from 



-J = fc B lnZ = fc B ln^e-^" (1.19) 

n 

= k B ln^2n(E)e- 0E (1.20) 

E 

= k B \nJ2e lnn{E) e^ E (1.21) 

E 

= k B \nJ2e {S{E) - E/T)/kB (1.22) 

E 

Whereas the microcanonical entropy is a function of the energy of microstates, its 
Legendre transform is summed over energy and is a function of 1/T, the coefficient 
of E in both the thermodynamical expression of the Legendre transform Eq.(1.17) 
and the statistical one Eq.(1.22). 
We know that 



(E) = 91nZ 



d/3 
d{-F/T) 
d {1/T) 



(1.23) 



which clarifies the connection between the statistical mechanics and thermody- 
namical definitions of Legendre transform. E in the case of thermodynamics is 
really the average energy from the statistical mechanical point of view. The last 
equation could have been written down directly from the statistical expression for 
-FjT. 



1.3 Second quantization 



1.3.1 Describing symmetrized or antisymmetrized states 



States that describe identical particles must be either symmetrized, for bosons, 
or antisymmetrized, for fermions. To simplify the calculations, it is useful to 
use second quantization. As its name suggest, there is also an axiomatic way to 
introduce this method as a quantization of fields but here we will just introduce it 
as a calculational tool. The approach will be familiar already if you master ladder 
operators for the harmonic oscillator. 

For definiteness, let us concentrate on fermions. This can be translated for 
bosons. Define the operator ip^ (r) that creates a particle in a position eigenstate 
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|r) and antisymmetrizes the resulting wave function. Define also the vacuum state 
|0) that is destroyed by the adjoint, namely tp (r) |0) = 0. In this language 

^(r)V f (r')|0) = -^(|r)|r')-|r'>|r>) ( L2 4) 

= |r,r') = -|r',r). (1.25) 

The state to the right is clearly normalized and antisymmetric. There are two 
copies of the one-particle Hilbert space. In one component of the wave function, 
the particle in the first copy is at |r) , in the other component the particle in the 
first copy is at |r') . Clearly, that can become quite complicated. The two body- 
wave function (r, r' \ip) is antisymmetric and in the case where there are only two 
one-particle states occupied it is a Slater determinant. Clearly, that becomes a 
mess. In terms of the creation-annihilation operators however, all we need to know 
is that by definition of these operators, 

^ (r) V f (r') + i> ] (r') ^ (r) = 0. (1.26) 

We use the short-hand for anticommutation 

{^(r),^(r')} = 0. (1.27) 

Taking the adjoint, 

W(r)>(r')} = 0. (1.28) 

The only thing missing is that 

{V(r),V t (r')}=<5(r-r'). (1.29) 

That is a bit more complicated to show, but let us take it for granted. It is clear 
that if (r) creates a particle, then ip (r) removes one (or destroys it). If the 
particles are at different positions, that can be done in any order. If r = r', then 
it will matter if we create a particle before destroying it. If the creation occurs 
before the destruction, there will be one more particle to destroy. The Dirac delta 
function comes from normalization in the continuum. For discrete basis, we would 
have unity on the right. 



1.3.2 Change of basis 



A key formula for the "field" operators ip' (r) is the formula for basis change. 
Suppose that one has a new complete basis of one-particle states |a). Then, we 
can change basis as follows: 

|r)=5»<a|r> (1.30) 

a 

Given the definition of creation operators, the creation operator ip^ (r) for a par- 
ticle in state |r) is related to the creation operator for a particle in state |a) by 
the analogous formula, namely 

^(r)=$>t>|r>. (1.31a) 

a 

This formula is quite useful. 
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1.3.3 Second quantized version of operators 



One-body operators 

If we know the matrix elements of an operator in the one-particle basis, the cal- 
culation of any observable can be reduced to some algebra with the creation- 
annihilation operators. In other words, not only states, but also operators cor- 
responding to observables can be written using creation-annihilation operators. 
The expression for these operators is independent of the number of particles and 
formally analogous to the calculation of averages of operators in first quantized 
notation. 

To be more specific, consider the operator for the density of particles at position 
|r) . It can be written as ip^ (r) ip (r) as we prove now. Since ABC — CAB = 
ABC + ACB - ACB - CAB the commutator ot this operator with tp^ (r') is, 

^(v)i,(r),^(r')\ = ^(r){i,(r),^(r')}-{^(v),^ (r')} V (#-32) 

= S (r - r') V f (r) (1.33) 

We can now use the following little "theorem" on commutator of ladder operators: 

Theorem 1 If[A,B] = f3B and \a) is an eigenstate of A with eigenvalue a, then 
B \a) is an eigenstate of A with eigenvalue a+(3, as follows from AB \a)—BA \a) = 
A(B\a))-a(B\a)) = f3(B\a)). 

Since ip j (r) ip (r) |0) = 0, the above implies that ^ (r) ip (r) (V (r x ) |0)) = 
5 (r — ri) (ip^ ( r i) |0)^j , and generally a state ip^ (ri) ip^ (r 2 ) . . . |0) is an eigenstate 

of ip^ (r) ip (r) with eigenvalue S (r — ri) + S (r — r 2 ) + ... Clearly, the potential 
energy of identical electrons in a potential V (r) can ve written 

1^ (r)V(r)^(r)rf 3 r. (1.34) 

The same reasoning leads to the kinetic energy in the momentum basis, where it 
is diagonal 

h 2 k 2 , . d 3 k 



1 



c t(k)l^c(k)— 4r. (1-35) 



Returning to the position-space basis, we obtain 

/^(r)(-^)<Mr)d 3 r. (1.36) 

In other words, for any one-body operator, we can always obtain its second- 
quantized form in the one-particle basis |a) where it is diagonal: 

J2 cl (a| O \a) c a =J2 4 ("| O cp. (1.37) 

a a.{i 

If we change to an arbitrary basis 

|a)=X>)<i|a> (1.38) 

i 

the operator takes the form 

53 cl (a\ O c = H i) <*l O \j) (j\ a) c a = £ 4 (i\ O \j) Cj . (1.39) 

a,/3 oc.i.j i.j 
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Example 2 Let ip^ a (r) be the creation operator for the position state |r) with the 
spin a =|, I . We know the matrix elements of all component of the spin operators 
in the basis where S z is diagonal. Thus, from the last formula, we see that the 
three components of the spin operator are 

f ^(r)Q^^(r)d 3 r (1.40) 

where, as usual, the Pauli matrices are given by o z = ^ ^ ^ ^ , a v = ^ ? ^ ) . r ' = 

1 

1 



Two-body operators 

Let us now consider a two-body operator such as the potential energy. It is 
diagonal in position-space. The Coulomb interaction 

V c (r-r') = ] -^ 7] (1.41) 

is an example. The second quantized Coulomb energy takes the form 

J V c (r - r') \ (p (r) p (r') - 5 (r - r') p (r)) d 3 rdV (1.42) 

where the 1/2 comes from avoiding double-counting and 5 (r — r') p (r) is necessary 
not to count the interaction of an electron with itself. Including spin, the density 
operator is 

p(r)=£>£(r)^(r). (1.43) 

(7 

Substituting in the expression for the Coulomb interaction and using anti-commutation 
relations, we obtain 

\ E / V c ( r - r ') ^ ( r ) ^< ( r ') ^, ( r ) ( L44 ) 

It is an interesting and not very long exercise to prove that formula (which happens 
to have the same form for bosons and fermions). 

Let us change to some arbitrary basis. First notice that in terms of the potential 
energy operator V c 

V c (r-r') = (r|(r'|V c |r)|r'). (1.45) 

Then, the change of basis 

4(r) = ^4(i|r). (1.46a) 

i 

leads to the following two-particle analog of the one-body operator Eq.(1.39) in 
an arbitrary basis 

^EE(»i oi v c \k) io 44,c^ Cfe(7 . (i.47) 

a,a' ijkl 
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1.4 Hartree-Fock approximation 



The Hartree-Fock approximation is the simplest approximation to the many-body 
problem. It is a mean-field theory of the full Hamiltonian, that we will call "The 
theory of everything". We will begin by writing it explicitely then proceed with 
two theorems that form the basis of this approximation. 



1.4.1 The theory of everything 

Gathering the results of the previous section, an electron gas interacting with a 
static lattice takes the form 

Hioe = E / (r) + Vc,e-i (r)) * a (r) d 3 r 

+ ^E/ K(r-rVt(r)^(r')t,(r')t(r)dW (1.48) 

a, a' 

where V c ,e-i (r) is the electron-ion Coulomb potential. The dynamics of the ions 
(phonons) can be added to this problem, but until the rest of these introductory 
notes, we shall take the lattice as static. We need the to allow the lattice to move 
to have the complete "theory of everything" we want to solve in this School. But 
the above is certainly a non-trivial start. 



1.4.2 Variational theorem 



The Ritz variational principle states that any normalized wave function satisfies 

^\nW)>^\n\^) (i.49) 

where |^ ) is the ground state wave function. 

Proof. That follows easily by expanding \ip) = J2 i a, \ipi) , where H = 
Eilipi) , and using E < E t : 

> E ^H 2 = <V | |Vo> (1-50) 



In the Hartree Fock approximation, we use the variational principle to look for 
the best one-body Green function for H toe - In other words, we use our formula 
for a change of basis (there is no sum on repeated spin index here) 

^t( r ) = ^cL(i,a|r,,} = ^4X(r) (1.51) 

i i 

4 = / (r) (r,a| i, a) = J d 3 r^ ( r ) ^ ( r ) (i. 52 ) 
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and write our ground state wave function as 

iv-hf) = AAAAi ■ ■ ■ 4/2t4/ 2 i i°) • ( L53 ) 

Our variational parameters are the one-particle Green functions 0* (r) . Note that 
the most general wave function would be a linear combination of wave functions 
of the type \^ HF ) , each with different one-particle states occupied. 



1.4.3 Wick's theorem 

To compute (ip HF \ H IV'hf) > we expand each of the creation-annihilation opera- 
tors in the Hamiltonian Eq.(1.48) in the basis we are looking for, using the change 
of basis formula Eq.(1.51). Consider first the quadratic term and focus on the 
second quantized operators. We need to know 

\1>hf) (1-54) 

The key to compute such matrix elements is to simply use the anticommutation 
relations for the creation-annihilation operators and the fact that annihilation 
operators acting on the vacuum give zero. Let us do this slowly. 

The anticommutation relations for the operators c^J are as follows: 

{**,4*>} = /rf 3 r|d 3 r'C(r){^(r)>t,(r')}^(r') (1-55) 

= |d 3 rC(r)^ v (r) = ^W (1-56) 

so 

<0| c iT ct T |0> = 1 - <0| ct T c iT |0> = 1. (1.57) 

Generalizing this reasoning, we see that (ip HF \ i/>hf) = 1- Now, (tp HF \ q^-Cj-f \ipjjp) 
will vanish if either % or j are not in the list of occupied states in \tp HF ) since 
also annihilates the vacuum in the bra. If i and j are both in the list of occupied 
states, (iPhf\ c l] c j1 IV'hf) — since Cj-r will remove a particle in state j in \tp HF ) 
while c-| will remove a particle in state i in (ip HF | . If the list of particles is not the 
same in the bra and in the ket, the annihilation operators can be anticommuted 
directly to the vacuum and will destroy it. With this, we have that 

N/2 r f fi 2 v 2 \ 

= EE - 1 - + ^-,(r)U„(r)d 3 r. (1.59) 

a i=l ■' \ ' 

To compute the expectation value of the interacting part of H toe we need 

(^Hf\ c\ a c]a' C k<y'Cla \^Hf) ■ (1-60) 

Since \iPhf) ^ s a direct product of wave functions for up and down spins, if the 
spins are different, we obtain 

(tPhfI clAja'Cka'Cla \lp HF ) = 5^6 j}k . (1.61) 

If the spins are identical, something new happens. If the conditions k — I or 
i = j are satisfied, the expectation value vanishes because of the anticommutation 
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relations (Pauli principle). Consider k different from I. Since all we need is that 
the list of states created be the same as the list of states destroyed there are two 
possibilities 

{^HF\ 4a c )a c k<yCla \^ H f) = S i,l S j,k ~ S ilk Sj,l. (1.62) 

The last contribution is known as the exchange contribution. The difference in 
sign comes from the anticommutation. All these results, including the cases k = I 
or i = j for same spin, can be summarized by 

O'hfI c\a C ]a' C kcr'Cl<T \^Phf) = ^i,l^j,k ~ $i,kS j,l5„ ,a> ■ (1.63) 

The last result can be written as 

(4>HF\ 4a C ]a' C ka'Cla |V>ffF> = {^Hf\ 4a C lv \^Hf) {^Hf\ c]«t' C ^' I^HfK 1 -^) 

- (i>HF\4* C k<r' \^HF) {i ) HF\c)„>Clo 1^^-65) 

A four point correlation function has been factored into a product of two-point 
correlation functions. For states such as \iPhf) that are single-particle states, 
creation operators are "contracted" in all possible ways with the destruction op- 
erators. This elegant form is a special case of Wick's theorem. It applies to 
expectation values of any number of creation and annihilation operators. The 
signs follow from anticommutation. 



1.4.4 Minimization and Hartree-Fock equations 

Using Wick's theorem Eq.(1.63) and proceeding with the Coulomb interaction 
between electrons as we did with the one-body part of the Hamiltonian in Eq.(1.59) 
we obtain 

N/2 r ( h 2 v 2 \ 

(1> HF \ H toe \$ HF ) =^^2 <j>l (r) ( -— + V c , e -i (r)J <j> ia (r)d 3 r 

N/2 N/2 

+ E E E 2 / Vc (r r,) (r) *** (r) (r ' } (r ' } (L66) 

a, a' i—1 j — 1 

- S a , a ,<f>t (r) 4> ia (r') <f>* Ja , (r') <f> ja , (r)] d 3 rdV. (1.67) 

To find our variational parameters, namely the functions <p irj (r), we minimize 
the above, subject to the constraint that the wave functions must be orthonor- 
malized. This means that we take partial derivatives with respect to all variables 
in the above expression. We satisfy the constraints 

J C (r) <f> ja , (r) cPr-SijS^, = (1.68) 

using Lagrange multipliers. We have to think of (j)* a (r) and <\> ia (r) as independents 
variable defined at each different position r and for each index i, a. To take the 
partial derivatives carefully, one should discretize space and take the limit but the 
final result is pretty obvious. All we need to know is that what replaces the partial 
derivative in the continuum version is the functional derivative 

Wur (r) 
*<t> ja > (r') 

HI (r) 
H jtr > (r') 



5(r - r') SijSv^i 
0. 



(1.69) 
(1.70) 
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The result of the minimization with respect of <j)* a (r) is straightforward. One 
obtains 




d 3 r'V ex (r,r')<^ (r') 



N/2 



Ys^ii^a ( r ) 



(1.71) 




N/2 



dVV c (r-r') (r) 



(1.72) 



a' j=l 

N/2 



V ex (r,r') 



V c (t-t>)Y,4>% {*')<}> jo (r). 



(1.73) 



The matrix 7^ is a real symmetric matrix of Lagrange multipliers. Diagonalizing 
7jj and writing the eigenvalues £j, the above equation looks like a Schrodinger 
equation. The Hartree contribution Vh (r) has the physical interpretation that 
each electron interacts with the average density of the other electrons 



The exchange contribution V ex (r,r') has no classical analog. It comes from the 
anticommutation of indistinguishible particles. The E{ can be interpreted as single- 
particle excitation energies only if removing a particle does not modify too much 
the effective potentials. 



Suppose we have one-body states, obtained either from Hartree-Fock or from 
Density Functional Theory (DFT). The latter is a much better approach than 
Hartree-Fock. Nevertheless, it does not diagonalize the Hamiltonian. If the prob- 
lem has been solved for a translationally invariant lattice, the one-particle states 
will be Bloch states indexed by crystal momentum k and band index n. If we 
expand the creation-annihilation operators in that basis using the general formu- 
las for one-particle Eq.(1.39) and two-particle Eq.(1.47) parts of the Hamiltonian, 
clearly it will not be diagonal. Suppose that a material has s and p electrons, 
for which DFT does a good job. In addition, suppose that there are only a few 
bands of d character near the Fermi surface. Assuming that the only part of the 
Hamiltonian that is not diagonal in the DFT basis concerns the states in those d 
band, it is possible to write a much simpler form of the Hamiltonian. We will see 
that nevertheless, solving such "model" Hamiltonians is non-trivial, despite their 
simple-looking form. 

Model Hamiltonians can now explicitly be constructed using cold atoms in 
optical traps. A laser interference pattern can be used to create an optical lattice 
potential using the AC Stark effect. One can control tunneling between potential 
minima as well as the interation of atoms between them. 



N/2 




(1.74) 



1.5 Model Hamiltonians 
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1.5.1 The Hubbard model 



Restricting ourselves to a single band and expanding in the Wannier basis associ- 
ated with the Bloch states, the Hamiltonian takes the form 

H = E E 4 <*l K \j) c 3a + 01 V c \k) \l) cic\ a ,c la ,c ka . (1.75) 

a i,j a, a' ijkl 

where K contains all the one-body parts of the Hamiltonin, namely kinetic energy 
and lattice potential energy. The operator c^J annihilate (create) a particle in 
a Wannier state centered at lattice site i and with spin a. The one-body part 
by itself is essentially the DFT band structure. In 1964, Hubbard, Kanamori 
and Gutzwiller did the most dramatic of approximations, hoping to have a model 
simple enough to solve. They assumed that (i\ (j \ V c \k) \l) would be much larger 
than all other interaction matrix elements when all lattice sites are equal. Defining 
tij = (i\ K \j) and U = (i\ (i\ V c \i) \i) , and using Ci a Ci a — they were left with 



h — ^2 ^E **j4 c j> + 2 ^44 /C *CT' c i(7 
= EE *« 4 cjo- + E u 4 4 Ci i c ^T ( 1 • 76 ) 

<7 i , j i 
(7 i,j i 

Most of the time, one considers hopping only to nearest neighbors. The model can 
be solved exactly only in one dimension using the Bethe ansatz, and in infinite 
dimension. The latter solution is the basis for Dynamical Mean Field Theory 
(DMFT) that will be discussed at this School. Despite that the Hubbard model 
is the simplest model of interacting electrons, it is far from simple to solve. 

Atoms in optical lattices can be used to artificially create a system described 
by the Hubbard model with parameters that are tunable. The laser intensity of 
the trapping potential and the magnetic field are the control parameters. The 
derivation given in the case of solids is phenomenological and the parameters 
entering the Hamiltonian are not known precisely. In the case of cold atoms, one 
can find conditions where the Hubbard model description is very accurate. By 
the way, interesting physics occurs only in the nano Kelvin range. Discussing how 
such low temperatures are achieved would distract us to much. 

Important physics is contained in the Hubbard model. For example, the in- 
teraction piece is diagonal in the localized Wannier basis, while the kinetic energy 
is diagonal in the momentum basis. Depending on filling and on the strength of 
U compared with band parameters, the true eigenstates will be localized or ex- 
tended. The localized solution is called a Mott insulator. The Hubbard model 
can describe ferromagnetism, antiferromagnetism (commensurate and incommen- 
surate) and it is also believed to describe high-temperature superconductivity, 
depending on lattice and range of interaction parameters. 



1.5.2 Heisenberg and t-J model 



Suppose we are in the limit where U is much larger than the bandwidth. One 
expects that in low energy eigenstates, single-particle Wannier states will be either 
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empty or occupied by a spin up or a spin down electron and that double occupation 
will be small. If we could write an effective Hamiltonian valid at low energy, that 
means that we would reduce the size of the Hilbert space from roughly A N to 3^ 
for an N site lattice. This is possible. The effective Hamiltonian that one obtains 
in this case is the t — J model, which becomes the Heisenberg model at half-filling. 

To obtain this model, one can use canonical transformations or equivalently 
degenerate perturbation theory. Although both approches are equivalent, the one 
that is most systematic is the canonical transformation approach. Nevertheless, 
we will see a simplified version of the degenerate perturbation theory approach 
since it is sufficient for our purpose and simpler to use. 

We start from the point of view that the unperturbed part of the Hamiltonian 
is the potential energy. If there is no hopping, the ground state has no double 
occupancy and it is highly degenerate since the spins can take any orientation. 
Hopping will split this degeneracy. Let us write the eigenvalue problem for the 
Hubbard Hamiltonian in the block form 




(1.78) 



where Hu contains only terms that stay within the singly occupied subspace, H\2 
and H21 contains hopping that links the singly occupied subspace with the other 
ones and H22 contains terms that connect states where there is double occupancy. 
Formally, this separation can be achieved using projection operators. To project 
a state in the singly occupied subspace, one uses Hu = PHP where the projector 
P is 

N 

P = 11 (1 - n iT n a ) . (1.79) 

i=l 

Returning to the block form of the Hamiltonian, we can solve for Y = (E — H 22 y H 21 X 
and write 

(Hu + tfi2 (E - H22V 1 #21) X = EX. (1.80) 

What save us here is that the eigenstates we are looking for are near E = whereas 
H 2 2 will act on states where there is one singly occupied state since the hopping 
term in H\2 can at most create one doubly occupied state from a state with no 
double occupation. The leading term in H22 will thus simbply give a contribution 
U which is large compared to E. We are left with the eigenvalue problem 

<H n ^Hl) X = EX. (1.81) 

The first part of the Hamiltonian Hu contains only hopping between states 
where no site is doubly occupied. The potential energy in those states vanishes. 

The quantity H12H21 can be computed as follows. The only term of the original 
Hamiltonian that links singly and doubly occupied states is the hopping part. Let 
us consider only nearest neighbor hopping with = —t. Then 

#12^21 = ^ E E ( C U- + h.c) (cl,c la , + h.c) (1.82) 

(ij)a (kl)a' 

where each nearest-neighbor bond (ij) is counted only once in the sum. Since 
we leave from a state with singly occupied sites and return to a state with singly 
occupied sites, (hi) = (ij) survives as well as cases such as (kl) = (il) if one of 
the sites i is empty in the initial state. The latter contribution is called correlated 
hopping. It describes second-neighbor hopping through a doubly occupied state. 
In the t—J model, this term is often neglected on the grounds that it is proportional 
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to t 2 /U whereas Hn is of order t. That is not necessarily a good reason to neglect 
this term. Details of the derivation will follow in a subsequent version of these 
notes that will be on the archive of PITP. At half-filling we obtain 



= 4 



Uh 2 



^2 s * • s j + cst - 



{ij) 



This is the form of the Heisenberg Hamiltonian. In the case where we are not 
at half filling, with everything including the correlated hopping term, the t — J 
Hamiltonian takes the following form 



H 



+P 



<U> (ij) V 7 



4 i y C X-a C j-^ C \+k,a C 0+k',cr+C^ +k _ a c\_ a Cj^ cr CjJ r ki t a^j 



(1.83) 



where J = it 2 / (UH 2 ) and the last term is the three-site hopping term that is 
usually neglected. 

It is remarkable, but expected, that at half-filling the effective Hamiltonian is 
a spin-only Hamiltonian (The first term in the above equation does not contribute 
when there is no hole because of the projection operators). From the point of 
view of perturbation theory, the potential energy is the large term. We are in an 
insulating phase and hopping has split the spin degeneracy. 

Classically, the ground state on a hypercubic lattice would be an antiferromag- 
net. This mechanism for antiferromagnetism is known as super exchange. 

In closing, one should remember that to compute the expectation value of 
any operator in the singly occupied space, one must first write it in block form, 
in other words, one should not forget the contribution from the Y component 
of the wave function. For example, the potential energy of the Hubbard model 
calculated in the low energy subspace will be proportional to the exchange energy 
of the Heisenberg model. 



1.5.3 Anderson lattice model 



In the Anderson lattice model, on purely phenomenological grounds one considers 
localized states (/Jo-) wrtn a Hubbard U, hybridized with a conduction band (c^a^j 
of non-interacting electrons. This model is particularly useful for heavy fermions, 
for example, where one can think of the localized states as being / electrons: 

H A = Hf + H c + Hf c (1.84) 

(T i i 

H < = EE £ k4 CT c k . (1.86) 

a k 

E s- = EE^ +h - ft ( L87 ) 

a i 

In the case where there is only one site with / electrons, one speaks of the Anderson 
impurity model. When U is large, one can proceed as for the t— J Hamiltonian and 
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obtain an effective model where there is no double occupancy of the impurity and 
where the spin of the conduction electrons interacts with the spin of the impurity. 
The transformation is called the Schrieffer-Wolf transformation and the effective 
Hamiltonian is the Kondo Hamiltonian. 



1.6 Broken symmetry and canonical transformations 



The occurence of broken symmetry can be obtained from mathematical arguments 
only in very few situations, such as the Ising model in two dimensions. A simple 
paramagnetic state and a state with broken symmetry are separated by a phase 
transition, in other words by singularities in the free energy. Hence, the broken 
symmetry state cannot be obtained perturbatively. One postulates a one-body 
Hamiltonian where the symmetry is broken its stability verified using variational 
arguments. In this and many other contexts, canonical transformations are key 
tools to understand and solve the problem. We have seen examples above. Basis 
changes obtained from unitary transformations preserve the (anti) commutation 
relations. Such transformations are called canonical. We will illustrate these 
concepts with the example of superconductivity. 



1.6.1 The BCS Hamiltonian 

The general idea of Cooper pairs is that c P jct p ^ almost plays the role of a boson 
6 p . Commutation relations are not the same, but we want to use the general 
idea that superconductivity will be described by a non-zero expectation value of 
6 P by analogy to superfluidity. The expectation value ^ c pj c - p j.^) occurs in the 
Ginzburg-Landau theory as a pair wave function. The mean-field state will be 
described by a coherent state. 

We first write the general Hamiltonian in momentum space and, in the spirit 
of Weiss, the trial Hamiltonian for the mean-field takes the form 

H E -^N = H -tiN+yJ2 U (P-P')( c U c -Pl) c -p'l c p'l 

p,p' 

+ }7 E U (P - P') c pT c -p4 ( c -p'i c P'T> 
p>p' 

= # -^ + E( A P c -Pi c PT+ c PT ct -Pi A p) (1.88) 
p 

where we defined 

A P = ^E C/ (P-P')( c -P'i c P'T)- (1-89) 
p' 

The potential U (p — p') is an effective attraction that comes from phonons in 
standard BCS theory. We take this for granted. The states within an energy shell 
of size hw D around the Fermi level are those that are subject to that attraction. The 
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kinetic part of the Hamiltonian is given by 

H -iiN = ^(ffp-M)c P;CT Cp iCT (1.90) 

p,CT 

= ECp c p^ c p,- (i-9i) 

p,CT 

In the so-called jellium model, e p = fi 2 p 2 /2m e but one can take a more general 
dispersion relation. In matrix form, the combination of all these terms gives, 
within a constant 

£(4, -p0(1; (L92) 

One is looking for a canonical transformation that diagonalize the Hamiltonian. 
When this will be done, the c^ p ^ will be linear combinations of eigenoperators. 
These linear combinations will involve A p . To find the value of A p , it will suffice to 
substitute the eigenoperator expression for c pCT in the definition of A p , Eq.(1.89). 
This will give a self-consistent expression for A p . 

Let us define the Nambu spinor 

*p = ( T ) (1-93) 



c -pI 



whose anticommutator is 



{*P,i.*^. J -}=W*i,i (1-94) 



where i and j identiby the components of the Nambu spinor. Any unitary trans- 
formation of the Nambu spinors will satisfy the anticommutation relations, as one 
can easily check. Since the Hamiltonian matrix is Hermitian, it can be diagonalized 
by a unitary transformation. 

Eigenvalues E p are obtained from the characteristic equation 

(A P -Cp)(Ap + Cp)-|Ap| 2 =0 (1.95) 

where one used ( p = £_ p valid for a lattice with inversion symmetry. The solutions 
are 

A p = ±i? p = ± V / C P + |A p | 2 (1.96) 

and the eigenvectors obey 

-a; ±e p + c p J V «2 P J y ' 

whose solution is 

(±E p - C p ) ai P = A p a 2p (1.98) 
The constraint of normalization for a unitary transformation is 

|ai P | 2 + |a 2p | 2 = l. (1-99) 
The unitary transformation U 

u = (:i ) (i.ioo) 

u ] = ( u i y ) (i.ioi) 
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where 



p - 



V2 



(-a 



1/2 



1/2 



diagonalizes the Hamiltonian 



-E„ 



g»V2p 



17. 



Using this result, we can write 

H E = £ ( c pT c_ pi ) C/C/t ( C| ) [/t/t 

= E ( a PT a -Pi ) ( ^ 
p v 







c -p4 



(JL.102) 

(1.103) 
(1.104) 



where the new operators are related to the old by the Bogoliubov- Valentin (1958) 
transformation 



,t 

-pi 



a 



= c/t 



C -Pl 



-v* u„ 



c -p4 



(1.105) 



The ground state is the state that is annihilated by these new operators 

a pa \BCS) = 0. 

The new operators are linear combination of creation-annihilation operators since 
the eigenstate is a linear combination of states having different numbers of parti- 
cles. At zero temperature for example, one can check explicitely that the following 
state is indeed annihilated by a pa 



\BCS) = J] (l + VLt) l°>- 



The value of the gap A p is obtained from the self-consistency equation Eq.(1.89). 
It suffices to write the c p j en as a function of the diagonal operators a p(T . Inverting 
the Bogoliubov transformation Eq.(1.105) gives 



c -pl 



«pT 



whose adjoint is 
We also note that 



>i ) = ( a P T a -pl ) ( 



n(E p ) = (a pT a pT 



1 



p + i 



(1.106) 



(1.107) 



(1.108) 



The Fermi-Dirac distribution arises from the fact the the Hamiltonian is diagonal 
and quadratic when written as a function of fermionic operators a p ' . We can now 
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compute the mean value of the pair operator. 




( c -p'4 c p't) = ((v a P 'T + u P ,a - P 'l) ( u p'«p'T - v P ,a -p'l)) I 1 - 109 ) 

«-p'4«-p U ) (1-110) 

(1.111) 

l ^ P ' (l-2n(E p ,)) (1.112) 

(1 - 2n (E p ,)) . (1.113) 
2n(£ P 0) (1.114) 
Substituting in self-consistency equation, we Eq.(1.89) on obtain 

a p = ~w £ U (p - p0 §T (1 - 271 (i?p,)) ' (L115) 
p' p 

where A p is in general complex. This is known as the BCS equation. 

Remark 3 Even when the interaction depends on p - p', the phase is necessarily 
independent of p. Indeed, the gap equation can be rewritten in the form 

[C P A P ] = --L £ Cp U (p - p') C p , [C p , A p ,] . (1.116) 

w/iere 

27ie gap equation can then be reinterpreted as an eigenvalue equation. The eigen- 
vectors are in brackets and the eigenvalue is unity. Since the matrix — C P U (p — p') C p > / (2V) 
whose eigenvalues we are looking for is real and symmetric, the eigenvector is real 
within a global phase, i.e. a complex number e 1 ^ that multiplies all components 
of the eigenvector. This independence of p of the phase is known as "phase co- 
herence". It is key to superconductivity, If the eigenvalue of the gap equation is 
degenerate, something new can happen. One obvious degeneracy is associated with 
time-reversal symmetry. When this symmetry is broken, there is still an overall 
p independent phase, but the order parameter is complex in a way that does not 
correspond to a global phase. This in general gives, for example, a non-trivial 
value of the orbital angular momentum. 

Remark 4 Coherence: Since (f> lp + 4> 2p = 4> for all values of p, all the pairs are 
added to the wave function with exactly the same phase. This can be seen from 
the BCS wave function at zero temperature 



It is the interactions that impose that phase coherence that is at the origin of the 
phenomenon of superconductivity. Only the overall p independent phase of A is 
arbitrary. The global gauge symmetry is broken by fixing the phase since phase 
and number obey an uncertainty relation. Fixing the phase thus corresponds to 
making the total number of particles uncertain. 
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1.7 Elementary quantum mechanics and path inte- 
grals 



There is one formulation of quantum mechanics that is particularly useful for 
numerical simulations. It is the path integral formulation. The path integral 
formulation can be postulated. Quantum mechanics follows. Here we start from 
more familiar concepts and obtain the path integral formulation, as Feynman did. 

Let us go back to a single particle. The amplitude for a particle to go from 
position Xi to position Xf in a time t is given by 

G(x f ,t;xi,0) = (xf\e- iHt/n \x t ) . (1.118) 

This amplitude G (xf, t; Xi, 0) is known as the propagator. To simplify the notation 
we work in one-dimension and we work in units where H — 1. It is the basic object 
of this section. 



1.7.1 Physical interpretation of the propagator 

There are several ways to physically understand the propagator. From the basic 
postulates of quantum mechanics, squaring G (xj, t; Xi, 0) gives the probability 
that we are in eigenstate of position Xf at time t if the starting state is a position 
eigenstate a^. Also, if we know the propagator we know the amplitude to go from 
any state to any other one. Indeed, inserting complete sets of position eigenstates 
we find that 

(V/|e-""/ ft hk>= JdxidxfrfixriAixJixfle-^lxi) (1.119) 

Another way to see how to use the propagator is to define the retarded prop- 
agator 

G R (x f ,t;x i ,0) = (x f \e- im \x i )e(t) (1.120) 

where 9 (t) is the heaviside step function. Inserting a complete set of energy 
eigenstates, we find 

G R (x f ,t; Xi ,0) = Y,( x f\ n ) e ~ iEnt ( n \ x i) e{ f) 

n 

= J2^(xf)K(x*)e- lE " t 0(t) (1-121) 

n 

The Fourier transform of this quantity with r\ a positive real number is 

[°° dte^+^G* (x f ,t; Xi , 0) = i £ K (1.122) 
J -co n z + vrj hi n 

The poles of this function give the eigenenergies and the residues are related to 
the wave functions. In the many-body context, a generalization of the propagator 
occurs very naturally in perturbation theory. 

In statistical physics, (xf\p\xi) is a quantity of interest. Using the known form 
of the density matrix, we have (xf\p\xi) — {xf\e~P H \xi) /Z. Hence, computing 
these matrix elements is like computing the propagator in imaginary time r with 
the substitution t — > —it. This analogy holds also in the many-bady context. 
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1.7.2 Computing the propagator with the path integral 



In general, H contains non-commuting pieces. The potential energy V is diagonal 
in position space, but the kinetic energy K is diagonal in momentum space. Hence, 
computing the action of e~ lHt on \xi) is non-trivial since we need to diagonalize 
the Hamiltonian to compute the value of the exponential of an operator and that 
Hamiltonian contains two non-commuting pieces that are diagonal in different 
basis. The key observation is that if the time interval t is very small, say e, then 
the error that we do in writing the exponential as a product of exponentials is of 
order e 1 since it depends on the commutator of Ke with Ve 



e 



-iHs -iKe -iVs , r> (A 



+ 0(e 2 ). (1.123) 

In fact the error of order e 2 is in the argument of the exponential, as one can see 
from the Baker-Campbell-Hausdorff formula e A e B = e M with 

M = A + B + X - [A, B] + a 2 [A, [A, B]] + ... (1.124) 

In numerical calculations it is important to keep the exponential form since this 
garantees unitarity. 

Other factorizations give errors of even higher order. For example, 

e -iHe _ e -iVe/2 e -iKe e -iVe/2 ^ ^ 

gives an error of order e 3 . In practice, for numerical simulations it is quite useful 
to use factorizations that lead to higher order errors. To continue analytically 
however, the simplest factorization suffices. 

In the factorized form, we can take advantage of the fact that we can introduce 
complete sets of states where the various pieces of the Hamiltonian are diagonal 
to compute the propagator for an infinitesimal time 

(x f \e- iKe e- iVe \ Xi ) = J ^-(x f \e-^\p)(p\e-^\x t ) (1.126) 

= f ^P_^[-^+p^f-^)-^v{x i )] (i.i27) 

J 27T 

where we used (x\ p) = e lpx . The last formula can be rewritten 

e i(p±-H)s (1.128) 



/ 



2tt 

where 



x f 



(1.129) 

e 

The argument of the exponential is the Lagrangian times the time interval. It thus 
has the units of action and is made dimensionless by dividing by the quantum of 
action h that we have set to unity. 

For a finite time interval, we simply split the time evolution operator into 
evolution pieces that involve over an infinitesimal time interval 

N 

-iHt 



e 



= ]Je H£ (1.130) 



where e = t/N. There is no approximation here. Inserting N — 1 complete sets of 
states, we have 



, JV-l 



(x f \e lHt \xi) = / Yl dx j( x f\ e lHE \xN-i}(x N -i\e lHe \x N - 2 ) {x 2 \ ■ ■ ■ 

laiiXzile-"^!**). (1.131) 



26 QUANTUM MECHANICS FOR INTERACTING SYSTEMS IN CONDENSED MATTER 



Each of the N matrix elements can be evaluated now using the previous trick so 
that the propagator is given by the formally exact expression 



-iHt 



Xi) 



.N-l N 

a!™, / n dx > / n^ ex p 

J .7=1 J .7 = 1 



Xf - x N -i p% , ,\ 



X N -i-X N - 2 Pn-1 t/ ( 
~ -^^(V, 

+ ... 

/ Xx-Xi p\ \ 

To do the calculation, this is what one has to do. Formally however, the final 
expression is quite nice. It can be written as a path integral in phase space 

(x f \e- iHt \xi) = J [DxDp]exp{i J dt\px - H (p,x)]} (1.133) 

[DxDp] expiS (x, p) (1. 134) 



where the definition of the measure [DxDp] is clear by comparison and where S is 
the action. 

It is more natural to work in configuration space where the Lagrangian is 
normally defined. This comes out automatically by doing the integral over all 
the intermediate momenta. They can be done exactly since they are all Gaussian 
integrals that are easily obtained by completing the square 



f dp . ( X N -! - X N - 2 Pn-i , 

y-exp,^ — |: 



2-Kis 



2-Kis 



exp 



exp 



(1.132) 



im [ xn-i — xn-2 
irn 



l N-l L 



(1.135) 



Remark 5 The above is the propagator for a free particle. In that case, the time 
interval could be arbitrary and the result could also be obtained using our earlier 
decomposition on energy eigenstates since 



^2<f> n (x f )<j>* n ( Xi )e- iE ^= [ 



dPjp(x f -Xi)-it£ 

2tt 



(1.136) 



Once the integrals over momenta have been done, we are left with 



(xf\e lHt \xi) = lim 



r N-l 

■/n- 

J 3=1 



dxj 



mN 
2nit 



N 



exp 



+z£ Ul — * — 



V(x N _ 2 ) 



V x f 



+ is 



m ( x\ — Xi 



V(x t ) 



mx 2 — V (x) 



r x f 

JXi 



Dxe iS(x,x) 



(1.137) 
(1.138) 



where the formal expression makes clear only that it is the integral of the La- 
grangian, hence the action, that comes in the argument of the exponential. The 
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integration measure here is different from the one we had before. This form is 
particularly useful for statistical physics where all the integrals are clearly con- 
vergent, as opposed to the present case where they oscillate rapidly and do not 
always have a clear meaning. 

The physical interpretation of this result is quite interesting. It says that the 
amplitude for going from one point to another in a given time is given by the 
sum amplitydes for all possible ways of going between these two points in the 
given time, each path, or trajectory, being weighted by an exponential whose 
phase is the classical action measured in units of the action quantum H. The 
classical limit is obtained when the action is large compared with the quantum of 
action. Indeed, in that case the integral can be evaluated in the stationary phase 
approximation. In that approximation, one expands the action to quadratic order 
around the trajectory that minimizes the action. That trajectory, given by the 
Euler-Lagrange equation, is the classical trajectory according to the principle of 
least action. 



1.7.3 Coherent-state path integrals 

In the many-body context, the amplitudes that are interesting are of the form 



In the special case where only the ground state contibutes and that state is the 
vacuum state (i.e. no particle present), the above reduces precisely to our pre- 
vious definition since c\ |0) = \x,) and (0| e~ pH Cf (t) = (0| e~ 0H e iHt c f e- iHt = 
(0\c f e- iHt = (x f \e- iHt . 

To derive a path integral formulation for that type of amplitude, we note that 
destruction operators in H always appear first on the right. Hence, if we replace 
the position eigenstates in the one-particle case by eigenstates of the destruction 
operator, we will be able to derive a path integral formulation in the many-body 
case by following an analogous route. We will not do the full derivation here. The 
final result is that both for bosons and fermions, the path integral also involves 
exponentials of the action. For fermions, one must introduce Grassmann algebra 
with non-commuting numbers to define coherent states. For bosons the situation 
is simpler. 

Let us see how boson coherent states are constructed. Let [a, a*] = 1, then 
define the coherent state \z) by 



To show that this is an eigenstate of a, note first that one can easily show by 
induction that 



Tr pc f (t)4 



(1.139) 



\z) = e 



e za |0). 



(1.140) 




(1.141) 



which formally looks like 




(1.142) 



and since the exponential is defined in terms of its power series 



a, e 



z a 



Q e za 



= ze 



za 



(1.143) 
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Using our little theorem on commutators of ladder operators (1.3.2), we have that 
since a |0) = then a (e za] |0)^ = z [e zo ^ |0)^ and \z) is an eigenstate of a. 
To show that \z) is normalized, consider 

(z \z) = e-^ 2 (0| e z ' a e zaf |0) = (0| e za " |0) 

= 1 (1.144) 

In the last step, one has simply expanded the exponential in a power series and 
used the normalization of the vacuum. 
Finally we the closure relation 

I=\ J dzdz*\z)(z\ (1.145) 

that can be proven by taking matrix elements with states with arbitrary number 
of bosons \n) = (V) |0) j\fn\ and doing the integral in polar coordinates. 
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2. RESPONSE FUNCTIONS AND 
THEIR GENERAL PROPERTIES 



Whether perturbation theory is applicable or not, we rarely need all the informa- 
tion contained in the wave-function. A reduced description in terms of only a few 
variables suffices if it allows us to explain what can be observed by experimental 
probes. Correlation functions offer us such a description. As for any physical 
theory, we thus first discuss which quantities are observable, or in other words, 
what it is that we want to compute. 

Whatever the appropriate microscopic description of the system, or whatever 
the underlying broken symmetry, the result of any given type of experiment can 
be expressed as a specific correlation function. 

We will need to treat two different aspects of correlation functions. 

First, general properties, which are independent from the specific manner in 
which we compute correlation functions. For example 

• Symmetries 

• Positivity 

• Fluctuation-dissipation theorems relating linear response and equilibrium 
fluctuations 

• Kramers-Kronig transformations, which follow from causality 

• Kubo relations, such as that relating linear response to a specific correlation 
function. 

• Sum rules 

• Goldstone theorem, which follows from Bogoliubov inequalities 

Second, we will need to develop techniques to compute specific correlation func- 
tions. Sometimes, phenomenological considerations suffice to find, with unknown 
parameters, the functional dependence of correlations functions on say wave-vector 
and frequency. These phenomenological considerations apply in particular in the 
hydrodynamic regime, and whenever projection operator techniques are used. 

Microscopic approaches will lead us to use another type of correlation func- 
tions, namely Green's functions. In fact, Green's function are just one type of 
correlation function. They will appear very naturally. Furthermore, many of the 
general properties of correlation functions which we discuss in the present chapter 
will transpose directly to these functions. Much of this section is inspired from 
Foster. [1] 

In the present chapter, we intend to 

• Recall that scattering experiments are a measure of equilibrium fluctuations 

• Linear response to an external perturbation can be expressed as an equilib- 
rium correlation function 

And this correlation function can be related to equilibrium fluctuations by the 
fluctuation-dissipation theorem. 

• Then we discuss general properties of correlation functions 

• Give a specific example of sum-rule calculation. 
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2.1 Relation between correlation functions and ex- 
periments 



Scattering experiments with weak probes usually measure various equilibrium cor- 
relation functions of a system. What we mean by "weak probes" is simply that 
Fermi's Golden rule and the Born approximation are all that we need to describe 
the effect of the system on the external probe, and vice- versa. As an example, con- 
sider inelastic electron scattering but it should be clear that similar considerations 
apply to a large number of cases: inelastic light scattering, neutron scattering, 
etc... The plan is simply to use Fermi's Golden Rule to compute the cross section. 
One obtains 



da 
defdQf 



(2tt)3S5 ki 



\V1 



/d*e***(p q (t)p_ q (0)) 



(2.1) 



Forgetting for all the details, the key point is that the cross section is related to the 
Fourier transform of the density-density correlation function. The trick, due to Van 
Hove, to derive this formula from the Golden rule is to use the Dirac representation 
of the delta function for energy conservation and the Heisenberg representation to 
express the final result as a correlation function. Since in the Born approximation, 
incident and final states of the probe are plane waves, everything about the probe is 
known. The only reference to it will be through explicitly known matrix elements 
and quantum numbers, such as momentum, energy, spin etc... 

Definition 3 Note in passing that we use the following definitions for Fourier 
transforms in the continuum 



/ q = fd 3 r/(r)c 



-iqr 



(2.2) 



/(r)= f-Ms / q e iq ' r 



/ dt g{t)i 



g(t) = J %g u e-™ 



(2.3) 
(2.4) 
(2.5) 



2.2 Li near- response theory 



We are interested in the response of a system to a weak external perturbation. 
The electrical conductivity is the response to a weak applied field, the thermal 
conductivity the response to a thermal gradient etc... The result will be again 
an equilibrium correlation function. We will be able to relate this correlation 
function to equilibrium correlation functions of the type just calculated at the end 
of the last section by developing the so-called "fluctuation-dissipation theorem". 
The plan to compute the effect of an external perturbation is to add it to the 
Hamiltonian and then to treat it as a perturbation, taking the full interacting 
Hamiltonian of the system H as the unperturbed Hamiltonian. 
Let 

H (i) =H + SH(t) (2.6) 
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where Hq is the Hamiltonian of the system under study and SH(t) is the pertur- 
bation given by the time-dependent Hamiltonian 



5H(t)= - f d 3 r^(r)a,(r,i). (2.7) 



In this expression, Ai is some observable of the system (excluding external per- 
turbation) in the Schrodinger representation, while a,i(r,t) is the external field. 
Examples of such couplings to external fields include the coupling to a magnetic 
field h through the magnetization M, (Ai (r) = M z (r) ; aj(r,i) = h z (r, t)) or 
the coupling to an electromagnetic vector potential A/c through a current j, 
(Ai (r) = j x (r)5i,x', a>i( r ,t) = A x (r,t)5i iX /c) or that of a scalar potential <fi through 
the density p (Ai (r) = p(r) ; di(r,t) = (f>(r,t)). In this approach, it is clear that 
the external perturbation is represented in the semi-classical approximation, in 
other words it is not quantized. We first pause to recall the various representa- 
tions, or pictures, of quantum mechanics, introducing the interaction representa- 
tion as the framework where perturbation theory is most easily formulated. Then 
we go on to derive linear response theory. 



2.2.1 Schrodinger and Heisenberg pictures 

Since the Hamiltonian is the infinitesimal generator of time translations, Schrodinger's 
equation for a time-dependent Hamiltonian takes the same form as usual, 

iH lh = Hms - (2 ' 8) 

Using the fact that H(t) is Hermitian, one can easily prove that d (tp s \tp s ) /dt = 0, 
in other words that probability is conserved. Hence, the solution of this equation 
will be given by 

i> s (t) = U(t,t )ips(to) (2.9) 
where U(t,to) is a unitary operator satisfying 

U(to,to) = l (2.10) 

while by time-reversal symmetry 

U(t ,t)U(t,t ) = l. (2.11a) 

Conservation of probability gives 

U(t,t )^U(t,t ) = l (2.12) 

so that combining the last result with the definition of the inverse, we have, 

U(t,t r l =U(t,t Q ?. (2.13) 

Furthermore, when we can use time-reversal invariance, Eq.(2.11a), we also have 

U^to)- 1 = U(t,t ) ] =U(t ,t). (2.14) 

By definition, for all values of t, the expectation value of an operator is the 
same in either the Schrodinger, or the Heisenberg picture. 



Ms (t)\ O s |Vs (t)> = 0>h\ Oh (t) |<M . (2.15) 
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In the Heisenberg picture the operators are time-dependent while in the Schrodinger 
picture, only the wave functions are time dependent. Let us choose t = to be 
the time where both representations coincide. The choice of this time is arbitrary, 
but taking t = simplifies greatly the notation. We have then that 

O s (t = 0) = O H (t = 0) = O s (2.16) 

iP s (t = 0) = i> H (t = 0) = ^ s (2.17) 

Using the expression for the time-dependent wave function, and the equality of 
matrix elements Eq.(2.15), we obtain 

O H (t) = U\t,0)O s U(t,0). (2.18) 

One recovers all the usual results for time-independent Hamiltonians by noting 
that in this case, the solution of Schrodinger's equation is, 



U(t,t ) = e -i«(t-to)/fi. (2.19) 



Remark 6 When there is time-reversal invariance, then it is useful to replace the 
adjoint by the time-reversed operator, so that the connection between Heisenberg 
and Schrodinger picture Eq.(2.18) becomes 

O H (t) = U(0,t)O s U(t,0). (2.20) 

Because we do not want to assume for the time being that there is time-reversal 
invariance, we shall stick here with the usual expression Eq.(2.18) but in much of 
the later chapters, the above representation will be used. 



2.2.2 Interaction picture and perturbation theory 



Perturbation theory is best formulated in the "interaction representation" . In this 
picture, one can think of both operators and wave functions as evolving, as we 
will see. We take 

(2.21) 



H (t) = H + 8H(t) 



where Ho is time-independent as above, but the proof can be generalized to time- 
dependent Ho simply by replacing e tH o t / h everywhere below by the appropriate 
evolution operator. 

The definition of the evolution operator in the interaction representation Ui(t, 0) 
is given by 

U(t, 0) = e- iHnt / h Ui(t, 0). (2.22) 



and 

so that for example 



17(0, t) = [//(Me 



iHat/h 



U(t,t Q ) = e- lH " t / h U I (t,to)e lH ° t °/ h 



(2.23) 



(2.24) 



We have used the fact that Uj(t, to) obeys the same general properties of unitarity 
as an ordinary evolution operator, as can easily be checked. Again the interaction 
representation will coincide with the other two at t = 0. The justification for 
the definition of Uj above is that when the external perturbation 5H(t) is small, 
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Ui(t, to) is close to unity. If we write again the equality of matrix elements in the 
general case, we obtain 

(1> s (t)\O s |Vs (*)} = (ip s \ U\t, 0)O s U(t,0) \ip s ) (2.25) 

= (iP s \ Uj(t,0)e m " t / h O s e- iHot / h U I (t,0) |V S ) (2.26) 

= ^ s \U\{t,Q)0 I {t)U I {t,Q)\i; s ) (2.27) 

This last result is important. It can be interpreted as saying that the operators 
in the interaction representation evolve with 

(2.28) 



Oi (t) = e iH Q t/n Ose -iH t/n 



while the wave functions obey 



|Vz (t)) = Uj(t,0)\^ s ) 



(2.29) 



In other words, in the interaction picture both the operators and the wave function 
evolve. We still have to find the equation of motion for £//(f, to)- The result will 
justify why we introduced the interaction representation. Start from Schrodinger's 
equation, 



ih- 



dt 



(2.30) 



which gives the equation of motion for Ui(t, 0), namely 

#oe~^ ot/s £//(i,0) +e- tHot / h ih^U I (t,0) = n(t) e - lH ° t / h U I (t 1 0) 

ih^Uj^O) = e iHot / h 5H(t)e- iHot / n Ui(t,0). 

so that using the definition of time evolution of an arbitrary operator in the inter- 
action representation as above (2.28) the equation for the time evolution operator 
Ui(t,0) in the interaction representation may be written, 



(2.31) 
(2.32) 



ift^C/}(t,0) = 6Hi{t)U I {t,Q) 



(2.33) 



with the initial condition 

£7/(0,0) = 1. (2.34) 

As expected, Eq.(2.33) tells us that, if there is no perturbation, Uj is equal to 
unity for all times and only the operators and not the wave function evolve. The 
interaction representation then reduces to the Heisenberg representation. Multi- 
plying the equation of motion from the right by Ui(0, to) we have for an arbitrary 
initial time 

'ifc£E//(Mo) = $W/(t)£Mt,to)| ( 2 - 35 ) 



We will come back later to a formal solution of this equation. To linear order 
in the external perturbation, it is an easy equation to solve by iteration using 
the initial condition as the initial guess. Indeed, integrating on both sides of the 
equation of motion (2.35) and using the initial condition, we have 



Uj(t,t ) = 1 - * /* dt> SHjit'W,^) 



which, iterated to first order, gives, 

17/(4, to) = 1 - I I dt' 5Hi(t') + 0{SHj) 

» J to 



(2.36) 



(2.37) 



and correspondingly 



■ ,-t 

Uj(t,t ) = 1 + - / dt' 8Hi(t') + O(SHj) (2.38) 

" J to 
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2.2.3 Linear response 



Returning to our general considerations, in the case of interest to us the external 
perturbation in the interaction representation is of the form, 



(2.39) 



where for short we wrote A®(r, t) to represent a system's observable evolving in 
the system's Heisenberg representation, 



A?(r,t) =e lHat ' h A l {v)e- lHat/h . 



(2.40) 



Suppose we want the expectation value of the observable B in the presence of 
the external perturbation turned on at time to- Then, starting from a thermal 
equilibrium state at time t , it suffices to evolve B (r) with the full evolution 
operator, including the external perturbation. 1 

(B(r,t)) = (U\t,t )B(r)U(t 7 t )) (2.41) 

Using the interaction representation Eq.(2.24), the last equation becomes 

(B(r,t)> = (^(Mo)£°M)tMMo)) • (2.42) 

In this last expression, B°(r, t) is now in the system's Heisenberg representation of 
the system without the external perturbation. We also used the fact that e ~ lHata l h 
commutes with the density matrix and that the trace has a cyclic property to 
cancel the extra e -* H o*o/ fi anc j e iH t /h appearing in the equation for the evolution 
operator in Eq.(2.24). 

Using the explicit expression Eq.(2.39) for the external perturbation in the 
equation for the evolution operator in the interaction representation (2.37), we 
have that the term linear in applied field is then given by, 



S(B(r,t)) = 



h 



to 



dt' / d 3 r'([B a (r,t),A^r',t')]) ai (r',t'). 



(2.43) 



It is customary to take to = — oo, assuming that the perturbation is turned-on 
adiabatically slowly. One then defines a "retarded" response function, or suscep- 
tibility x R , by 





S(B(v,t)} = 


J^dt'Jd'r' X £4,M;r',i'K(r'/) 


with, 










t') = j-([B%r,t),A^r',t')})6(t-t'). 



(2.44) 



(2.45) 



This response function is called "retarded" because the response always comes after 
the perturbation, as expected in a causal system. The function 6(t — t') ensures 
this causality. One can also define anti-causal response functions. We come back 
to this later. For the moment, recall that the superscript here means to zeroth 
order in the external probe. In other words, the linear response is given by an 
equilibrium correlation function. One normally does not write the superscript 
which is usually kept to mean non-interacting system. From now on, we drop this 
superscript. 



*We let the density matrix take its initial equilibrium value. This is physically appealing. But 
we could have as well started from a representation where it is the density matrix that evolves 
in time and the operators that are constant. 
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Remark 7 Translationally invariant case: Since we compute equilibrium aver- 
ages, the susceptibility XBAi ( r ' *i r '' O can depend only on the time difference. 
In the translationally invariant case, the susceptibility is also a function of only 
r r' so iftai Fourier transforming the expression for the linear response (2.44), 
we obtain from the convolution theorem in this case, 



5{B( q ,u)) =xi Ai (q,wK(q,w). (2.46) 



Remark 8 Frequency of the response: The response is at the same frequency as 
the external field, a feature which does not survive in non-linear response. 

Remark 9 Onsager reciprocity relations: Given the expression for the response 
function in terms of a commutator of Hermitian operators, it is clear that the 
response of the operator B to an external perturbation that couples to A is simply 
related to the response of A to a perturbation that couples to B, in other words 
where the operators have reversed roles. These are "Onsager 's reciprocity rela- 
tions". 

Remark 10 Validity of linear response and heating: Finally, we can ask whether 
it is really justified to linearize the response. Not always since the external pertur- 
bation can be large. But certain arguments suggest that it is basically never correct 
in practice to linearize the response. Indeed, assume we apply an external electric 
field E. As long as the energy gained by the action of the field is smaller than 
k B T, the linearization should be correct. In other words, linear response theory 
should be valid for a time 

' < < 2 -«» 

This is unfortunately a ridiculously small time. Taking v « ^Jk B T jm the condi- 
tion becomes t < y/mk B T/eE with E = IV '/cm, ^mk B T/eE « VlO- 30 10- 23 10 2 /KT 
10~ 6 s. Indeed, one finds that unless there is a temperature gradient, or an explicit 
interaction with a system in equilibrium (such as phonons), the second order term 
in perturbation theory is secular, i.e. it grows linearly with time. This is nothing 
more than the phenomenon of Joule heating. [2] We are then forced to conclude 
that linear response theory applies, only as long as the system is maintained in 
equilibrium by some means: for example by explicitly including interactions with 
phonons which are by force taken to be in thermal equilibrium, or by allowing for a 
thermal gradient in the system that carries heat to the boundaries. In a Boltzmann 
picture, one can see explicitly that if the second-order term in E is kept small by 
collisions with a system in thermal equilibrium, then the linear term is basically 
equal to what we would have obtained by never going to second-order in the first 
place. [2] 

Remark 11 Reversibility and linear response: Other arguments against linear 
response theory center on the fact that a correlation function where operators all 
evolve reversibly cannot describe irreversible processes. [3] We will see explicitly 
later that it is possible to compute irreversible absorption with this approach. We 
will also see how irreversibility comes in the infinite-volume limit. 



2.3 General properties of correlation functions 



It is useful to know analytic properties that do not depend on the microscopic 
model considered. This has at least two advantages: a) to check whether approx- 
imation schemes satisfy these exact relations b) to formulate phenomenological 
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relations which are consistent. We will see that approximate calculations cannot 
satisfy all known exact relations for correlation functions, but it will be obvious 
that violating certain relations is more harmful than violating others. Many of 
the general properties which we will discuss in the present context have trivial 
generalizations for Green's function. Working on these general properties now 
will make them look more natural later when we introduce the curious Green's 
function beast! 



2.3.1 Notations and definitions 



To start with, recall the definition 



X ^(r, t; r', t>) = - ([B(r, t),A(r', f )]> 6{t t>). 



(2.48) 



We have removed the superscript which was only used as a crutch in the deriva- 
tion of linear-response theory to indicate that the operators were evolving with the 
unperturbed Hamiltonian. Since the unperturbed Hamiltonian in the present con- 
text is the full Hamiltonian of the system, including interactions, we will drop the 
superscript from now on. It will be used later in a context where the unperturbed 
Hamiltonian is that of non-interacting particles. 

We define one more correlation function which will, in most cases of physical 
interest, play the role of the quantity which describes absorption. Welcome \" 



X 'Uv,fy,t') = ^([B(r,t),A(r',t')]). 



(2.49) 



The two in the denominator looks strange, but it will allow x" to generally be 
the imaginary part of a response function without extra factors of 2. With this 
definition, we have 



xf A (r, t; r', t') = 2*x£ A (r, t; r', t')6{t - t'). 



To shorten the notation, we will also use the notation 



xl A it- t ') = i{[Mt),A ] {t')})e{t-t>). 



(2.50) 



(2.51) 



In this notation we include in the indices i and j the positions as well as any other 
label of the operator such as vector or spin component. In this notation, we have 
not assumed translational invariance. We did however assume time-translation 
invariance. Since we are working with equilibrium averages above, this is always 
true. 



Exercise 2.3.1 Check time-translational invariance explicitly by using Heisen- 
berg's representation, the cyclic property of the trace and the fact that the den- 
sity matrix (Z~ 1 e~ l3H in the canonical ensemble, or »" 1 e _ ^ ff_ ' lJV ) in the grand- 
canonical) commutes with the time- evolution operator e - lHt / h . 

Corresponding to the short-hand notation, we have 



(t-t) = ±([Mt),M*)]) 



xlA j (t-t') = 2ix , AiAj (t-t')e(t-t'). 



(2.52) 
(2.53) 
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2.3.2 Symmetry properties of the response functions 



The quantity Xaa- contains all the non-trivial information on the response. 

Indeed, the causal response is simply obtained by multiplying by a trivial 9(t — t') 
function. Certain symmetries of this response function depend on the particular 
symmetry of the Hamiltonian, others are quite general. Let us consider them in 
turn.[l] 

Properties that depend on the symmetry of the Hamiltonian 

Let S be a symmetry of the Hamiltonian. By this we mean that the operator S 
representing the symmetry commutes with the Hamiltonian 

[H, S} = (2.54) 

To be more precise, in the context of statistical mechanics we say that S is a 
symmetry of the system when it commutes with the density matrix 

[p,S]=0 (2.55) 

In other words, S~ 1 pS — p thus the spectrum of the density matrix is unaffected 
by the symmetry operation. The operator S is in general unitary or antiunitary 
as we will see below. 

To extract non-trivial consequences of the existence of a symmetry, one first 
takes advantage of the fact that the trace can be computed in any complete basis 
set. This means that the thermal average of any operator O is equal to its thermal 
average in a basis where the symmetry operation S has been applied to every 
basis function. Since the symmetry operation commutes with the density matrix 
by assumption, one can then let the symmetry operations act on the operators 
instead of on the basis functions. In other words, we have 



(S^OS) = (O) (2.56) 



It is because S and O in general do not commute that the above equation leads 
to non-trivial consequences. 

Let us look in turn at the consequences of translational invariance and of 
invariance under a parity transformation r a -^ — r a . We skip the proofs. 

• Translational invariance: 



X B A (r, t; r',t') = X ba^ + R, *; + R, O 



(2.57) 



so that Xba i s a function of r — r' only. Since we already know that Xba 18 
a function only oit — t', in such cases we write 



X BA (r,t;r',t') = XBA(r-r';t-t') 



(2.58) 



• Parity: Under a parity transformation, operators transform as follows 

p- 1 0(r)P = e p O(-r) (2.59) 

where e p = ±1. This number is known as the "signature" under parity 
transformation. That e p = ±1 is the only possibility for simple operators 
like density and momentum follows from the fact that applying the parity 
operation twice is the same as doing nothing. In other words, P 2 = 1. 
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When we also have translational invariance, the last result means that x'ba( v ~ 
t') is even or odd in r — r' depending on whether the operators have the same 
or opposite signatures under parity. Correspondingly, the Fourier transform 
in the translationally invariant case is odd or even, as can easily be proven 
by a change of integration variables in the Fourier transform 



• Time-reversal symmetry in the absence of spin: From the Schrodinger equa- 
tion in the absence of spin, one can see that when the Hamiltonian is real, 
then complex conjugation leads to an equation that evolves the complex con- 
jugate wave function as if t — > —t. We thus take time-reversed states as just 
this operation of complex conjugation. A system in equilibrium obeys time- 
inversion symmetry, unless an external magnetic field is applied. This means 
that equilibrium averages evaluated with time-reversed states are equal to 
equilibrium averages evaluated with the original bases. In fact time-inversion 
symmetry is a very subtle subject. A very complete discussion may be found 
in Gottfried [4] and Sakurai [8] . We present an oversimplified discussion. Let 
us call T t the operator that time-reverses a state. This is the operation of 
complex conjugation that we will call K. The first thing to notice it that it 
is unlike any other operator in quantum mechanics. In particular, the Dirac 
notation must be used with extreme care. Indeed, for standard operators, 
say X, we have the associative axiom 



This is clearly incorrect if X is the complex conjugation operator. Hence, 
we must absolutely specify if it acts on the right or on the left. Hence, we 
will write K when we want to take the complex conjugate of a ket, and K 
to take the complex conjugate of a bra. 

Remark 12 Antiunitary operators: Time reversal is an antiunitary oper- 
ation. The key property that differentiates an anti-unitary operator from a 
unitary one is its action on a linear combination 



In general such an operator is called antilinear. Antiunitarity comes in when 
we restrict ourselves to antilinear operators that preserve the norm. The time 
reversal operator is such an operator. Under time reversal, an arbitrary ma- 
trix element preserves its norm, but not its phase. This is easy to see from 
the fact that for an arbitrary matrix element 1 KK \tp 2 ) — (ip 2 IV'i) 
(ip 1 \ip 2 ) the phase changes sign under complex conjugation while the square 
modulus (if) 2 IV'i) (^i 1^2) * s invariant. Gottfried^] shows that only discrete 
transformations (not continuous ones) can be described by anti-unitary op- 
erators. This reference also discusses the theorem by Wigner that states 
that if we declare that two descriptions of quantum mechanics are equiva- 
lent if \(ip 2 |i/> 1 )| = \(i>2 (equality of "rays") then both unitary and 
anti-unitary transformations are allowed. 



Xba(<1; t-t')= e p B s p A x' BA (-V, t - t') 



(2.61) 



(a\X\(3) = (a\ (X\(3)) = ((a\X) \[3) 



(2.62) 



T t (oi |Vi> + a 2 IV2)) = a\T t |^> + a* 2 T t \^ 2 ) 



(2.63) 
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Remark 13 The adjoint is not the inverse. Note that T^T t = KK , so this 
last quantity is not the identity because the rightmost complex conjugation 
operator acts to the right, and the leftmost one to the left. Again, it is not 
convenient to talk about time-reversal in the usual Dirac notation. 



• Returning to the action of the time reversal operation on a Schrodinger 
operator, we see that the expectation value of an arbitrary operator between 
time reversed states is 

(i\ KOK \j) = (<i| K) [KO* \j)) = «i| O* \j))* = (j\ 0t* \i) . (2.64) 

In the above expression, we used one of the properties of the hermitian 
product, namely (k \l)* — (I \k)* , as well as the definition of the adjoint of 
an operator A: (k\A \l)* = (l\A^ \k) . Applying this expression Eq.(2.64) for 
diagonal expectation values, and recalling that the density matrix is real, we 
find for equilibrium averages, 

^KOK^ = {0^)=e t {0^) (2.65) 

The last equality defines the signature of the time-reversal operation for 
operators. One easily finds that e 4 = +1 for position while e* = — 1 for 
velocity or momentum, etc... We can use this last result to find the effect of 
the time-reversal invariance on general correlation functions. The action of 
time reversal Eq.(2.65) gives, when A and B are self-adjoint operators, and 
KB = HK 

I^KA(t)BK^ = ^B*e- lHt / n A*e lHt / h ^ 

= e^iBAi-t)) (2.66) 

In addition to the signature, the order of operators is changed as well as the 
sign of time. For X%a (* - O this immediately leads to 

X%a 3 it t>) = e\e]^ AjAi (-f - (-*)) (2.67) 
and for the corresponding Fourier transform in frequency, 



XW W ) = e HxW") • (2-68) 



Remark 14 In the case of an equilibrium average where both the density 
matrix and the Hamiltonian commute with the time-reversal operation, we 
have, as in Eq. ( 2. 56), 

(Tf 1 OT t ) = (O) . (2.69) 

Hence as expected, Eqs.(2.65) and (2.69) together imply that operators that 
have an odd signature with respect to time reversal symmetry have a vanish- 
ing expectation value in equilibrium. 



• Time-reversal symmetry in the presence of spin: Spin should transform un- 
der time reversal like angular momentum r x p, in other words it should 
change sign since r does not while p does. Complex conjugation has this 
property for r x p but not for spin represented by Pauli matrices. We should 
really wait for the section where we treat fermions to discuss this problem 
but we can start to address it here. To come out from the problem that 
complex conjugation does not suffice anymore, it suffices to notice that in 
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general the time reversal operator has to be represented by a unitary oper- 
ator times complex conjugation. The resulting operator is still anti-unitary, 
as can easily be proven. Let us thus write 



KU 



(2.70) 



where K is complex conjugation again and U is a unitary operator U^U = 1 
in spin space that we need to find. Note that the action on a bra is given by 



(2.71) 



Let us first repeat the steps of calculating expectation values in time-reversed 
states, as in Eq.(2.64), but for the more general case 

(i\ U^KOKU \j) = (<i| u^k) [KO*U\j)) = «i| rfcru \j))* = (j\ rfo^u \i) 

(2.72) 

Computing the equilibrium trace with U'0'*U is thus equivalent to com- 
puting the equilibrium trace in time-reversed states but with O. If we take 
for O the spin <x, the net effect of the time-reversal operation should be to 
change the direction of the spin, in other words, we want 



(2.73) 
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The expression for U will depend on the basis states for spin. Using the 
Pauli matrix basis 

(2.74) 

we have = <x, and a* = a x , a* = —cr y , cr* z = a z so that Eq.(2.73) for 
time reversal gives us the following set of equations for the unitary operator 
U 

U^a x U = -a x (2.75) 
tf(T y U = Uy (2.76) 
U^a z U = -a z (2.77) 
Given the fundamental properties of Pauli matrices 

(JiCfj + Uj<Ji = for i i= j 

-1 = 1 
ia k 



'• u 3 



(2.78) 
(2.79) 

where i,j,k are cyclic permutations of x,y,z, the solution to the set of 
equations for U is 

U = e l5 a y (2.80) 

where S is an arbitrary real phase. This is like a tt rotation along the y 
axis so that already we can expect that up will be transformed into down 
as we were hoping intuitively. In summary, the time reversal operator in 
the presence of spin multiplies the spin part by e lS <j y and takes the complex 
conjugate. 

(2.81) 



T t = Ke iS a 



Note the action of this operator on real spinors quantized along the z direc- 
tion 



Tt |T) 



-le 



14) 



(2.82) 
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T t |4) - ie- lS |T) (2.83) 

The time reversal operator thus transforms up into down and vice versa but 
with a phase. Even if we can choose e~ lS = i to make the phase real, the 
prefactor cannot be +1 for both of the above equations. In particular, note 
that T t T t |t) = — |t) j another strange property of spinors. The application 
of two time reversal operations on spinors is like a 2ir rotation around y so 
that it changes the phase of the spinor. It can be proven that this result is 
independent of the choice of quantization axis, as we can expect. [4] As far 
as the main topic of the present section is concerned, observables such as 
angular momentum will have a simple signature under time reversal (they 
are always two spinors that come in for each observable A{) so that the 
results of the previous section are basically unmodified. 

When Xa a X 1 ^) i s rea ^' the properties of being a commutator (2.85) and of 
Hermiticity (2.87) allow us to further show that x" AiAi (w) is also an odd function 
of frequency, an important result that we show in the following section. 

Properties that follow from the definition. 

Let us thus write down the general symmetry properties of Xa a M ~ that 
simply follow from its definition (2.52). 

• Commutator: Since it is a commutator, we have 

X Ai A j (t-t') = -X AjAi (t'-t) (2.84) 

which in frequency space reads, 



Hermiticity: Taking the observables as Hermitian, as is usually the case, one 
can use the cyclic property of the trace and the Hermiticity of the density 
matrix to show that 



(2.86) 



(Proof for Hermitian operators: {[A i7 Aj])* = Tr {pAiAj — pAjAi}* 

= Tr {AjAip — AiAjp} = Tr {p [Aj,Ai]} with p the density matrix.) 

In Fourier space, this becomes, 



x'LaAu) = x'LaAu) 



(2.87) 



Remark 15 Non-hermitian operators: It is important to note that the operators 
Ai may be non- Hermitian, as is the case for superconductivity. In such cases, one 
should remember that the above property may not be satisfied. 

Remark 16 Most useful property: The most important consequence of this sec- 
tion that we will often use is that correlation functions such as x" p P ( w ) are °dd 
in frequency and real 



V PqP- 



(2.88) 
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To prove this, we first use Hermiticity Eq.(2.87) in the form 

r i * 

to show that x"„ n (lj) is real 



*"PqP- 



i (w) = |d 3 r|dVe-^( r - r ') x '; rPr/ (a;) 



J d s r J d 3 r't 



iq-(r-r') // 



x'LM 



X. p 



(2.89) 

(2.90) 
(2.91) 
(2.92) 



The commutator property Eq.(2.85), Xp qP _ q ( w ) = — Xp_ p q ( — w ) arl ^ symmetry 
under parity transformation Eq.(2.61), x" p q p q ( — w ) = x'p P q (~ w ) ^ en suffice to 
show that Xp p q ( w ) * s a ' so * n frequency x'p nP ( w ) = ~ Xp qP ( — w )- Instead 
of parity, one could have invoked time-reversal symmetry Eq.(2.68) and the com- 
mutator property Eq.(2.85) to show that x'p p ( Ll) ) * s odd since then x'p p , ( w ) = 
Xp r>Pr ( u ) = -Xp rPr ,(-^) immediately implies that Xp qP _ q M = -Xp qP _ q (-w). 



2.3.3 Kramers-Kronig relations and causality 



These Kramers-Kronig relations are by far the best known and most useful rela- 
tions. They relate real and imaginary parts of response functions and they come 
simply from causality. Causality is insured by the presence of the 6 function in the 
expression for the response functions Eq.(2.53). Causality simply states that the 
response to an applied field at time t' occurs only at time t later. This is satisfied 
in general in our formalism, as can be seen by looking back at the formula for 
the linear response Eq.(2.44). Kramers-Kronig relations are the same causality 
statement as above, seen from the perspective of Fourier transforms. To be more 
specific, in this section we will derive the following results: 



Re 



(2.93) 



Im 



X% Aj M 



v r du,' H*W"') 

J 7T UJ'—LO 



(2.94) 



They come from analytic properties of the response functions in the complex 
frequency plane. We give two derivations. 



The straightforward manner: 

Let us first derive the relations the easy way. Suppose that we know the Fourier 
transform in frequency Xa-a-( u ) °f the response function. We call it the retarded 
function because the response comes after the perturbation. It is causal. One way 
to make sure that its real time version X^A,-^ ~ O contains 9{t — t') is to have 
XAiAji^) analytic in the upper half-plane. To see that analyticity in the upper 
half-plane is a sufficient condition to have 9(t — t'), consider 



x% A] (t-t')= / 



2V 



-iw(t-t') R 



x%aM- 



(2.95) 
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If t — t' is negative, then it is possible to close the contour in the upper half 
plane since the exponential will decrease at positive imaginary frequencies. Since 
XAiAj ( w ) i s analytic in that half-plane, the result will be zero, which is just another 
way to say that Xa 4 Aj(* — i s proportional to 9(t — t'), as we had planned to 
show. In the next subsection, we will show that analyticity in the upper half plane 
is also a necessary condition to have 6(t — t'). 

Assuming that XA t Aj ( w ) i s analytic in the upper half plane, it is then easy to 
derive the Kramers-Kronig relations. It now suffices to use 



/ 



1 



7T LU — LU — IT] 



XA i A j {"') = 2iXA i A^ + ir i) 



(2.96) 



which is easy to prove by applying the residue theorem on a contour closed in the 
upper half plane where X^A.,- ( w ) is analytic. This also assumes that X^A,^') 
falls off at least like a small power of so that there is no contribution from 
the part at oo. We then need the following identity, 



lim 



LOzLirj 
lim — ^ ^ 

7J->0 LU Z + T} 1 



lim 



%r\ 



LU 2 + T] 2 ^~ LU 1 + J] 1 



= V-±iir5(w) (2.97) 

LU 



where S is Dirac's delta function and V means principal part integral. — Suppose 
the factor 1/ (oj + iff) on the left is in an integral that can be done by contour 
integration. Then, knowing the definition of the delta function, this can be used 
as the definition of principal part. — Using this identity and setting equal the real 
parts of our contour integral (2.96) we obtain, upon taking the lim 77 — > 0, 



7" 



Re 



Im 



ur — u> 

while from the imaginary part, 



XA iA M 



-2Im 



XA iAj H 



(2.98) 



V 



/v 



Im 



LU' — Cd 



+ Re 



xIaA") 



2 Re 



xIa, (w) 



(2.99) 



This is precisely what we mean by the Kramers-Kronig relations, namely we re- 
cover the results Eqs. (2.93) (2.94) at the beginning of this section. From the proof 
just given, Kramers-Kronig relations will apply if 

• XAiA* ( w ) i s analytic, as a function of complex frequency, in the upper half- 
plane. 



Xa a ( u ) f aus °ff a * least as a small power of u at infinity. 



Spectral representation and alternate derivation. 

It is instructive to perform a derivation which starts from what we found earlier. 
We will gain as a bonus an explicit expression for real and imaginary parts in 
terms of correlation functions, as well as a derivation of the analyticity properties 
from scratch. In fact this will also complete the proof that analyticity in the upper 
half-plane is both necessary and sufficient to have causality. 

Using the convolution theorem, we would write for the frequency-space version 
of the response functions, (2.53) 

Xt4>) = 2z j ^(^^V). (2.100) 
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This looks nice, but it does not really mean anything yet because we encounter a 
serious problem when we try to evaluate the Fourier transform of the 9 function. 
Indeed, 



I 

J — C 



dte wt 6{t) = 



ilo 



DO 





(2.101) 



and we have no idea what e %ulco means. To remedy this, we have to return to 
the expression for the linear response (2.44). Assuming that the external field 
di is turned-on adiabatically from t = — oo, we multiply whatever we had before 
by e vt , taking the limit of vanishing n at the end of the calculation. We also 
adiabatically turn off the response at t — > oo by using a factor e _J?t .The equation 
for the response in time (2.53) is then simply multiplied by e''(* so that it 
still depends only on the time difference. Furthermore, when we take its Fourier 
transform, d(t — t')e lLd ( t ~ t \ everything proceeds as before, except that we 

can use the extra convergence factor e _J) ( t_ * ' , to make sense out of the Fourier 
transform of the Heaviside theta function. To be more specific, the equation for 
the response (2.53) now reads, 



X% Aj {t f)e-"(*-*') = 2i X X Aj (t t')6(t Oe-"(*-° 
so that in the calculation of the response (2.100) we have, 

/oo p i(uj-\-ir/— u/)(t— t') 
d(t - t'y^+^-u'w-^eu - 1') = — : —\f = - 



(2.102) 



1 



i(uj' — u> — in) ' 

(2.103) 

Everything behaves as if we had computed the Fourier transform for oj + ir/ instead 

of LJ, 



duo' X%a,V) 



(2.104) 



(2.105) 

7T Ld' — (LU + 17]) 

This function is called the "retarded response" to distinguish it from what we 
would have obtained with 6{t! — t) instead of 6(t — t'). The retarded response is 
causal, in other words, the response occurs only after the perturbation. In the anti- 
causal case ( "advanced response" ) the response all occurs before the perturbation 
is applied. In the latter case, the convergence factor is e - ^* instead of e 7 ^* 
Introducing a new function 



X Al A, (z) = I 



7T OJ' —Z 



we can write for the retarded response, 



and for the advanced one, 



z—LO— : vr] • 



(2.106) 
(2.107) 
(2.108) 



Using the above results, it is easy to see that XAiAj ( w ) ls analytic in the upper half 
plane, while Xa\a- ( tL) ) i s analytic in the lower-half plane. One can even explicitly 
see from the equation (2.105) for the function which is analytic in the upper-half 
plane x\ iAj ( w ) that the poles or the integrand in the lower-half frequency plane 
are just below the real axis, a distance 77 along the imaginary direction. The 
residue at a given pole will depend on the value of Xa A at the corresponding 
value of the real coordinate of the pole. 
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Remark 17 If x" A . A .{oj') is given by a sum of delta function, then XAiA-i z ) has 
a series of poles just below the real axis. This is what occurs with the so-called 
Lehman representation that we will discuss later. On the other hand, if Xa.^ (<*>') 
is a continuous function, as occurs in an infinite system, then the poles ofxAiAj ( z ) 
are below the real axis, but not ncessarily close to it. The passage from a series 
of poles to a continuous function is what introduces irreversibility in many-body 
systems. 

Definition 4 Equations such as (2.106) are called spectral representations. 

Remark 18 Why spectral representation: The reason for this name is that, as we 
discuss in the next section below, Xa A contains information on dissipation 
or, alternatively, on the spectrum of excitations. Hence, in that kind of equations, 
the response is expressed in terms of the spectrum of excitations. We will also 
have spectral representations for Green's functions. 

XAiAji 2 ) 1S a function which is equal to Xa\a 3 f° r z infinitesimally above 
the real axis, and to Xa^a ( u ) f° r z infinitesimally below the real axis. On the 
real axis of the complex z plane XAiA 3 { z ) has a cut whenever XAiA,( w ) 7^ since 



X A% A 3 ( w + «?) - X Ai A,- ( w - ir l) = 2i XA t A, M- 



(2.109) 



So much for taking the Fourier transform of a response which is so simple 
looking in its ordinary time version Eq.(2.53). Time-reversal invariance (2.68) and 
Hermiticity in Eq.(2.87) imply, for two operators with the same signature under 
time-reversal, that X^A-C^O 1S a rea l function. Hence, from the mathematical 
identity for principal part Eq.(2.97) and from the spectral representation (2.106) 
we have, for two hermitian operators Ai , Aj with the same signature under time 
reversal, that 



Im 



Xa.a.M =x" AiAj { UJ ) 



(2.110) 



so that from the spectral representation we recover the first of the Kramers-Kronig 
relation (2.93). The other one can be derived following the same route as in 



the simpler derivation, namely apply J 



duj 



7T Ul — UJ ! +17] 



on both sides of the spectral 



representation. For two hermitian operators Ai, Aj with opposite signatures under 
time reversal Eqs.(2.68) and (2.87) imply that XAiAji^') ls purely imaginary. In 



this case, 



Re 



X%aM 



1 X%a, 



(2.111) 



Remark 19 Kramers-Kronig and time reversal: The Kramers Kronig relations 
do not depend on these subtleties of signatures under time-reversal. However 
the relation between real and imaginary parts of the response and commutator 
Eq.(2.110) does. If we can compute either the real or imaginary part of the re- 
sponse, the Kramers Kronig relations give us the part we do not know. In any 
case, everything about the system is in XAjA-( w )- 



2.3.4 Positivity of w%"(w) and dissipation 



Without proof, we state that positivity of dissipation implies that 



X%aM)u>0. 



(2.112) 
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We have seen that for Hermitian operators with the same signature under time 
reversal, x'a t A » ( w ) ^ s a rea ^ an d °dd function of frequency so the above equation is 
satisfied. The positive definiteness of x\ i A i by itself however does not suffice 
to prove that Xa A ( w ) i s an °dd function of frequency. 

One can check explicitely that x'a-aX 00 ) contains spectral information about 
excited states by doing backwards the steps that lead us from Fermi's golden rule 
to correlation functions. 



2.3.5 Fluctuation-dissipation theorem 



This very useful theorem relates linear response to equilibrium fluctuations mea- 
sured in scattering experiments. It takes the form, 



S Ai Aj (w) 



2h 



l-e 



(2.113) 



where n B (cj) = 1/ (e' 9 ''" - l) is the Bose factor while the "structure factor" or 
correlation function is defined by, 

S AiAj {t) = (Mt)Aj) (A,) (A,) = ((Mt) (A,)) (A,(0) - (Aj))) (2.114) 

= (SAi(t)SAj) . (2.115) 

We have already encountered the charge structure factor in the context of inelas- 
tic neutron scattering. Clearly, the left-hand side of the fluctuation-dissipation 
theorem Eq.(2.113) is a correlation function for dissipation while the right-hand 
side contains the dissipation function x" just discussed. This is a key theorem of 
statistical physics. 

To prove the theorem, it suffices to trivially relate the definitions, 



then to use the key identity, 



2h 



2h 



s Aj A i (-t) = s Ai A j (t-ih/3) 



(2.116) 



(2.117) 



This kind of periodicity of equilibrium correlation functions will be used over and 
over in the context of Green's functions. It will allow to define Fourier expansions 
in terms of so-called Matsubara frequencies. The proof of the identity simply 
uses the definition of the time evolution operator and the cyclic property of the 
trace. More specifically, using the cyclic property of the trace, time-translation 
invariance follows, and 

S AjAi {-t) = Z~ X Tr [e-^SAji-tyAi] = Z~ x Tr [e~ pH SAjSM*)] ( 2 - 118 ) 

To reverse the order of SAi and SAj, it suffices to use the cyclic property, so that 

S AjAi (-t) = Z- x Tr [{ii(i)e- ffl m . (2.119) 

Simple manipulations and Heisenberg's representation for the time-evolution of 
the operators gives, 



S A M-t) = Z ~ lTr [e- 0H e^ H SA i (t)e^ H SA j 



(2.120) 
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= Z~ x Tr [e-P H 5Ai(t - ihp)SA j ] = S AtAi {t - ih/3). (2.121) 

This is precisely what we wanted to prove. The rest is an exercise in Fourier 
transforms, 

J dte iut S AiAj {t - ihfj) = J dte^+^S^ (t) = e~^S AiAj (w). (2.122) 

To prove the last result, we had to move the integration contour from t to t + ihf3, 
in other words in the imaginary time direction. Because of the convergence factor 
e -P H j n the traces, expectations of any number of operators of the type e lHt Ae~ lHt 
are analytic in the imaginary time direction for — ih(3 < t < ih/3, hence it is 
permissible to displace the integration contour as we did. Fourier transforming 
the relation between X AiAj {t) an d susceptibility Eq.(2.116), one then recovers the 
fluctuation-dissipation theorem (2.113). 
A few remarks before concluding. 

Remark 20 Alternate derivation: Formally, the Fourier transform gives the same 
result as what we found above if we use the exponential representation of the Taylor 
series, 

s AiAj (t-ih/3) = e- ih e&s AiAj (t). 

Remark 21 Relation to detailed balance: The Fourier-space version of the peri- 
odicity condition (2.117) is a statement of detailed balance: 



S AjAi (-u) = e-^S AiAj {u) . (2.123) 



Indeed, in one case the energy Hu is absorbed in the process, while in the other case 
it has the opposite sign (is emitted). This is one way of seeing the basic physical 
reason for the existence of the fluctuation-dissipation theorem: Even though the 
response apparently had two different orders for the operators, the order of the 
operators in thermal equilibrium can be reversed using the cyclic property of the 
trace.. 

Remark 22 Physical explanation of fluctuation- dissipation theorem: Physically, 
the fluctuation- dissipation theorem is a statement that the return to equilibrium is 
governed by the same laws, whether the perturbation was created by an external 
field or by a spontaneous fluctuation. 



2.3.6 Sum rules 



All the many-body Physics of the response or scattering experiments is in the cal- 
culation of unequal-time commutators. These commutators in general involve the 
time evolution of the systems and thus they are non-trivial to evaluate. However, 
equal-time commutators are easy to evaluate in general using the usual commu- 
tation relations. Equal-time corresponds to integral over frequency as seen from 
Fourier space. Hence the name sum rules. We will not in general be able to 
satisfy all possible sum-rules since this would mean basically an exact solution 
to the problem, or computing infinite-order high-frequency expansion. In brief, 
sum-rules are useful to 

• Relate different experiments to each other. 

• Establish high frequency limits of correlation functions. 

• Provide constraints on phenomenological parameters or on approximate the- 
ories. 
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Thermodynamic sum-rules. 

Suppose we compute the linear response to a time-independent perturbation. For 
example, compute the response of the magnetization to a time-independent mag- 
netic field. This should give us the susceptibility. Naturally, we have to leave the 
adiabatic switching-on, i.e. the infinitesimal i]. In general then, 

6(A t (uj = 0)) = x% Aj (w = 0)o> = 0). (2.124) 

Returning to the notation where q is explicitly written, 

<5<^(q,w = 0)) = X %aM^ = 0K(q>^ = 0)- (2.125) 

Using the spectral representation (2.106) and the usual relation between in and 
principal parts, Eq.(2.97), we also have, 

X L 4 ,< q ,, _ o) = r = v r ( 2 .i26) 

J -co K UJ-11] IT UJ 

There is no contribution from the imaginary part on the grounds that there can 
be no dissipation zero- frequency dissipation in a stable system. In fact, as long 
as the thermodynamic derivatives involve operators that have the same symmetry 
under time reversal, then Xa-aXq^) i s °dd, as proven at the end of the section 
on symmetry properties, so that ~x" AiA .{q^,ui = 0) = 0. Note that in practice, 
the principal part in the above equation is not necessary since Xa-a (Qi u ) usually 
vanishes linearly in to for small uj. To be completly general however, it is preferable 
to keep the principal part. 

Recalling that the thermodynamic derivatives are in general for uniform (q = 
0) applied probes, the above formula become, 

BA- 

lim X % Aj (q, w = 0) = ^ = X A iAj ■ (2.127) 



q 



XA i Aj= lim q^oJ_ OC c 



(2.128) 



This is called a thermodynamic sum-rule. As an example, consider the density n 
response. It obeys the so-called compressibility sum rule, 

lim XL(W = 0) = lim r d ^*^l = (»t) . (2 . 129a) 
As usual, a few remarks are in order: 

Remark 23 Order of limits: It is extremely important to note that for thermo- 
dynamic sume rules, the uj — > limit is taken first, before the q — > limit. 
The other limit describes transport properties as we shall see. Another impor- 
tant question is that of the principal part integral. As follows from a problem 
set, in the long wave length limit we have x^„(q,w) = ^^^ Xbb where D is 
the diffusion constant. One can check explicitly with that form that at any fi- 
nite q, it does not matter whether we take or not the principal part integral. We 
did not take it in Eq.(2.129a). If we take the limit q^ before doing the in- 
tegral however, lim q ^o Xnn (l^) takes the form ui5 (uj) so it is important NOT 
to take the principal part integral to get the correct result (in other words, un- 
der the integral sign the 77 — > limit must be taken before the q — > limit). We 
also see this as follows. If we return to the original form lim r) ^ l/( w — iff) = 
lim^o uj/(uj 2 + n 2 ) + in/ (uj 2 + rj 2 ), and then do the integral of the first term (real 
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part), we can check that we have to take the 77 — > limit under the integral sign 
before the q ^ limit to recover the result obtained by doing the integral at finite 
q and then taking the q ^ limit ( the latter is unambiguous and does not depend 
on the presence of the principal part in the integral). Physically, this means that 
the adiabatic turning-on time must be longer than the diffusion time to allow the 
conserved quantity to relax. This is summarized by the following set of equations 



limP r ^x'lm = lim r (2129b) 



q 



= r *± lim ^n(q^) (2 . 129c) 

* V T- Mm X ' LM (2.129d) 
J-00 * w 



Remark 24 Thermodynamic sum-rule and moments: Thermodynamic sum-rules 
are in a sense the inverse first moment over frequency of x' l A i A i (lAO (the latter 
being analogous to the weight). Other sum-rules are over positive moments, as we 
now demonstrate. 

Alternate derivation: Here is another way to derive the thermodynamic sum 
rules. First note that thermodynamic variables involve conserved quantities, 
namely quantities that commute with the Hamiltonian. Take for example N, 
the total number of particles. Since N commutes with the Hamiltonian, in the 
grand-canonical ensemble we have the classical result 

<™>-<*>S(£L 

By definition, 

/OO J 
^S NN {u) (2.129e) 
-00 ^ 

Using the general fluctuation-dissipation theorem, we now relate this quantity to 
Xnn ( w ) as follows. Because n q for q = is simply the total number of particles 
N and hence is conserved, (n q=0 (t) n q=0 ) is time independent. In frequency 
space then, this correlation function is a delta function in frequency. For such a 
conserved quantity, the fluctuation-dissipation theorem Eq.(2.113) then becomes 

oh 9 
S nn (lo) = lim 1 _ e _ phu 7d{ ffi (u) = -^Xnn(") (2.12M) 

from which we obtain what is basically the thermodynamic sum-rule Eq.(2.129a) 

(NN) - (N) 2 = r pS NN (uj) (2.129g) 

1 ( dn \ ( 2 .129h) 



[3lj (3 \dn 



T,V 



This is then the classical form of the fluctuation-dissipation theorem. In this form, 
the density fluctuations are related to the response (dn/d/j,) TV (itself related to 
the compressibility). 

Moments, sum rules, and high-frequency expansions. 

Odd derivatives of x'^A, at equal-time are easy to compute and provide us with 
moments: 



(2.129i) 
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1 

h 



M*)- i,-(0) 



H 



,a,-(o) 



t=0 \ L L L "J "J J / t=0 

(2.129J) 

which may all easily be computed through n equal-time commutations with the 
Hamiltonian. 

Suppose the spectrum of excitations is bounded, as usually happens when 
the input momentum q is finite. Then, x'^A = for > D where D is 
some large frequency. Then, for oj > D, we can expand the denominator since 
the condition cj'/u <C 1 will always be satisfied. This gives us a high-frequency 
expansion, 



2^„=i ^ J_ oc — ^ ) XAiAj <A W ) 



(2.129k) 
(2.1291) 



where we have explicitly taken into account the fact that only odd moments of 
XAiA d° n °t van ish because it is an odd function. Clearly, in the uj — > oo limit, 
the susceptibilities in general scale as 1/w 2 , a property we will use later in the 
context of analytic continuations. 



The f sum-rule as an example. 

When the potential-energy part of the Hamiltonian commutes with the density 
operator, while the kinetic-energy part is that of free electrons (not true for tight- 
binding electrons) we find that 



nq 



(2.129m) 



This is the f sum-rule. It is valid for an arbitrary value of the wave vector q. It is a 
direct consequence of the commutation-relation between momentum and position, 
and has been first discussed in the context of electronic transitions in atoms. The 
proof is as follows. We first use the above results for moments 



dt 



n_ q (t) 



-— ([[#,«,(*)], n_ q (t)]) 



(2.129n) 
(2.129o) 



In the first equality, we have also used translational invariance to write, 



d(r 



iq l r— r 



1 



/(r - r') = - / dre~ lclr / dr' 'e' 1 ^ /(r - r') (2.129p) 



where V is the integration volume. The computation of the equal-time commutator 
is self-explanatory, 



= J dre-^J2^-^) = Y,' 



-«q-r Q 



(2.129q) 



= -hq x e~ i<ir P. 



(2.129r) 



Assuming that the interactions commute with the density operator, and using 
[p • p, n] = p [p, n] + [p, n] p we have 



[H,n q (*)]=E 



'4 

2m 



,n c 



^E(p^-(-^- iqr 



-Rqe"^) ■ Pj8 ) 

(2.129s) 
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[[tf,n q (t)] ,n_ q (t)] = -- ^fi 2 q 2 e- iq - r "e^^ = (2.129t) 

i8=l 

which proves the result (2.129m) when substituted in the expression in terms of 
commutator (2.129o) with n = N/V. The result of the commutators is a number 
not an operator, so the thermodynamic average is trivial in this case! (Things will 
be different with tight-binding models.) 
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3. GREEN FUNCTIONS 



We are now ready to start working with the real many-body problem. This chapter 
will be rather formal but we will make the link with the previous chapter and also 
we will try to do applications as soon as possible in the following chapter. 

The single-particle basis state is a complete basis that is used most often. Note 
however that a simple wave-function such as 

ip(x,y) = {x-y)Ne-\ x - y \' a (3.1) 

for two electrons in one dimension, with N and a constants, is a perfectly ac- 
ceptable antisymmetric wave function. To expand it in a single-particle basis 
state however requires a sum over many (in general an infinite number of) anti- 
symmetrized one-particle states. 

The plan then is as follows. We will then motivate the definition of the Green's 
function in the many-body context by analogy with what we just saw. Then, we 
return to perturbation theory to show that it is most natural to work in imaginary 
time. This leads us to the Matsubara Green's function. We show that if we know 
this Green's function, we also happen to know the retarded one, as well as all 
the one-body quantities of physical interest. We will once more spend some time 
on the interpretation of the spectral weight, develop some formulas for working 
with the Fourier series representation of the imaginary time functions (Matsubara 
frequencies). This should put us in a good position to start doing perturbation 
theory, which is all based on Wick's theorem. Hence, we will spend some time 
proving this theorem as well as the very general linked-cluster theorem that is very 
useful in practice. 



3.0.7 Second quantized operators in the Heisenberg picture 



In a previous chapter, we showed how to translate one- and two-body operators in 
the Schrodinger picture into the language of second quantization. The Heisenberg 
picture is defined as usual. In this section, we derive a few useful identities and 
study the case of quadratic Hamiltonians as an example. 
In the Heisenberg picture 



Ok(t) 



e iHt c k e- lHt 



cl(t) 



JHtJ 



iHt 



(3.2) 



It is easy to compute the time evolution in the case where the Hamiltonian is 
quadratic in creation and annihilation operators. Take for example 



(3.3) 



The time evolution may be found from the Heisenberg equation of motion, which 
follows from differentiating the definition of the Heisenberg operators 



.dCk(t) 



Ck(t),H 



(3.4) 
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To evaluate the commutator, we note that since H commutes with itself, 



Ht u„-iHt 



Ck(t),H = c k {t),e lMt He 



c k (t) , e lHt e k'4' c k' e" 



iHt 



(3.5) 



c k (*), (5Z c k'4' Wck' (*) ) = 2 Ck ' [pk(*),4' (*)Pk' («)] ■ (3-6) 

V k' / J k' 

Commutator identities: The following are very useful identities to get equa- 
tions of motions, and in general equal-time commutators. 

[A, BC] = ABC - BCA = ABC - BAC + BAC - BCA (3.7) 

(3.8) 



[A, BC] = [A,B}C + B [A, C] 



[A, BC] = {A,B}C-B {A, C} 



(3.9) 



The first commutator identity is familiar from elementary quantum mechan- 
ics. The last one can be memorized by noting that it behaves as if the B 
had anticommuted with the A. 

The above identities can then be used to evaluate the needed commutator 
either for fermions 



c k (t) , cl (t) c w (t)] = {c k (t) , cl (t)} c k , (t) + = (5 k , k ,c k (t) 
or for bosons 

c k (t) , 4 (t) Ck/ (t)] = [ck (t) , c k , (t)] c k / (i) + = <5 k , k 'c k (t) 

in either case then, the equation of motion becomes 
,9c k (t) 



dt 



c k (t) , H = e k c k (i) 



whose solution is 
Taking the adjoint, 



c k (t) = e" ieki c k 



4 (*) = cLe* 



(3.10) 
(3.11) 

(3.12) 
(3.13) 
(3.14) 



If we had been working in a basis where H was not diagonal, then repeating 
the steps above, 

. da a (t) 



dt 



a a (t) ,H]=J2(P\H | 7 ) [oa (t) , 4 (t) a 7 (*)] = ^(a\H | 7 ) a 7 (t) 

0,1 7 

(3.15) 

whose solution is found by diagonalizing, integrating, and changing back the basis. 



3.1 Motivation of the definition of the second quan- 
tized Green's function G R 



When the Hamiltonian is quadratic in creation-annihilation operators, in other 
words when we have a one-body problem, the retarded single-particle Green's 
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function we are about to define does reduce to the Green's function we studied in 
the one-body Schrodinger equation. Its actual definition is however better suited 
for many-body problems as we shall see in the present section. 
Consider the definition we had before 

G R (r,t; r', t 1 ) = -i (r| e^*"*') | r ') 9 (t - t') . (3.16) 

Since in second-quantization the operator ip^ (r) creates a particle at point r, the 
following definition seems natural 

G R (r,t; r', t') = -i (GS\ </> (r) e^*"' V (r') \GS) 9 (t - t') (3.17) 

In this expression, \GS) is a many-body vacuum (ground-state). Choosing appro- 
priately the zero of energy, H \GS) = \GS) = so that the above result could 
be written 

G R (r,t; r',t') = -i (GS\ tp (r,t) V f (r', t') \GS) 9 (t - t') . (3.18) 

This is not quite what we want except in the case where there is a single parti- 
cle propagating. Indeed, to keep the physical definition of the propagator, it is 
convenient to have at time t = t' + + 

G R (r,t + 0+; r', t) = -i5 (r - r') (3.19) 

reflecting the fact that the wave-function does not have the time to evolve in an 
infinitesimal time. However, in the present case, the many-body vacuum \GS) is 
a linear combination of Slater determinants, 

\GS) = J dn... J dr N * (n.-.rjv)^ (ri) ...^ (r N )\0) (3.20) 

where ^ (ri...rjv) reduces to the Schrodinger wave function and |0) is a real vac- 
uum. This means that (GS\ ip (r,t) ip^ (r', t) \GS) is not in general a delta function. 
This is a manifestation of the fact that we have a many-body problem and that 
particles are indistinguishable. 

Nevertheless, we can recover the desired simple initial condition Eq.(3.19) even 
in the Many-Body case by adopting the following definition, which in a way takes 
into account the fact that not only electrons, but also holes can now propagate: 

G R (r,i;r',t') = -i(GS\ {ip(r,t),^ (r',i')} \GS)6(t-t') ; for fermions 

(3.21) 

G R (r,i; r', t') = -i (GS\ (r,t) , V f (r', t')] \GS) 9 (t - t') ; for bosons 

(3.22) 

This is the zero-temperature definition. At finite temperature, the ground-state 
expectation value is replaced by a thermodynamic average. Hence we shall in 
general work with 

Definition 5 

(3.23) 



G R (r,t;r',t') = -i/U(r,t),^ (r' ,t')\\ 9 (t - t') : for fe 



G R (r,t;r',t') = -i{ i> (r,t) ,^ {r 1 ,t') )9 (t - t') ; for bosons 



(3.24) 



These definitions have the desired property that at t = t' + + , we have that 
G R (r,t + + ; r', t) = —iS (r — r') as follows from commutation or anti-commutation 
relations 
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Remark 25 Analogies: This definition is now analogous to \ R = 2ix"0(t — tf) 
which we had in linear response. The imaginary part of the Green's function will 
again be a commutator or an anticommutator and hence will obey sum-rules. 

Remark 26 Green's function as a response function: Physically, this definition 
makes obvious that the Green's function is the response to an external probe which 
couples linearly to creation- annihilation operators. In the case of fermions, the 
external probe has to be an anticommuting number (a Grassmann variable, as we 
shall discuss later). 



3.1.1 Examples with quadratic Hamiltonians: 



When the Hamiltonian is quadratic in creation-annihilation operators, the equa- 
tion of motion obeyed by this Green's function is the same as in the one-body 
case. An example of quadratic Hamiltonian is that for free particles 

(r|F|r 1 )=-g(r|r 1 ) = -g*(r-P 1 ). (3.25) 
In the general second quantized case, we write 

H = jdn j dv 2 ^ {v 2 ,t){Y 2 \H\v 1 )4,{Y 1 ,t) (3.26) 

We give two calculations of the Green's function, one directly from the definition 
and one from the equations of motion (Schrodinger's equation). 

Calculation from the definition. For a quadratic Hamiltonian, one can also 
compute directly the Green's function from its definition since, if \n) is an 
eigenbasis, <j> n (r) = (r |n), (n'\ H \n) = E n 5 n , n < 

^ (r,t) = J2 < r l n > a « (*) = £ e ~* Ent < r I") a « = J2 e ~ %Ent( t>n (r) a n 

n n n 

(3.27) 

(r,t) , V f (r', 0) } = J2 E e ~ iEnt( t>n (r) {«„, «L} <t>* m (O = £ e~ iE ^ n (r) C (r') 

n m n 

(3.28) 

G R (r,t; r', 0) = -i ( { </> (r,i) , V f (r', 0) } } (t) = -i £ e" iB » V„ (r) C (r') (t) 



(r, r» = / dfe*^)* (-i) ^ e- <B » V„ (r) C (r') 6> (t) = £ 



(3.29) 

0) C 0') 
w + in - E n 

(3.30) 



Calculation from the equations of motion In general, the equation of mo- 
tion can be obtained as follows 

i^G R (r,i; r', f ) = [-i ({^ (r,t) , ^ (r', f )}) f (* - *')] (3-31) 

{ij (r,t) , V f (r', f)}) * (* " + * ({ V (r,t)] , ^ (r', O}) (* " O 

(3.32) 
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Following the steps analogous to those in Eq.(3.10) above, using the anti- 
commutation relations Eqs. (??)(??) it is clear that 



ff,V(r,i) 



r\H |n) V(ri^) 



so that 



i^G* (r,i;r',0 



(3.33) 



(3.34) 



= 5{v-v')5{t-t')-i J dv 1 {v\H\v 1 ){[^{v 1 fy,^ {v',t')))e{t-t') 

= S(r-r')5(t- t') + J dn (r| H |n) G K (n,t; r', i') (3.35) 

This last expression may be rewritten as 

J dr 1 (r\i^-H\r 1 )G R (r 1 ,t;r',t') = 5 (r - r') 5 (t - t') (3.36) 

= (r \r')S{t-t') (3.37) 

where we recognize the equation (??) found in the previous Chapter. For- 
mally then 

(r| (i^ Hj G R (t-f) \r') = (r |r') S (t t') (3.38) 

so that the operator form of the Green's function is the same as that found 
before, namely 



( " ■ <Ht-t') 



(3.39) 



It is convenient to rewrite the result for the equation of motion Eq.(3.36) in 
the following form that is more symmetrical in space and time. 

J dr, J dh (r| - H |n) 5 (t - h) G R (n,ti; r', t') = 5(r-r')5(t- t 1 ) 

(3.40) 

We may as well let time play a more important role since in the many-body 
case it will be essential, as we have already argued in the context of the 
frequency dependence of the self-energy. The inverse of the Green's function 
in this notation is just like above, 

G R (r,t; r l5 hy 1 = (r| i^-H |n) 6 (t - h) . (3.41) 

Seen from this point of view, the integrals over time and space are the 
continuum generalization of matrix multiplication. The delta function is 
like the identity matrix. 

Definition 6 The following short-hand notation is often used 



G R {1,V) = G R {r,t;r',t') 



G R (1,T) 1 G R (1,1') =5{l-V) 



(3.42) 



(3.43) 



where the index with the overbar stands for an integral. 



MOTIVATION OF THE DEFINITION OF THE SECOND QUANTIZED GREEN'S FUNC- 
TION G R 59 



3.2 Interaction representation and time-ordered prod- 
uct 



Perturbation theory in the many-body case is less trivial than in the one-body 
case. In the many-body case time and frequency dependence are unavoidable. To 
construct perturbation theory we will follow the same steps as those used in the 
derivation of linear response theory. The only difference is that we will write a 
formally exact solution for the evolution operator in the interaction representation 
instead of using only the first order result. The important concept of time-ordered 
product comes out naturally from this exercise. 

The plan is to recall the Heisenberg and Schrodinger pictures, and then to 
introduce the interaction representation in the case where the Hamiltonian can be 
written in the form 

H = H + V (3.44) 

where 

[H o ,V]^0 (3.45) 

Let us begin. We assume that H is time independent. Typical matrix elements 
we want to compute at finite temperature are of the form 

(i|e-^(*)^(f)K> (3-46) 

We do not write explicitly indices other than time to keep the notation simple. 
Recall the Heisenberg and Schrodinger picture 

i, H (t) = e im i> s e- im (3.47) 

We define the time evolution operator 

U (t, 0) = e- lHt (3.48) 

so that 

1> H (t) = U(0,t)1> s U(t,0) (3.49) 

Because from now on we assume time-reversal symmetry, we will always make the 
replacement 

C/ f (t, 0) = U (0, t) (3.50) 
as we just did. The differential equation for the time-evolution operator is 

.dU (t,0) 



dt 



= HU(t,0) (3.51) 



With the initial condition U (0, 0) = 1 it has U (t, 0) = e lHt as its solution. It 
obeys the semi-group property 

U (t, t') = U (t, 0) U (0, t') = e- iH (*"*') (3.52) 

U- 1 (t,0)=U(0,t) (3.53) 
U(to,to) = l (3.54) 

for arbitrary t 
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We are now ready to introduce the interaction representation. In this repre- 
sentation, the fields evolve with the unperturbed Hamiltonian 



a iH t 



-iHgt 



(3.55) 



Note that we now use the caret (hat) to mean "interaction picture". We hope this 
change of notation causes no confusion. To introduce these interaction represen- 
tation fields in a general matrix element, 

(i\ e- m ^ H (t) tf H (t') \i) = (i\ e-? H U (0, t) i> s U (t, 0) U (0, t') 4*7 (t', 0) \i) 

(3.56) 

it suffices to notice that it is easy to remove the extra e lHot coming from the 
replacement of ip s by e~ lHot ip (t) e tHat simply by including them in the definition 
of the evolution operator in the interaction representation 



U(t,0) = e lHot U(t,0) 

U(0,t) = U(0,t)e- lHot 
U (t, 0) U (0, t) = U (0, t) U (t, 0) = 1 



(3.57) 

(3.58) 
(3.59) 



With these definitions, we have that our general matrix element takes the form 

(i\ e-? H ^ H it) ^ H (f) \i) = (i\ e^ H U (0, t) i> (t) U (t, 0) U (0, t') $ (f ) U (f, 0) \i) 

(3.60) 

The purpose of the exercise is evidently to find a perturbation expansion for the 
evolution operator in the interaction representation. It will be built starting from 
its equation of motion 

i 0) = e* Hot {-H + H)U (t, 0) = e iHot V (e~ iHot e iHot ) U (t, 0) (3.61) 

Since a general operator is a product of tp fields, it will also evolve with time in 
the same way so it is natural to define the interaction representation for V as well. 
Our final result for the equation of motion for U (t, 0) is then 

Multiplying on the right by U {0,to) we have a more general equation 



; dU(t,t ) _ 



at 



= V(t)U(t,to) 



(3.62) 



Remark 27 Difficulties associated with the fact that we have non-commuting op- 
erators: The solution of this equation is not e - l S v i t ) dt ^ We will see momentarily 
how the real solution looks formally like an exponential while at the same time 
being very different from it. To write the solution as a simple exponential is wrong 
because it assumes that we can manipulate U (t,to) as if it was a number. In re- 
ality it is an operator so that du ^ to ^ JJ (t, f ) _1 ^ In U (t,t ) . Indeed, note the 
ambiguity in writing the definition of this derivative: Should we write 

^-\nU (t, t ) = lim U (t, to)' 1 \G (t + At, t ) - U (t, t )] /At 



or 



lim 

At^O 



U(t + At, t ) - U (t, to) U (t, to)' 1 /At ? 



(3.63) 
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The two limits cannot be identical since in general 



lim 

At^O 



U(t + At,t ),U(t,t ) *1 ^0. (3.64) 



because U (t, t ) is made up of operators such as V and e lH(,t that do not commute 
with each other. 

To solve the equation for the evolution operator Eq.(3.62), it is more conve- 
nient to write the equivalent integral equation that is then solved by iteration. 
Integration on both sides of the equation and use of the initial condition Eq.(3.54) 
gives immediately 



(3.65) 
(3.66) 



f ^tl^ dt 1 = -i f dt'V (f ) U (t',t ) 

Jt at Jto 

U(t,t ) = l-i [ dt'V (f) U {t', t ) 

Jto 

Solving by iteration, we find 

U(t,t ) = 1 -i [ dt'V (t')U (t' ,t ) = (3.67) 

Jto 

= 1 - i j dt'V (t') + {-if j dt'V (*') j dt"V {t") (3.68) 

J to J to J to 

t t' t" 

+ (-t) 3 f dt'V(t') f dt"V{t") f dt"'V (t'") + ... (3.69) 

J tQ J t() J t() 

Suppose t > to and consider a typical term in this series. By suitably defining a 
contour C and time-ordering operator along this contour T c , it can be rearranged 
as follows 



ft 

H) 3 



f dt'V{t') f dt"V(t") ( dt"'V{t"') (3.70) 

J tQ JtQ JtQ 



= H) 3 ^t c 



/ dh&ih) f dt 2 V{t 2 ) f dt 3 V(t 3 ) 
Jc Jc Jc 



(3.71) 



where 

• C is a contour that is here just a real line segment going from to to t. 

• T c is the "time-ordering operator" . Assuming t > t , it places the operator 
which appear later on the contour C to the left. For the time being, T c orders 
operators that are bosonic in nature. A generalization will appear soon with 
fermionic Green's functions. 

• The integral on the lef-hand side of the last equation covers all possible times 
such that the operators with the time that is largest (latest) are to the left. 
The || comes from the fact that for a general V (ii) V (t 2 ) V (t 3 ) there are 
3! ways of ordering the operators. All these possible orders appear in the 
integrals on the right-hand side of the last equation. The operator T c always 
orders them in the order corresponding to the left-hand side, but this means 
that the integral on the left-hand side appears 3! times on the right-hand 
side, hence the overall factor of ^. 

• A product of operators on which T c acts is called a time-ordered product. 
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One also needs U (0, t). In this case, with t > 0, the operators at the earliest 
time are on the left. This means that the contour on which the T c is defined is 
ordered along the opposite direction. 

A general term of the series may thus be written as 



OO 

^(Mo) = ]TH) fe TTT„ 



kV 



which we can in turn write in the convenient notation 



U(t,t ) = T c 



exp (^-i J c dtiV (tij^j 



(3.72) 



(3.73) 



where the contour is as defined above. In other words, operators are ordered right 
to left from to to t whether t, as a real number, is larger or smaller than t . 

We can check the limiting case [H , V] = 0. Then V is independent of time and 
we recover the expected exponential expression for the time evolution operator. 

The definition of the time-ordering operator is extremely useful in practice not 
only as a formal device that allows the time evolution to still look like an expo- 
nential operator (which is explicitly unitary) but also because in many instances it 
will allow us to treat operators on which it acts as if they were ordinary numbers. 

In the zero-temperature formalism, the analog of U (t, t ) is the so-called 5 
matrix. The time-ordering concept is due to Feynman and Dyson. 

Remark 28 Non-quadratic unperturbed Hamiltonians: It is important to notice 
that in everything above, Ho does not need to be quadratic in creation- annihilation 
operators. With very few exceptions however, [10] it is quadratic since we want 
the "unperturbed" Hamiltonian to be easily solvable. Note that the case where Hq 
is time dependent can also be treated but in this case we would have an evolution 
operator U (t, 0) instead of e~ lHot . The only property of the exponential that we 
really use in the above derivation is the composition law obeyed by time- evolution 
operators in general, namely Uq (t,t') Uo (f ,t") = Uq (t,t") . 

Remark 29 The general case of time- dependent Hamiltonians: The problem we 
just solved for the time evolution in the interaction picture Eq.(3.62) is a much 
more general problem that poses itself whenever the Hamiltonian is time- dependent. 



3.3 Kadanoff-Baym and Keldysh-Schwinger contours 



While we have discussed only the time evolution of the operators in the interaction 
representation, it is clear that we should also take into account the fact that the 
density matrix e~^ H should also be calculated with perturbative methods. The 
results of the previous section can trivially be extended to the density matrix by 
a simple analytic continuation t — > —it. In doing so in the present section, we will 
discover the many advantages of imaginary time for statistical mechanics. 

Let us define evolution operators and the interaction representation for the 
density matrix in basically the same way as before 

e~ m = U (-if3, 0) = e - iff o(-ifl U {-if3, 0) = e- (m «U {-if3, 0) (3.74) 

The solution of the imaginary time evolution equation 

. dU(it",0) r r ( ..„,fj ,..„ n , 
1 a/,-„A = v ( lt > u ( lt >°) 
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is then 



U(-i0,O) = T c 



exp — i 



d{it")V{it") 



where 



Im (t) 



(3.75) 
(3.76) 

V(it") = e- t " H °Ve t " H ° (3.77) 

and the contour C now proceeds from t" = to t" = —f3. 

Overall now, the matrix elements that we need to evaluate can be expressed in 
such a way that the trace will be performed over the unperturbed density matrix. 
Indeed, using our above results, we find 

(i\ e-? H i, H (t) V+ (*') |t) = (i\ e-P H °U{-0, 0) U (0, t) j> (t) U (t, 0) U (0, t') ^ + (f ) U (f, 0) \i) 

(3.78) 

We want to take initial states at a time to so that in practical calculations where 
the system is out of equilibrium we can choose to = — oo where we can assume that 
the system is in equilibrium at this initial time. Hence, we are here considering 
a more general case than we really need but that is not more difficult so let us 
continue. Since we are evaluating a trace, we are free to take 



\i) = 17(0, to) |i (to)) 

then we have 

(i\ e-? H = (i (t )\ U (t , 0) e-f m = (i (t )| ( e -^ 2t ' 



(3.79) 



)(' 



,iH„t 0lD -iHto^ p -PH 



(3.80) 



= (i (t )| e-PHo^Haito^^Hito-rP) = ^ e ~PH fj ^ _ ^ Q) ( 3 . 8 1) 

This allows us to write an arbitrary matrix element entering the thermodynamic 
trace as the evolution along a contour in complex time 

(i\ e-? H i> H (t) (t') \i) = (i (to)\e-? H °U(t - i[3, 0) U (0, t) ^ (t) U (t, 0) U (0, f) ^ (f) [7 (f, 0) |i) 

= (i (t )| e-^C/ (t - i& t ) & (t , t) ^ (t) rj (t, t') $ (t') L> (t', t ) \i (to)) (3.82) 

How would we evaluate the retarded Green's function in practice using this 
approach? Take the case of fermions. It is convenient to define G > (t — t') and 
G< (t - t') by 



G> (t - 1' 


) = -i(^ H it)^ H {t')) 




G<(t-t 


>) = i^ H {t>)i> H {t)) 



(3.83) 
(3.84) 



in such a way that 



G R (t U) = i ({^ H (t) , < (t')})8 (t - t') = [G> (t - t') - G< (t - f )] (t - t') 

(3.85) 

To evaluate G > (t — t') for example, we would expand the evolution operators 
such as U (t! ', to) as a power series in V, each power of V being associated with an 
integral of a time ordered product that would start from to to go to the creation 

operator ip (t') , then go to the destruction operator ip (t) until it returns to to — ij3. 
This contour is illustrated in Fig. (3-1). It is this contour that determines the order 
of the operators, so that even if t' is a larger number than t, as illustrated on the 

right panel of this figure, the operator tp (t) always occur after (t') on the 
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Im(t) 



i 




to 


A + 

V(t') 

7 
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to-iP 


^ a', , 

¥(t) 



Im(t) 
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Re(t) 



to 



t -iP 



¥(t) 



3 A + 

V(t') 



► Re(t) 



Figure 3-1 Kadanoff-Baym contour to compute G > (i — t') 



contour, i.e. tp (t) is on the left of ip (if) in the algebraic expression. The parts 
of the contour that follow the real axis are displaced slightly along the imaginary 
direction for clarity. 

We will see momentarily that it is possible to avoid this complicated con- 
tour to make calculations of equilibrium quantities. However, in non-equilibrium 
situations, such contours are unavoidable. In practice however, what is used 
by most authors is the Keldysh-Schwinger contour that is obtained by inserting 

U (t',oo) U (oo,t') = 1 to the left of (£') in the algebraic expression Eq.(3.82). 
In practice this greatly simplifies the calculations since the contour, illustrated in 
Fig. (3-2), is such that integrals always go from — oo to oo. To specify if a given 
creation or annihilation operator is on the upper or the lower contour, a simple 
2x2 matrix suffices since there are only four possibilities.. 



Im(t) 







to 


A + 

r(f) 


— > 






' 4>'<t) 


to-iP 





* Re(t) 



Figure 3-2 Keldysh-Schwinger contour. 

In equilibrium, the analog of the fluctuation dissipation theorem in the form 
of Eq.(2.117) for correlation functions, allows us to relate G > and G < , which 
means that we can simplify matters greatly and work with a single Green function. 
Fundamentally, this is what allows us to introduce in the next section a simpler 
contour that is extremely more convenient for systems in equilibrium, and hence 
for linear response. 
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3.4 Matsubara Green's function and its relation to 
usual Green's functions. (The case of fermions) 



In thermodynamic equilibrium the time evolution operator as well as the density 
matrix are exponentials of H times a complex number. To evaluate these operators 
perturbatively, one needs to calculate time-ordered products along a contour in 
the complex time domain that is relatively complicated, as we saw in the previous 
section. In the present section, we introduce a Green's function that is itself a time- 
ordered product but along the imaginary time axis only, as illustrated in Fig. (3-3) 
below. This slight generalization of the Green's function is a mathematical device 
that is simple, elegant and extremely convenient since the integration contour is 
now simple. In a sense, we take advantage of the fact that we are free to define 
Green functions as we wish, as long as we connect them to observable quantities 
in the corresponding manner at the end of the calculation. What makes this 
Green function extremely useful for calculations is the fact that in evaluating 
time-ordered products that occur in the perturbation series, a theorem (Wick's 
theorem) tells us that all correlations functions are related to producs of time- 
ordered Green's functions. So we might as well focus on this quantity from the 
start. For thermodynamic quantities, since only equal-time correlation functions 
are needed, it is clear that evaluation in imaginary time or in real time should 
be equivalent since only t = is relevant. More generally, for time-dependent 
correlation functions we will see that in frequency space the analytic continuation 
to the physically relevant object, namely the retarded function, is trivial. Also, 
the same tricks apply not only to Green's functions but also to general response 
functions such as the density-density correlation function. 

After introducing the so-called Matsubara Green's function itself, we will study 
its properties. First, using essentially the same trick as for the fluctuation-dissipation 
theorem for correlation functions, we prove that these functions are antiperiodic 
in imaginary time. This allows us to expand these functions in a Fourier series. 
The spectral representation and the so-called Lehman representation then allow 
us to make a clear connection between the Matsubara Green's function and the 
retarded function through analytic continuation. As usual, the spectral represen- 
tation also allows us to do high-frequency expansions. We give specific examples 
of Matsubara Green's functions for non-interacting particles and show in general 
how to treat their Fourier series expansions, i.e. how to do sums over Matsubara 
frequencies. 



3.4.1 Definition 

The Matsubara Green's function is defined by 



The definition of Ref.([ll]) has an overall minus sign difference with the definition 
given here. 

Definition 7 The last equation above defines the time ordering operator for fermi- 
ons. It is very important to notice the minus sign associated with interchanging 




(3.86) 




(3.87) 
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two fermion operators. This time-ordering operator is thus a slight generalization 
of the time- ordering operator we encountered before. There was no minus sign in 
this case associated with the interchange of operators. The time-ordering operator 
for bosonic quantities, such as V that appeared in the perturbation expansion, will 
never have a minus sign associated with the exchange of bosonic operators. 

We still need to specify a few things. First, the thermodynamic average is in 
the grand-canonical ensemble 



_ Tr [ e -P(H-^)(D] 

(°} = Tr [ e -/3(Jf- M A0] 



(3.88) 



with \i the chemical potential and N is the total number of particle operator, while 
the time evolution of the operators is defined by 



i> (r,r) = 


e T{H-»N)^ s ( r ) e -r(H-^N) 




■0 f (r,r) = 


e r(H-^N )jfj + (r) e -r(H-KN) 



For convenience, it is useful to define 



K = H-nN 



(3.89) 
(3.90) 

(3.91) 



Several points should attract our attention: 



The correspondence with the real time evolution operators e lHt is done by 
noting that 

T = -lm(t) (3.92) 

or, in general for complex time 



t = it 



Strictly speaking, we should use tp (r, —ir) if we want the symbol tp (r,t) for 
t complex to mean the same thing as before. That is why several authors 
write tp (r,r) for the Matsubara field operator. We will stick with ip (r,r) 
since this lack of rigor does not usually lead to confusion. We have already 
given enough different meanings to^in previous sections! Furthermore, this 
type of change of "confusion" in the notation is very common in Physics. 
For example, we should never write / (k) to denote the Fourier transform of 
/(r). 

ip^ (r,r) is not the adjoint of ip(r,r). However, its analytic continuation 
t — > it is the adjoint of tp (r,t). 

Using as usual the cyclic property of the trace, it is clear that Q depends 
only on t — r' and not on t or t' separately. 

It suffices to define the Matsubara Green's function Q (r, r'; r) in the interval 
—0 < t < P- We do not need it outside of this interval. The perturbation 



expansion of U (— if3, 0) = T c 



exp 



J c drV 



(r)) 



evidently necessitates 



that we study at least the interval < t < j3 but the other part of the 
interval, namely — p3 < t < is also necessary if we want the time ordering 
operator to lead to both of the possible orders of tp and tp^: namely tp^ to 
the left of tp and tp^ to the right of tp. Both possibilities appear in G R . If 
we had only r > 0, only one possibility would appear in the Matsubara 
Green's function. We will see however in the next section that, in practice, 
antiperiodicity allows us to trivially take into account what happens in the 
interval —f3 < r < if we know what happens in the interval < t < f3. 
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The last contour considered in the previous section for U (—ift, 0) = T c exp J c drV (r) 

tells us that the time-ordering operator T T orders along the contour (Im (t) — — ft) > 
(Im(t') = ft) which corresponds to (r = ft) > (r' = —ft). The present con- 
tour is illustrated in Fig. (3-3). 



Im(t) 
A 



V(0) 



V|/(T) 



t 



T = -P 



Re(t) 



x = P 



Figure 3-3 Contour for time ordering in imaginary time. Only the time difference is 
important. The contour is translated slightly along the real-time axis for clarity. 



Remark 30 Role of extra chemical potential in time evolution: The extra chem- 
ical potential in the evolution operator e r ^~i iN ) is convenient to make all oper- 
ators, including the density matrix, evolve in the same way. It corresponds to 
measuring energies with respect to the chemical potential as we will see with the 
Lehman representation below. The extra e~ T ^ N disappears for equal-time quanti- 
ties (thermodynamics) and in the calculation of expectation values (0 + (t)0 (£')) 
for operators O which are bilinear in fermions of the form (ip + ?p) at equal time. 
Indeed in that case one has 0+(t) = e iHt O+e- lHt = e i{H-»N)tQ+ e -i{H-»N)t _ 
When Wick's theorem is used to compute expectation values, the creation and an- 
nihilation operators evolve then as above. In any case, as we just said, the addition 
of the chemical potential in the evolution operator just amounts to measuring the 
single-particle energies with respect to the chemical potential. 



3.4.2 Antiperiodicity and Fourier expansion (Matsubara frequencies) 



Suppose t < 0. Then 



G{t,t'-t) = {^ + (r',0)^(r,r)) 



(3.93) 



Using the cyclic property of the trace twice, as in the demonstration of the 
fluctuation-dissipation theorem it is easy to show that 



g(r,r';T) = -g(r,r'; T + ft) ; r<0 



(3.94) 
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This boundary condition is sometimes known as the Kubo-Martin-Schwinger (KMS) 
boundary condition. 



Proof: Let 

then 



Tr [e-? K ] 



G (r, r'; r) = e fm Tr [ e -^V + (O {e Kr ^ (r) e~ KT )] 
The cyclic property of the trace then tells us that 

Q (r, r'; r) = e^ 2 Tr [(e Kr <0 (r) e- KT ) e"' 3 * V + (r')] 



s^Tr[(e 



-0KRK\ 



(e KT i> (r) e~ KT ) e"^V + (r')] 



(^(r,r + /3)^+ (r',0)) 
= -S(r,r';r + /3) 



(3.95) 
(3.96) 

(3.97) 
(3.98) 
(3.99) 
(3.100) 



The last line follows because given that — (3 < t, we necessarily have t+/3 > 
so that the other function must be used in the definition of the Matsubara 
Green's function. 



If r > 0, the above arguments can be repeated to yie 



e(r,r';T-/3) = -S(r,r';r) ; t>0 



However, for r > note that 

g(r,r';T)^-g(r,r';r + p) 



T > 



(3.101) 



(3.102) 



While g (r, r'; r + (3) for r > is well defined, we never need this function. So we 
restrict ourselves to the interval — (3 < r < (3 described in the previous section. 

One can take advantage of the antiperiodicity property of the Green's function 
in the interval — /3 < r < (3 to expand it in a Fourier series that will automatically 
guaranty that the crucial antiperiodicity property is satisfied. More specifically, 
we write 

(3.103) 



g(r,r';r) = ^n=-oo^ lKT G(r,r';ik n ) 



where the so-called Matsubara frequencies for fermions are odd, namely 

n integer 



k n = (2n + l)7rT 



(2k+1)tt 





(3.104) 



The antiperiodicity property will be automatically fulfilled because e %knl3 — 

e -i(2n+l)7T _ 

The expansion coefficients are obtained as usual for Fourier series of antiperi- 
odic functions from 



g(ry-ik n ) = ^dTe ik ^g{vy-r) 



(3.105) 



Note that only the r > region of the domain of definition is needed, as promised. 

Remark 31 Domain of definition of the Matsubara Green's function: The value 
ofg(r,r';r) given by the Fourier series (3.103) for r outside the interval —(3 < 
t < (3, is in general different from the actual value of Eq.(3.86) g (r, r'; r — r') = 
— (T T ip (r,T) ^ + (r',r')). Indeed, to define a Fourier series one extends the func- 
tion defined in the interval —(3 < r < (3 so that it is periodic in t outside this inter- 



val with a period 2/3. The true function g (r, r'; t 



(T r ^(r,T)V; + (r',r')) 



has an envelope that is, instead, exponential outside the original interval. We 
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will see an explicit example in the case of the free particles. In perturbation ex- 
pansions, we never need Q{y,y';t) outside the interval where the series and the 
true definition give different answers. To avoid mathematical inconsistencies, it is 
nevertheless preferable in calculations to do Matsubara frequency sums before any 
other integral! It is possible to invert the order of integration and of summation 
but we must beware. 



3.4.3 Spectral representation, relation between G R and Q and analytic continuation 

By analogy with what we have done previously for response functions Xi it is 
useful to introduce the spectral representation for the retarded Green's function. 
We obtain explicitly Q (r, r';ik n ) by integration in the complex plane and find that 
is trivially related to G R (r, r';ui) . 
As before, we have 

G R (r, r';t) = -i ( { </> (r,t) (r', 0) } ) 9 (i) (3. 106) 

but this time, the evolution operator is defined to take into account the fact that 
we will work in the grand-canonical ensemble. By analogy with the definition of 
the Matsubara operators, we now have 

K = H-/jlN 

tP (r ,t) = e ltK iP s (r) e~ itK (3. 107) 

V> f (r,t) = e itK ^ (r) e~ ltK (3. 108) 

We now proceed by analogy with the response functions. On the left we show 
the definitions for response functions, and on the right the analogous definitions 
for response functions. Let 

G R {v,v' ] t)^-iA{v,v';t)9{t) ; X fj («) = 2i*& (*) 9 (t) (3.109) 

where the spectral weight is defined by 



A(r,r';^<{^)^ + (r',0)}> . (t) = ([A (r ,t) , A, (r' , 0)}) (3.110) 



Then taking the Fourier transform, one obtains the spectral representation 



G" (r,r;w) = J 



-oo 2-7T uj+irj—uj' 



X 



13 J -co I" u>' + irf) ' 

(3.111) 

The spectral weight will obey sum-rules, like x" did. For example 

(3.112) 



Ho A (r, r>') = ({ ^ (r,0) , ^ (r', 0)}> = 5 (r - r') 



From such sum rules, a high-frequency expansion can easily be found as usual. 
But that is not our subject for now. 

To establish the relation between the Matsubara Green's function and the 
retarded one, and by the same token establish the spectral representation for Q, 
consider 



Q (r, r'; r) = - (i/> (r,r) V+ (r', 0)) 9 (r) + ty+ (r', 0) $ (r,r)> 9 (-r) (3.113) 
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g(r,r';ik n ) = / dre lk ™ T Q (r, r': 
Jo 



r) 



dre 



ik n r 



(^(r,r)V + ( r ',0))] 



(3.114) 



(3.115) 



Assume that k n > 0. Then, as illustrated in Fig. (3-4), we can deform the contour 
of integration within the domain of analyticity along Re (t) — Im (r) > 0. (The 
analyticity of (if) (r,r) ■0 + (r',0)) in that domain comes from e~@ H in the trace. 
You will be able to prove this later by calculating Q (r, r'; r) with the help of the 
spectral representation Eq.(3.122) and tricks for evaluating sums on Matsubara 
frequencies). For Im (r) = oo there will be no contribution from the small segment 
since e lknT becomes a decaying exponential. The integral becomes 



Im(t) = 



Re(x) 



Re(x) = -p ^ 



Re(x) = p 



->> > — 

x = it 
— < <- 



f-Re(t) = Im(x) 



00 



Figure 3-4 Deformed contour used to relate the Matsubara and the retarded Green's 
functions. 



Q (r,r';ik n ) = 

pt — OO 

/ d(it)[-(e^ s (r)e-^ + s (r'))]e^ 
Jt=o 



(3.116) 



t=o 



d(it) 



] ip s (r) e 



-iK(t-i/3)„,,+ (r / 



D (ifc„)i(t— if. 



In the last integral, we then use the results 

e (ik n )i(-ip) = e (ik n )p = _2 



e? K e iKt 1> 3 (r) e^'e^^s (r') 



oo JO 



(3.117) 
(3.118) 



,iK{t-if))^ s (r) e -iK{t-iff)^ (r') 



Kt e -dK^ fr." 



(3.119) 

It then suffices to cancel the left most e@ K with the density matrix and to use the 
cyclic property of the trace to obtain for the integrand of the last integral, 

(3.120) 



(3.121) 



= -(V t (r',0)V(r,t) / 
Overall then, the integral in Eq.(3.116) is equal to 

G (r, r'; ik n ) = -i £° dt ({i/t (r,t) , V f (r', 0) } ) e^>< 



MATSUBARA GREEN'S FUNCTION AND ITS RELATION TO USUAL GREEN'S FUNC- 
TIONS. (THE CASE OF FERMIONS) 71 



c( r r'-ik \ - r°° do/ A ( r y->"') 

y (r, r , iK n ) - j_ oo ik -—, 



(3.122) 



All that we assumed to deform the contour was that k n > 0. Thus, ik n — > 
uj + in with 77 > is consistent with the hypothesis and allows us to deform 
the contour as advertized. Comparing the formula for Q(r,r';ik n ) for k n > 
with the expression for the retarded Green's function(3.111), we see that analytic 
continuation is possible. 



G R (r, r'; uj) = lim ifej ^ +it? Q (r, r'; ik n ) 



(3.123) 



If we had started with k n <C 0, analytic continuation ik n 
Green's function would have been possible. 



w — in to the advanced 



Remark 32 Connectedness and periodicity: For a general correlation function, 
similar spectral representations can also be defined for connected functions (see 
below) so that periodicity or anti-periodicity can be defined. 



3.4.4 Spectral weight and rules for analytical continuation 



In this section, we summarize what we have learned for the analytic properties of 
the Matsubara Green's function and we clarify the rules for analytic continuation. [12] 

The key result for understanding the analytical properties of Q is the spectral 
representation Eq.(3.122) 



G(r,r 



i ikn) = / 

J — I 



doo' A(r,r';uj') 
2tt ik n — uj' 



(3.124) 



The spectral weight A (r, r'; u') was discussed just in the previous subsection (See 
also Eq. (3.179)). 

The Matsubara Green's function and the retarded functions are special case of 
a more general function defined in the complex frequency plane by 



G(r r'-z) = f°° d^ A{r ' r '^') 



(3.125) 



This function is analytic everywhere except on the real axis. Physically interesting 
special cases are 

G (r,r';ik n ) = G (r,r';ik n ) 
G R (r, r'; lo) = lim G (r, r'; u + irf) 



G A (r, r'; uj) = lim G (r, r'; uj - in) 

j]— >o 

The function G (r, r'; z) has a jump on the real axis given by 



A (r, r'; u') = i lim^^ [G (r, r'; uj + in) — G (r, r'; uj — ir/)} 



(3.126) 
(3.127) 

(3.128) 



A (r, r'; uj') = i [G R (r, r'; u) - G A (r, r': 



co- 



in the special case where A(r,r';ui') is real (which is almost always the case in 
practice since we consider r = r' or k = k'), we have 



A(r,r';u') = —2 Im G R (r,r';w) 



(3.129) 
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like we have often used in the one-body case. 

The previous results are summarized in Fig. (3-5) which displays the analytic 
structure of G (r, r'; z) . This function is analytical everywhere except on the real 
axis where it has a branch cut leading to a jump Eq.(3.128) in the value of the 
function as we approach the real axis from either the upper or lower complex half- 
plane. The limit as we come from the upper half-plane is equal to G R (r, r'; co) 
whereas from the lower half-plane it is equal to G A (r, r'; uj) . The Matsubara 
Green's function is defined only on a discrete but infinite set of points along the 
imaginary frequency axis. 



G(z) = gV) 



E*- Re(z) 
G(z) = gV) 

Figure 3-5 Analytical structure of G(z) in the complex frequency plane. G{z) 
reduces to either G R (uj) , G A (uj) or Q (iuj n ) depending on the value of the complex 
frequency z. There is a branch cut along the real axis. 

The problem of finding G R (r,r'; uj) along the real-time axis from the knowl- 
edge of the Matsubara Green's function is a problem of analytical continuation. 
Unfortunately, G (z = ik n ) does not have a unique analytical continuation be- 
cause there is an infinite number of analytical functions that have the same value 
along this discrete set of points. For example, suppose we know G (z = ik n ) , then 
GO) (1 + {eP z + 1)) has the same value as G (z) for all points z = ik n because 
e ifc„/3 -|-i — o. Baym and Mermin[13], using results from the theory of complex 
functions, have obtained the following result. 

Theorem 8 // 

1. G (z) is analytical in the upper half-plane 

2. G (z) = Q (ik n ) for all Matsubara frequencies 

3. lim z _ >00 zG (z) = est 

then the analytical continuation is unique and 

G R (r,r';w) = lim Q(r,r';ik n ) (3.130) 

The key point is the third one on the asymptotic behavior at high frequency. 
That this is the correct asymptotic behavior at high frequency follows trivially from 
the spectral representation Eq.(3.125) as long as we remember that the spectral 
weight is bounded in frequency. The non-trivial statement is that this asymptotic 
behavior suffices to make the analytical continuation unique. In practice this rarely 
poses a problem. The simple replacement ik n — > uj +ir) suffices. Nevertheless, the 
asymptotic behavior reflects a very fundamental property of the physical system, 
namely the anticommutation relations! It is thus crucial to check that it is satisfied. 
More on the meaning of the asymptotic behavior in subsection (3.6.1). 
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3.4.5 Matsubara Green's function in momentum space and non-interacting case 



We first present the definition of the Matsubara Green's function in momentum 
space since this is where, in translationally invariant systems, it will be diagonal. 
With our definition of momentum and real space second quantized operators, and 
our normalization for momentum eigenstates Eq.(??) we have 



g (r, r'; t-t') = -( T t ^ (r,r) ^ (r' ; r')\ - 



= 1 ik.r-ik'-r' = 1 i(k-k').(^) +i (^).(r-r' 



(r |k) (k' |r' ) = -e~ - - - = - e 



T T ^(r|k)c k (r)^4, (r')(k'|r') 

k k' / 

(3.131) 
(3.132) 



Assuming space translation invariance, we can integrate over the center of mass 
coordinate y J d ( iL 2~ il ) = 1- Since 

1 



,i(k-k')-(^ 



V 



(27r) d 5(k-k') =5 K v 



we are left with 

g{v,v'-T-r') = 



V^E^(r)cL(r')e ik '-(-')^ 



(3.133) 



(3.134) 



g(k;r-r') = j d (r - r') e^-') - (t t ± £ c k , (r) <£, (/) e^-') 

(3.135) 

5 (k; t - r') = - (T T c k (r) 4 (r')\ (3.136) 



which could have been guessed from the start! Our definitions of Fourier trans- 
forms just make this work. 

Remark 33 Momentum indices and translational invariance: Note that the con- 
servation of total momentum corresponding to translational invariance corresponds 
to the sum of the momentum indices of the creation- annihilation operators being 
equal to zero. The sign of momentum is counted as negative when it appears on a 
creation operator. 

Example of no n- interacting particles 

For non-interacting particles let us consider a quadratic diagonal Hamiltonian 

K =J2 (e k - /x) c+c k = Ck4 c k (3.137) 

k k 

The result for the Green's function may be obtained either directly from the 
definition or by integrating the equations of motion. Both ways of obtaining 
the simple result 



£ (k;i/c„) 



are instructive, so let us do both. Assuming for one moment that the above result 
is correct, our rules for analytic continuation then immediately give us the retarded 
function 



G R (k;w) 



that has precisely the form we expect from our experience with the one-body case. 
The only difference with the one-body case is in the presence of the chemical 
potential in £ k . 
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From the definition To evaluate the Green's function from its definition, we 
need c k (r) . That quantity may be obtained by solving the Heisenberg equations 
of motion, 

^ = [#o,c k ] = -CkCk (3.140) 

The anticommutator was easy to evaluate using our standard trick Eq.(3.9). The 
resulting differential equation is easy to integrate given the initial condition on 
Heisenberg operators. We obtain, 

c k (r)=e-^ T c k (3.141) 

so that substituting in the definition, 

Go (k; r) = - (T T c k (r) c+) = -e"^ [(c k c+) 6 (r) - (c+c k ) 6 (-r)] (3.142) 

using the standard result from elementary statistical mechanics, 

<^c k ) = /(C k ) = ^ TT (3.143) 
and (c k cj) = 1 — (c k c k ) we obtain 



Go (k; r) = -e~^ [(1 - / (C k )) B (r) - / (C k ) (-r)] . (3.144) 



Remark 34 Inadequacy of Matsubara representation outside the domain of defi- 
nition: We see here clearly that if r < f/ie equality 

Go (k; r + 0) = -g (k; r) (3.145) 

is satisfied because e~^ k/3 (1 - / (£ k )) = / (£ k ) . On i/ie oi/ier /land, 

So (k; r + 3/3) ^ g (k; r + /3) (3. 146) 

as we might aaue believed if we had trusted the expansion 



oo 

G (k;T) = - e-^ T Go(k;ik n ) 



outside its domain of validity! The conclusion is that as long as the Matsubara 
frequency representation is used to compute functions inside the domain —(3 < t < 
f3, it is correct. The perturbation expansion of the interaction picture does not force 
us to use Green's functions outside this domain, so the Matsubara representation 
is safe! 

Remark 35 Alternate evaluation of time evolution: We could have obtained the 
time evolution also by using the identity 

e A C'e A = C + [A, C] + 1 [A, [A, C]] + i [A, [A, [A, C}}} + ... (3.147) 

that follows from expanding the exponential operators. This is less direct. 

Remark 36 Appearance of Go (k; r) : It is instructive to plot Go (k; r) as a func- 
tion of imaginary time. In some energy units, let us take /? = 5, and then consider 
three possible values of Ck- First £ k = 0.2, i.e. for a value of momentum above 
the Fermi surface, then a value right at the Fermi surface, ( k = and finally 
a value £ k = —0.2 corresponding to a momentum right below the Fermi surface. 
These cases are illustrated respectively in Figs. (3-6) to (3-8). Note that the jump 
at t = is always unity, reflecting the anticommutation relations. What is meant 
by antiperiodicity also becomes clear. The extremal values near ±/3 and ±0 are 
simply related to the occupation number, independently of interactions. 
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Figure 3-6 Go (P: T ) f° r a value of momentum above the Fermi surface. 
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Figure 3-7 Qo (p, r) for a value of momentum at the Fermi surface. 




Figure 3-8 Go (p, t) for a value of momentum below the Fermi surface. 
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Let us continue with the derivation of the Matsubara frequency result Go (k; ik n ). 
'u Jo 



Go 



(k; ik n ) = f dre^Go (k; r) = - (1 - / (( k )) f dre^e^ (3.148) 
Jo Jo 



= -(i-/(Ck)) 
(i-/(C k )) 



ifc„ - C k 

e -Ck/3 _ 1 1 



Ck Ck 



The last equality follows because 



e C k /3 1 

(1_/(Ck)) = ^Tl = e-C^ + l 
We thus have our final result Eq.(3.138) for non- interacting particles. 



(3.149) 
(3.150) 

(3.151) 



Prom the equations of motion In complete analogy with the derivation in 
subsection (3.1.1) we can obtain the equations of motion in the quadratic case. 



-S(r)({c k (r),ct})-/T T (^c k ( 



(3.152) 



(3.153) 



Using the equal-time anticommutation relations as well as the Heisenberg equa- 
tions of motion for free particles Eq.(3.140) the above equation becomes, 



|-So(k;T) = -5(T) + C k (T TCk (r)4 



so that the equation of motion for the Matsubara propagator is 



(&+Ck)0o(k;T) = -<J(r) 



(3.154) 



(3.155) 



To obtain the Matsubara-frequency result, we only need to integrate on both sides 
using the general expression to obtain Fourier coefficients Eq.(3.105) 



/' 

Jo- 



d_ 



+ C k So (k;r) 



„ik n r 



dr = -1 



so that integrating by parts, 



Go (k; T)|g_ - ik n G (k; ik n ) + ( k Go (k; ik n ) = -1 



(3.156) 



(3.157) 



Note that we had to specify that the domain of integration includes 0. The inte- 
grated term disappears because of the KMS boundary conditions (antiperiodicity) 
Eq.(3.94). Indeed, antiperiodicity implies that 



e**-' Go (k;r)|g_ = -Go (k;/T) - Go (k;0") = 



(3.158) 



Eq.(3.157) for the Matsubara Green's function then immediately gives us the de- 
sired result Eq. (3.138). 
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3.4.6 Sums over Matsubara frequencies 



In doing practical calculations, we will have to become familiar with sums over 
Matsubara frequencies. When we have products of Green's functions, we will use 
partial fractions in such a way that we will basically always have to evaluate sums 
such as 

T ?*^ ,3 ' 159) 

where T = . We have however to be careful since the result of this sum is 
ambiguous. Indeed, returning back to the motivation for these sums, recall that 



g(k-r)=Tj2 



(3.160) 



We already know that the Green's function has a jump at t = 0. In other words, 



lim G (k;r) 



(ok<£) 



7^ 



lim Q (k;r) = (c£c k ) 



This inequality in turn means that 

-ik n 0~ 



-ifc„cr 



1 



ik n Cit 



(3.161) 



(3.162) 



The sum does not converge uniformly in the interval including t = because the 
1/n decrease for n — > oo is too slow. Even if we can obtain a finite limit for the 
last sum by combining positive and negative Matsubara frequencies, what makes 
physical sense is only one or the other of the two limits r — > O*. 

Remark 37 Remark 38 The jump, lim r ^ - g (k;r) — lim T ^ + G (k;r) is always 
equal to unity because of the anticommutation relations. The slow convergence in 
1/ikn is thus a reflection of the anticommutation relations and will remain true 
even in the interacting case. If the (ik n ) 1 has a coefficient different from unity, 
the spectral weight is not normalized and the jump is not unity. This will be 
discussed shortly. 

Let us evaluate the Matsubara frequency sums. Considering again the case of 
fermions we will show as special cases that 



7i ik n £k 



= /(Ck) = 0o(k;O-) 



TE 



= -l + /(Ck) = So(k;0+) 



(3.163) 



(3.164) 



Obviously, the non-interacting Green's function has the correct jump Q (k;0 ) — 
Go (k;0+) = 1 

Proof: [14] To perform the sum over Matsubara frequencies, the standard trick is 
to go to the complex plane. The following function 



-(3 



1 



eP z + 1 



has poles for z equal to any fermionic Matsubara frequency: z = ik. 
residue at these poles is unity since for 



(3.165) 
Its 



z = ik n + 5z 



(3.166) 
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we have 



-/?- 



P z + 1 



= -/3 



= -/3- 



1 



e ik n /3+f35z i J — le' 3 * 2 + 1 

1 



lim <5z 



= 1 



+ 1 

Similarly the following function has the same poles and residues: 

1 



lim Sz 

z—ik n — >0 



= 1 



(3.167) 
(3.168) 

(3.169) 



To evaluate the r < case by contour integration, we use Cauchy's theorem on 
the contour C\, which is a sum of circles going counterclockwise around the 
points where z is equal to the Matsubara frequencies. Using Eq.(3.168) this 
allows us to establish the equality 



\ - e 



— ik^ 



1 



P ik n ^ k 



1 



clz 



2ni J Ci + lz-C, v 



(3.170) 



This contour can be deformed, as illustrated in Fig. (3-9), into C2+C3 (going 
through C[) with no contribution from the semi-circles at Re (z) — ±00 
because e ft} +1 insures convergence when Re (z) > despite e~ ZT in the 
numerator, and e~ ZT insures convergence when Re (z) < 0, r < 0. With 
the deformed contour C2 + C3, only the contribution from the pole in the 
clockwise direction is left so that we have 



/3 Z-m ifc„ — Ck 



e -Ck T 



-C k T 



/(C k ) 



(3.171) 



which agrees with the value of Go (k; r) in Eq.(3.144) when r < 0. In partic- 
ular, when r = 0~ we have proven the identity (3.164) . To evaluate the 

Im(z) 



C 



2 : 



C' 



1 



■ X:^i > Re(z) 



ft 

* 



Figure 3-9 Evaluation of fermionic Matsubara frequency sums in the complex plane. 

r > case we use the same contour but with the other form of auxiliary 
function Eq.(3.169). We then obtain, 



dz 



1 e- ik " T 1 



(3.172) 
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This contour can be deformed into C2 + C3 with no contribution from the 
semi-circles at Re (z) = ±00 because this time e~ ZT insures convergence 
when Re (z) > 0, t > and e - ( L +1 ensures convergence when Re (z) < 
despite e~ zr in the numerator. Again, from C2 + C3, only the contribution 
from the pole in the clockwise direction survives so that we have 



(3.173) 



which agrees with the value of Q (k; r) in Eq.(3.144) when t < 0. In partic- 
ular, when r = + we have proven the identity (3.163). 



3.5 Physical meaning of the spectral weight: Qua- 
si particles, effective mass, wave function renor- 
malization, momentum distribution. 



To discuss the Physical meaning of the spectral weight, we first find it in the 
non-interacting case, then write a formal general expression, the Lehman repre- 
sentation, that allows us to see its more general meaning. After our discussion 
of a photoemission experiment, we will be in a good position to understand the 
concepts of quasiparticles, wave-function renormalization, effective mass and mo- 
mentum distribution. We will even have a first look at Fermi liquid theory, and 
see how it helps us to understand photoemission experiments. 



3.5.1 Spectral weight for non-interacting particles 



The general result for the spectral weight in terms of the Green's function Eq.(3.128) 
gives us for non-interacting particles 



i4o(k,w) = 



1 



1 



2nS (to - C k ) 



(3.174) 
(3.175) 



In physical terms, this tells us that for non-interacting particles in a translationally 
invariant system, a single excited particle or hole of momentum k added to an 
eigenstate is an true excited eigenstate located an energy uj = £ k above or below 
the Fermi level. In the interacting case, the Lehman representation will show us 
clearly that what we just said is the correct interpretation 



3.5.2 Lehman representation 



For a general correlation function, not necessarily a Green's function, one estab- 
lishes the connection between Matsubara functions and retarded functions by using 
the Lehman representation. This representation is also extremely useful to extract 



80 



GREEN FUNCTIONS 



the physical significance of the poles of correlation functions so this is why we in- 
troduce it at this point. We have already seen examples of Lehman representation 
in the one-body case when we wrote in Eq.(3.30), 



G R (r,r» = £ 



<Ur)C(rQ 
to + iri - E n 



Let us consider the more general many-body case, starting from the Matsubara 
Green's function. It suffices to insert a complete set of energy eigenstates between 
each field operator in the expression for the spectral weight 

A(r,r';t) = ({if, (r,t) , ^ (r', 0)}) (3.176) 



[{n\ e-P K e iKt i> s (r) e~ iKt \m) (m| 4 (r') \n) 



m.n 



+ (n\ e~ pK ^ s (r') |m) (m\ e lKt ^J s (r) e~ iKt \n) 

We now use e~ %Kt \n) = e~ lKnt \n) with K n = E n — [iN n if there are N n par- 
ticles in the initial state |n) . In the first term above, (n| has one less parti- 
cle than | to) while the reverse is true in the second term so that K m — K n = 
(E m - (j, (N n + 1) - E n + nN n ) in the first term and ^ n --FT ro = (E n - fiN n -E m + (j, (N n - 1)) 
in the second. Taking the Fourier transform J dte lul * we have 

A(r,r';a/) = e 0Q x (3.177) 
[e- pK " (n\ iJ S (r) \m) (m\ ^ s (r') \n) 2n5 (a/ - (E m - M - E n )) 

mn 

+ e"^» (n| ^1 (r') \m) (m\ ^ s (r) \n) 2n5 (w' -(E n -fjL- E m )j 

One can interpret Physically the spectral weight as follows. It has two pieces, 
the first one for excited states with one more particle, and the second one for 
excited states with one more hole. Photoemission experiments (See Einstein's 
Nobel prize) access this last piece of the spectral weight, while Bremsstrahlung 
inverse spectroscopy (BIS) experiments measure the first piece. 1 Excited particle 
states contribute to positive frequencies oj' if their excitation energy is larger than 
the chemical potential, E m — E n > /i and to negative frequencies otherwise. Zero 
frequency means that the excitation energy is equal to the chemical potential. 
In other words, every excited single-particle or single-hole state corresponds to a 
delta function in the spectral weight whose weight depends on the overlap between 
initial states with one more particle at r' or one more hole at r, and the true excited 
states. 

Remark 39 At zero temperature, we have 

A{r,r';J) = Z' 1 £ [(0| V s W H (m\ 4 (r') |0> 2tt5 (w' - (E m - p - E )) 

m 

+ (0| V f s (r') \m) (m\ ^ s (r) |0) 2ir6 (w ; - (E - M - E m ))\ (3.178) 

Since to add a particle we need at least an energy fi, then E m — fi — E > and 
the first term, that adds particles, contributes only to positive frequencies. On the 
other hand, if we remove a particle, E — \i — E m < since the thermodynamic 
result dE/dN — \i implies, with dN = — 1, that E — E m < fi (the equality being 
satisfied only for the hole infinitesimally close to the Fermi surface) so that the 
second term, that adds holes, contributes only to the negative frequencies. 



1 To be more specific, these experiments add or remove particles in momentum, not position 

(t) 
S 



eigenstates. The only change that this implies in the discussion above is that i/>g' (r) should be 



replaced by Cp\ 
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The spectral representation Eq.(3.122) immediately tells us that the poles of 
the single-particle Green's functions are at the same position as delta functions 
in the spectral weight, in other words they are at the excited single-particle or 
single-hole states. Doing changes of dummy summation indices we can arrange so 
that it is always (n| that has one less particle. Then, 



A (r, r'; to') = eP n £„ m (e"^» + e"^) (n\ V> s (r) \m) (m\ ^ (r') |n) 2nS (to' - (K m - K n )) 



(3.179) 

Substituting in the spectral representation Eq.(3.122) we have, 



Q (r, r'; ik n ) = e? a E mn (e^ K - + e"^») W» s WM(H^(^)l"> 



ik„-(E m -E n -ti) 



(3.180) 



This is the Lehman representation. It tells us how to interpret the poles of the 
analytically continued Q (r, r'; ik n ) . 

Remark 40 Standard way of proving analytical continuation formula: The stan- 
dard way of proving that G R (u>) = limn. n ^ w +ir, Q (ik n ) is to first find the Lehman 
representation for both quantities. 



3.5.3 Probabilistic interpretation of the spectral weight 



For a different representation, for example for momentum, we have [15] in the 
translationally invariant case, by analogy with the above result for the spectral 
weight Eq.(3.179) 

A (k, a/) = e^ E m „ (e-0 K ™ + e~^") | (n| c k |m)| 2 2n6 (a/ - (K m - K n )) . 

(3.181) 

The overlap matrix element |(n|ck|m)| 2 that gives the magnitude of the delta 
function contribution to the spectral weight represents the overlap between the 
initial state with one more particle or hole in a momentum eigenstate and the 
true excited one-particle or one-hole state. The last equation clearly shows that 
A (k,a/) / (2n) is positive and we already know that it is normalized to unity, 

J ^A(k,a/) = ({c k ,4}) = l. (3.182) 

Hence it can be interpreted as the probability that a state formed from a true eigen- 
state \n) either by adding a particle in a single-particle state k, namely c\. \n) (or 
adding a hole Ck \ri) in a single- particle state k) is a true eigenstate whose energy 
is oj above or below the chemical potential. Clearly, adding a particle or a hole 
in a momentum eigenstate will lead to a true many-body eigenstate only if the 
momentum of each particle is individually conserved. This occurs only in the non- 
interacting case, so this is why the spectral weight is then a single delta function. 
In the more general case, many energy eigenstates will have a non-zero overlap 
with the state formed by simply adding a particle or a hole in a momentum eigen- 
state. While particle-like excitations will overlap mostly with eigenstates that are 
reached by adding positive ui, they can also overlap eigenstates that are reached 
by adding negative u. In an analogous manner, hole-like eigenstates will be mostly 
at negative cj. Let us see how this manifests itself in a specific experiment. 

Remark 41 Energy vs momentum in an interacting system: It is clear that in an 
interacting system one must distinguish the momentum and the energy variables. 
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The energy variable is uj. Knowing the momentum of a single added electron or 
hole is not enough to know the added energy. This added energy would be k 2 /2m 
only in the case of non-interacting electrons. 

Remark 42 Physical reason for high-frequency fall-off: The explicit expression 
for the spectral weight Eq.(3.181) suggests why the spectral weight falls off fast 
at large frequencies for a given k,as we have discussed in Subsection (3.6.1). A 
state formed by adding one particle ( or one hole ) of momentum k should have 
exponentially small overlap with the true eigenstates of the system that have one 
more particle ( or hole ) but an arbitrarily large energy difference uj with the initial 
state. 



3.5.4 Angle-resolved photoemission spectroscopy (ARPES) on a Fermi liquid com- 
pound. 

In a photoemission experiment, a photon ejects an electron from a solid. This is 
nothing but the old familiar photoelectric effect. In the angle resolved version of 
this experiment (ARPES), the energy and the direction of the outgoing electron 
are measured. This is illustrated in Fig. (3- 10). The outgoing electron energy can 
be measured. Because it is a free electron, this measurement gives the value of the 
wave vector through k 2 /2m. Using energy conservation, the energy of the outgoing 
electron is equal to the energy of the incident photon E p h, minus the work function 
W plus the energy of the electron in the system, ui, measured relative to the Fermi 
level. 




Figure 3-10 Schematic representation of an angle-resolved photoemission experi- 
ment. 

The energy of the electron in the system to will be mostly negative. The value 
of k|| may be extracted by simple geometric considerations from the value of k. 
Since in this experiment there is translational invariance only in the direction 
parallel to the plane, this means that in fact it is only the value of k|| that is 
conserved. Hence, it is only for layered systems that we really have access to 
both energy uj and total momentum k|| of the electron when it was in the system. 
Without going into details of the assumptions going into the derivation, Fermi's 
golden rule suggests, (see first section of Chapter 2) that the cross section for 
ejecting an electron of momentum k|| and energy uj (measured with respect to fx) 
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is proportional to 

oc Y^ e ~ PKm K n l c k,i \m)\ 2 5{u + n-{E m -E n )) (3.183) 



8 2 a 



ran 



« ^e-^" |(n|c k|| |m)| 2 5( W -(^ m -^„)) (3.184) 

cx J dte™' (cL uCku (i)) (3.185) 

which is "half of the spectral weight (the mostly negative-energy part). More 
specifically, we can rewrite this result in terms of the spectral weight as follows, 

cx f(u) A(k {l ,u) . (3.186) 



dQdcu 



Note the analogy with the cross section we found earlier for electron scattering. 
It was expressed in Eq.(2.1) in terms of a density-density correlation function. 
Through the fluctuation-dissipaton theorem Eq.(2.113), it was related to a spectral 
weight x" (analogous to A) through the Bose function (analogous to the Fermi 
function). 

Proof: The most direct and simple proof is from the Lehman representation 
Eq.(3.181). To get a few more general results about G< (k|| , to) and G > (k|| , w) 
we present the following alternate proof. The cross section is proportional 
to the Fourier transform of G < (k||,w)as defined in Eq.(3.84). 

<h(T a-iG < (k||,w) (3.187) 



One can relate G < and G > to the spectral weight in a very general way 
through the Fermi function. This is done using the usual cyclic property of 
the trace (fluctuation-dissipation theorem). From 

Ck„(*)i,) = Z-^le-^ie™^-™)^] (3-188) 

(eP K e-t K ) 4 e-' 3K (e im c k|| e-^)kl89) 



'II. 

= Z- lr Fr 



= (4, ,<*„(*+#)) (3- 19 °) 

one finds by simple use of definitions and change of integration variables, 

A(k|,, w ) = fdte^^c^ (i)+c k|| (i)c k|| ) (3.191) 

= / dte™ (4 ||Ck|| (t)) + J dte^+^V ( c t ||Ck|| (t + i/3)) 

= (l + e^) J die*** (4 |( <*,,(*)) (3.192) 

= /(w)- 1 (-iG < (k N ,w)) (3.193) 
Substituting in Eq.(3.187) proves Eq.(3.186). Note that since 

A(k||,w) =-i[G<(k ll ,uj) -G>(k||,a;)] (3.194) 
we also have the result 

iG>(k n ,u) =(l-/(w))^(k N ,w) (3.195) 
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The theoretical formula for the photoemission cross-section Eq.(3.186) neglects 
processes where energy is transferred from the outgoing electron to phonons or 
other excitations before it is detected (multiple scattering of outgoing electron). 
Such processes are referred to as "inelastic background" . 

The state of technology and historical coincidences have conspired so that the 
first class of layered (quasi-two-dimensional) compounds that became available 
for ARPES study around 1990 were high temperature superconductors. These 
materials have properties that make them non-conventional materials that are not 
yet understood using standard approaches of solid-state Physics. Hence, people 
started to look for two-dimensional materials that would behave as expected from 
standard models. Such a material, semimetallic TiTe-2 was finally found around 
1992. For our purposes, quasi-to-dimensional just means here that the Fermi 
velocity perpendicular to the planes is much smaller than the Fermi velocity in 
the planes. The results of this experiment [19] appear in Fig. (3-11). 



Fermi liquid line shape fits 




Energy relative to E F (meV) 

Figure 3-11 ARPES spectrum of 1 - T - TiTe 2 , after R. Claessen, R.O. Anderson, 
J.W. Allen, C.G. Olson, C. Janowitz, W.P. Ellis, S. Harm, M. Kalning, R. Manzke, 
and M. Skibowski, Phys. Rev. Lett 69, 808 (1992). 

We have to remember that the incident photon energy is 21.2eV while the 
variation of cu is on a scale of 200meV so that, for all practical purposes, the 
momentum vector in Fig. (3-10) is a fixed length vector. Hence, the angle with 
respect to the incident photon suffices to define the value of k||. Each curve in 
Fig. (3-11) is for a given kii, in other words for a given angle measured from the 
direction of incidence of the photon. The intensity is plotted as a function of the 
energy of the outgoing electron. The zero corresponds to an electron extracted 
from the Fermi level. Electrons with a smaller kinetic energy come from states 
with larger binding energy. In other words, each of the curves above is basically 
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a plot of the hole- like part of A (k||,w). From band structure calculations, one 
knows that the angle 9 — 14.75° corresponds to the Fermi level (marked Uf on the 
plot) of a Ti — 3d derived band. It is for this scattering angle that the agreement 
between experiment and Fermi liquid theory is best (see Sec. (3.5.6) below). The 
plots for angles 8 < 14.75° corresponds to wave vectors above the Fermi level. 
There, the intensity is much smaller than for the other peaks. For = 13°, 
the experimental results are scaled up by a factor 16. The intensity observed for 
wave-vectors above the Fermi wave vector comes from the overlap of particle-like 
excitations with eigenstates below the Fermi surface, a phenomenon we alluded to 
in the previous section. 

The energy resolution is 35meV. Nevertheless, it is clear that the line shapes 
are larger than the energy resolution: Clearly the spectral weight is not a delta 
function and the electrons in the system are not free particles. Nevertheless, 
there is a definite maximum in the spectra whose position changes with k||. It 
is tempting to associate the width of the line to a lifetime. In other words, a 
natural explanation of these spectra is that the electrons inside the system are 
"quasiparticles" whose energy disperses with wave vector and that have a lifetime. 
We try to make these concepts more precise below. 



3.5.5 Quasiparticles[17] 



For a general interacting system, the one-particle Green's function takes the form, 

1 



G (k,w) = 



w + ir \ - Ck - T, R ( k > w ) 

The corresponding spectral weight is, 
A(k,oj) = -21mG R (k,u) 

-2Im£ K (k, W ) 



(^-C k -ReE K (k,^)) + (lm£^(k,^)' 



(3.196) 

(3.197) 
(3.198) 



If the imaginary part of the self-energy, the scattering rate, is not too large and 
varies smoothly with frequency, the spectral weight will have a maximum whenever 



cj-C k -ReS fi (k,w) =0 



(3.199) 



Let Ek — \x be the value of u for which this equation is satisfied. E^ is the so-called 
quasiparticle energy. This energy is clearly in general different from the results of 
band structure calculations that are usually obtained by neglecting the frequency 
dependence of the self-energy. Expanding w — ( k — Re T, R (k, w) around oj — E^ + fi 
where A (k,w) is a maximum, we find 



w-C k -ReE fl (k,o;) sa 1 - 



<9Re£ fl (k,w) 



(oj-E k + n) + ... (3.200) 



Ek-fj,, 



If we define the "quasiparticle weight" or square of the wave function renormal- 
ization by 



■ ReS R (k,w)| 



(3.201) 



= E k-f 
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then in the vicinity of the maximum, the spectral weight takes the following simple 
form in the vicinity of the Fermi level, where the peak is sharpest 



A(k,w) 



2irZ k - 



-4ImE (k,w) 



27rZ k 



£k + M) 2 + (z k ImE K (k,^))' 

r k (w) 



i 



(W - S k 



M ) 2 + (r k H) 2 



+ mc (3.202) 



(3.203) 



The last equation needs some explanation. First, it is clear that we have defined 
the scattering rate 

(3.204) 



F k (u) = -Z k ImE B (k,w) 



Second, the quantity in square brackets looks, as a function of frequency, like a 
Lorentzian. At least if we can neglect the frequency dependence of the scattering 
rate. The integral over frequency of the square bracket is unity. Since A (k,w) /2ir 
is normalized to unity, this means both that 

Z k < 1 (3.205) 

and that there are additional contributions to the spectral weight that we have 
denoted inc in accord with the usual terminology of "incoherent background" . 
The equality in the last equation holds only if the real part of the self-energy is 
frequency independent. 

It is also natural to ask how the quasiparticle disperses, in other words, what is 
its effective Fermi velocity compared with that of the bare particle. Let us define 
the bare velocity by 

vu = V k C k (3.206) 

and the renormalized velocity by 

< = V k £ k (3.207) 

Then the relation between both quantities is easily obtained by taking the gradient 
of the quasiparticle equation Eq.(3.199) 

C k -ReE«(k,.E k - M )=0] 

dReZ R (k,w) 



V k [E y 
vt-vk- V k ReE fi (k,£ k 



M) 



vt = 



(3.208) 
(3.209) 



E k -n 



where V k in the last equation acts only on the first argument of Re T, R (k, £? k — fj,). 
The last equation is easily solved if we can write that k dependence of E fl as 
a function of £ k instead, something that is always possible for spherical Fermi 
surfaces. In such a case, we have 



l+ s | r ReS R (k,B k - j u) 



- Re£ R (k,i 



(3.210) 



In cases where the band structure has correctly treated the k dependence of the 
self-energy or when the latter is negligible, then the renormalized Fermi velocity 
differs from the bare one only through the famous quasiparticle renormalization 
factor. In other words, u£ = Z k v k . The equation for the renormalized velocity 
is also often written in terms of a mass renormalization instead. Indeed, we will 
discuss later the fact that the Fermi wave vector fcp is unmodified by interactions 
for spherical Fermi surfaces (Luttinger's theorem). Defining then m*v* kF = k F — 
mVk F means that our equation for the renormalized velocity gives us 



(3.211) 
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3.5.6 Fermi liquid interpretation of ARPES 



Let us see how to interpret the experiments of the previous subsection in light 
of the quasiparticle model just described. First of all, the wave vectors studied 
are all close to the Fermi surface as measured on the scale of kp. Hence, every 
quantity appearing in the quasiparticle spectral weight Eq. (3.203) is evaluated 
for k = kp so that only the frequency dependence of the remaining quantities is 
important. The experiments were carried out at T = 20K where the resistivity has 
a T 2 temperature dependence. This is the regime dominated by electron-electron 
interactions, where so-called Fermi liquid theory applies. What is Fermi liquid 
theory? It would require more than the few lines that we have to explain it, but 
roughly speaking, for our purposes, let us say that it uses the fact that phase space 
for electron-electron scattering vanishes at zero temperature and at the Fermi 
surface, to argue that the quasiparticle model applies to interacting electrons. 
Originally the model was developed by Landau for liquid 3 He which has fermionic 
properties, hence the name Fermi Liquid theory. It is a very deep theory that in 
a sense justifies all the successes of the almost-free electron picture of electrons in 
solids. We cannot do it justice here. A simple way to make its main ingredients 
plausible, [18] is to assume that near the Fermi surface, at frequencies much less 
than temperature, the self-energy is i) analytic and ii) has an imaginary part that 
vanishes at zero frequency. The latter result follows from general considerations 
on the Pauli principle and available phase space that we do not discuss here. 
Let us define real and imaginary parts of the retarded self-energy by 

E fl = £' + iE" (3.212) 

Our two hypothesis imply that E" has the Taylor expansion 

E" (k F ;w) = aw-7W 2 + ... (3.213) 

The imaginary part of the retarded self-energy must be negative to insure that the 
retarded Green's function has poles in the lower half-plane, as is clear from the 
general relation between Green function and self-energy Eq.(3.196). This means 
that we must have a = and 7 > 0. Fermi liquid theory keeps only the leading 
term 

E" = -JLU 2 

We will verify for simple models that this quadratic frequency dependence is es- 
sentially correct in d > 3. The real part is then obtained from the Kramers-Kronig 
relation Eq.(??), (Sec. 3. 6. 2) or from the spectral representation, 

hm[E'(k F ; W )-E'(k F ;oo)] = lim V f ^ ^ ^ (3.214) 

= <p / \JLl — L +LU p / v *l ' (3rgl5) 

i t oj' [ J f (a/) 2 ] 

The first term is the value of the real-part of the self-energy at zero-frequency. This 
constant contributes directly to the numerical value of the chemical potential (the 
Hartree-Fock shift E' (kp; 00) does not suffice to evaluate the chemical potential). 
The second term in the last equation tells us that 



|^'(k^) 



w=0 



V 



/dec 
7T 



dw' E " (k F ;h/) 



(3.216) 
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Since S" = — juj 2 the integral exists and is negative (if we assume a frequency 
cutoff as discussed below), hence 



< (3.217) 

w=0 



This in turn means that the corresponding value of Z kp is less than unity, as we 
had concluded in Eqs. (3.201) and (3.205) above. In summary, the analyticity hy- 
pothesis along with the vanishing of S" (0) implies the existence of quasiparticles. 

Remark 43 Warning: there are subtleties. The above results assume that there 
is a cutoff to S" (kp;u/) . The argument just mentioned in Eq. (3.216) fails when 
the integral diverges. Then, the low frequency expansion for the self-energy in 
Eq.(3.215) cannot be done. Expanding under the integral sign is no longer valid. 
One must do the principal part integral first. In fact, even for a Fermi liquid at 
finite temperature, T," (k F ;uj) ~ lo 2 + (nT) 2 so that the (ttT) 2 appears to lead 
to a divergent integral in Eq.(S.216). Returning to the original Kramer s-Kronig 
expression fo £' however, the principal part integral shows that the constant term 
(nT) 2 forT," (kp;w) does not contribute at all to £' if the cutoff in S" is symmetric 
at positive and negative frequencies. In practice one can encounter situations 
where dT,/dui > 0. In that case, we do not have a Fermi liquid since Z > 1 is 
inconsistent with the normalization of the spectral weight. One can work out an 
explicit example in the renormalized classical regime of spin fluctuations in two 
dimensions. (Appendix D of [28]). 

The solid lines in Fig. (3-11) are two-parameter fits that also take into account 
the wave vector and energy resolution of the experiment [19]. One parameter is 
Ef. — /j, while the other one is 7', a quantity defined by substituting the Fermi 
liquid approximation in the equation for damping Eq. (3.204) 

T kF (w) = Z kpl u 2 = iu 2 . (3.218) 

Contrary to E k , the damping parameter 7' is the same for all curves. The solid-line 
fits are obtained with 7' = AOeV^ 1 (/3 on the figure). The fits become increasingly 
worse as one moves away from the Fermi surface, as expected. It is important to 
notice, however, that even the small left-over weight for wave-vectors above the 
Fermi surface (0 < 14.75°) can be fitted with the same value of 7. This weight is 
the tail of a quasiparticle that could be observed at positive frequencies in inverse 
photoemission experiments (so-called BIS). The authors compared the results of 
their fits to the theoretical estimate, [20] 7 = 0.067cj p /e|,. Using uj p = 18.2eV, 
sf = 0.3eV and the extrapolated value of Z kp obtained by putting 2 r s — 10 in 
electron gas results, [21] they find 7' < 5(eV)^ 1 while their experimental results 
are consistent with 7' = 40 ± 5 (eV) 1 . The theoretical estimate is almost one 
order of magnitude smaller than the experimental result. This is not so bad given 
the crudeness of the theoretical model (electron gas with no lattice effect). In 
particular, this system is a semimetal so that there are other decay channels than 
just the one estimated from a single circular Fermi surface. Furthermore, electron 
gas calculations are formally correct only for small r s while there we have r s = 10. 

More recent experiments have been performed by Grioni's group [27]. Results 
are shown in Fig. (3-12). In this work, authors allow for a constant damping 
r = 17 meV coming from the temperature and from disorder and then they 
fit the rest with a Fermi velocity Hvp = 0.73 ± O.leV A close to band structure 
calculations, Hvp — 0.68 eV A and 7' that varies between 0.5 eV~ x (16°) and 
0.9 eV~ x (14.5°). The Fermi liquid fit is just as good, but the interpretation of 
the origin of the broadening terms is different. This shows that there is much 
uncertainty still in the interpretation of ARPES data, even for Fermi liquids. 

2 r s is the average electron spacing expressed in terms of the Bohr radius. 
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Figure 3-12 Figure 1 from Ref.[27] for the ARPES spectrum of 17-TiTe 2 measured 
near the Fermi surface crossing along the high-symmetry TM direction (6 = is 
normal emission). The lines are results of Fermi liquid fits and the inset shows a 
portion of the Brillouin zone with the relevant ellipsoidal electron pocket. 



Theoretical estimates for high-temperature superconductors are two orders of 
magnitude smaller than the observed result [19]. 

Remark 44 Asymmetry of the lineshape: The line shapes are asymmetrical, with 
a tail at energies far from the Fermi surface (large binding energies). This is 
consistent with the fact that the "inverse lifetime" T kF (to) = Zk F ^u) 2 is not a 
constant, but is instead larger at larger binding energies. 

Remark 45 Failure of Fermi liquid at high-frequency: Clearly the Fermi liquid 
expression for the self-energy fails at large frequencies since we know from its 
spectral representation that the real-part of the self-energy goes to a frequency- 
independent constant at large frequency, the first correction being proportional to 
1/ikn as discussed below in subsection (3.6.1). Conversely, there is always a cutoff 
in the imaginary part of the self-energy. This is not apparent in the Fermi liquid 
form above but we had to assume its existence for convergence. The cutoff on 
the imaginary part is analogous to the cutoff in x" ■ Absorption cannot occur at 
arbitrary high frequency. 

Remark 46 Destruction of quasiparticles by critical fluctuations in two dimen- 
sions: Note that it is only if E" vanishes fast enough with frequency that it is 
correct to expand the Kramers- Kronig expression in powers of the frequency to 
obtain Eq.(3.216). When E" (co) vanishes slower than ui 2 , then Eq.(3.216) for the 
slope of the real part is not valid. The integral does not converge uniformly and it 
is not possible to interchange the order of differentiation and integration. In such 
a case it is possible to have the opposite inequality for the slope of the real part 
^jE'(k,w)| > 0. This does not lead to any contradiction, such as Z kp > 1, 
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because there is no quasiparticle solution at uj = in this case. This situation 
occurs for example in two dimensions when classical thermal fluctuations create a 
pseudogap in the normal state before a zero-temperature phase transition is reached 
[22]. 



3.5.7 Momentum distribution in an interacting system 

In an interacting system, momentum is not a good quantum number so (c^cy^ is 

not equal to the Fermi distribution. On the other hand, (^c^c-^ can be computed 
from the spectral weight. Indeed, 



-(T T c k (r)ct 



c^Ck) = lim 

/ T— >0 _ 

= lim 0(k,r) (3.219) 

T— >0 _ 

To compute the latter quantity from the spectral weight, it suffices to use the 
spectral representation Eq.(3.122) 



lim G(k,T) = T lim V e- 4fc " T £(r,r';ifc„) 

71— — OO 

Tim t e-^r^'H^) 

n — — r*-i 



2tt ik„ — u: 1 



(3.220) 



Using the result Eq.(3.163) found above for the sum over Matsubara frequencies, 
we are left with 



c k c k ) = lim_ - Q (k,r) = (a/) A (k,o/) 



(3.221) 



with / (u/) the Fermi-Dirac distribution. 

This means that the momentum distribution is a Fermi-Dirac distribution only 
if the spectral weight is a delta function. This occurs for free particles or, more 
generally if the real-part of the self-energy is frequency independent since, in this 
case, the Kramers-Kronig relations imply that the imaginary part of the self-energy 
vanishes so that Eq.(3.198) for the spectral weight gives us a delta function. 

Remark 47 Jump of the momentum distribution at the Fermi level: Even if 
( c iL c k^) * s no-longer a Fermi-Dirac distribution in an interacting system, neverthe- 
less at zero-temperature in a system subject only to electron- electron interaction, 
there is a jump in (^c^Ck^ at the Fermi level. The existence of this jump can be 
seen as follows. At zero temperature, our last result gives us 

c W) = f_ % A ^ J ) ( 3 - 222 ) 

Let us take the quasiparticle form Eq. (3.203) of the spectral weight with the Fermi 
liquid expression Eq. (3.218) for the scattering rate. The incoherent background 
varies smoothly with k and hence cannot lead to any jump in occupation number. 
The quasiparticle piece on the other hand behaves when k — > hp, or in other words 
when E k — [i^ 0, as Z kF 8(uj). At least crudely speaking. When E k — [i^ 0~, this 
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delta function is inside the integration domain hence it contributes to the integral, 
while when E^ — fj, — > + the delta function is outside and does not contribute to 
the integral. This means that there is a big difference between these two nearby 
wave vectors, namely 

lirn_ (^Ck^ - lim + ^c k Ck^ = Z kp (3.223) 

In the above argument, we have done as if T k (w) was frequency independent and 
infinitesimally small in Eq. (3.203). This is not the case so our argument is rather 
crude. Nevertheless, if one uses the actual frequency- dependent forms and does 
the frequency integral explicitly, one can check that the above conclusion about the 
jump is true (although less trivial). 

Remark 48 Fermi surface and interactions: The conclusion of the previous re- 
mark is that even in an interacting system, there is a sharp Fermi surface as in 
the free electron model. For simplicity we have discussed the spinless case. A 
qualitative sketch of the zero-temperature momentum distribution in an interact- 
ing system appears in Fig. (3-13). Since momentum of a single particle is not 
a good quantum number anymore, some states above the Fermi momentum are 
now occupied while others below are empty. Nevertheless, the Fermi surface is 
unaffected. 



1 - 




Figure 3-13 Qualitative sketch of the zero-temperature momentum distribution in 
an interacting system. 



Remark 49 Luttinger's theorem: More generally, in a Fermi liquid the volume 
of reciprocal space contained within the Fermi surface defined by the jump, is inde- 
pendent of interactions. This is Luttinger's theorem. In the case where the Fermi 
surface is spherical, this means that kp is unaffected. 



3.6 A few more formal matters : asymptotic behav- 
ior and causality 

In designing approximations, we have to try to preserve as many as possible of the 
exact properties. Sum rules are such properties. They determine the structure 
of the high-frequency expansion and hence one can also check whether a given 
approximation preserves the sum rules by looking at its high-frequency expansion. 
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This is the first topic we will discuss. The second topic concerns restrictions 
imposed by causality. This has become a very important topic in the context of 
Dynamical Mean-field theory or other approaches that describe the physics that 
occurs at strong coupling, such as the Mott transition. We will come back on this 
in later chapters. 



3.6.1 Asymptotic behavior of Q (k:ik n ) and S (k;ik n ) 



As usual, the high-frequency asymptotic properties of the Green's function are 
determined by sum rules. From the spectral representation(3.122), we obtain, for 
the general interacting case 

JimJ^tt^Jim . r^^l (3.224) 
= lim — / _A(k;w')= Um — (l c k, <£}) = lim — (3.225) 

ife„^oo lk n v /_ 00 2lT ifc„^oo lk n ifc„^oo lk n 

Defining the self-energy as usual 

the correct asymptotic behavior for the Green's function implies that the self- 
energy at high frequency cannot diverge: It must go to a constant independent of 
frequency 

lim S (k, ik n ) = est. (3.227) 

ik n — >oo 

We will see later that the value of this constant is in fact given correctly by the 
Hartree-Fock approximation. 

The converse of the above result [18] for the Green's function, is that if 

lim G (k;ik n ) = lim — !— 

then that is all that is needed to obtain an approximation for the Green's function 
which obeys the anticommutation relation: 

G (k;(T) - g (k;0+) = (c+c k ) + (c k c+) = 1 (3.228) 

Proof :It suffices to notice that 

g (k;(T) - g (k;0+) = - [e- lk - {r - e -^- 0+ ] Q (k;ik n ) (3.229) 

^ n 

We can add and subtract the asymptotic behavior to obtain, 

\ ? - ( e **> "cK? ('^ - ) 56 

(3.230) 

In the first sum, g (k;ik n ) — ^- decays faster than ^- so that the convergence 
factors are not needed for the sum to converge. This means that this first 
sum vanishes. The last sum gives unity, as we easily see from the previous 
section. This proves our assertion. 
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Remark 50 High-frequency expansion for the Green's function and sum-rules: 
The coefficients of the high-frequency expansion ofQ(k;ik n ) in powers ofl/ik n 
are obtained from sum rules on the spectral weight, in complete analogy with what 
we have found in previous chapters. The fact that A (k,w) falls fast enough to 
allow us to expand under the integral sign follows from the fact that all frequency 
moments of A (k,w) , namely J (kou n A (k,w) , exist and are given by equal-time 
commutators. Explicit expressions for A (k,cj) in terms of matrix elements, as 
given in Subsection(S.5.3) above, show physically why A (k,cj) falls so fast at 
large frequencies. As an example, to show that the coefficient of the l/ik n term 
in the high frequency expansion is equal to ^-A(\l;uj') it is sufficient that 



J^oo Itn \ u ' A (k; exists. [16] This can be seen as follows, 



ik n Q (k; ik n ) 
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J — ( 



dm' 
27 
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doj' 
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A(k;oj') 



ik n 



—A (k;u/) 

2tt ik n — u)' 
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rf.231) 



(3.232) 



(3.233) 



AQt;u/)u/\ (3.234) 



If the integral exists then, it is a rigorous result that 

~27 



lim ik ri 

ik n — >oo 



'(k:ik n ) = / 

J — ( 



A(k;uj') 



(3.235) 



This is an important result. It suggests that approximate theories that give 1 as the 
coefficient of (ik n ) 1 in the high frequency expansion have a normalized spectral 
weight. However [16] the above proof assumes that there is indeed a spectral repre- 
sentation for Q (k; ik n ) . A Green's function for a theory that is not causal fails to 
have a spectral representation. If a spectral representation is possible, the analyti- 
cally continued approximate G R (k,w) is necessarily causal. Approximate theories 
may not be causal. This failure of causality may reflect a phase transition, as we 
will see later, or may simply be a sign that the approximation is bad. As an exam- 
ple, suppose that we obtain Q (k;ik n ) — (ik n — ia) 1 . This has the correct high- 
frequency behavior but its analytical continuation does not satisfy causality. It has 
no spectral representation. On the other hand, Q (k; ik n ) = (ik n + (k n / \k n \) ia) 1 
has a Lorentzian as a spectral weight and is causal. It may also occur that the 
approximate theory may have ^-A(\i.;La') = 1 but A(k;uj') < for some 
range of uJ . This unphysical result may again signal that the approximate theory 
fails because of a phase transition or because it is a bad approximation. 



3.6.2 Implications of causality for G R and 



Consider the retarded Green function as a matrix in r, r'. We will show that 
the real and imaginary parts of G R and of T, R are each Hermitian matrices, in 
addition, lmG R and ImS^ are both negative definite (except in the special case 
of non-interacting particles where ImT, R = 0). 

In analogy with the Matsubara Green function Eq.(3.180) G R has the Lehman 
representation 



G R (r, r» = e^ £ (e"^™ + e"^) 



MV> g (r) \m) (m[4 {r')\n 
uj + in- (E m - E n - n) 



(3.236) 
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Taking the complex conjugate of the transpose, we obtain 
ReG«(r',r;u,r = e*» £ (e~^ + e~^) H { ^}™ T ^ ( ? |n) * 

mn v r 7 

= ReG R (r,r';u) (3.237) 
ImG K (r',r; W ) 4 = -tt ^ (e"^"* + e"^") <n| Vs (r') |m)* <m| ^ s (r) \n)* 



mn 



xS (u — (E m — E n — (£)) (3.238) 
= ImG s (r,r';w) (3.239) 

which means that the real and imaginary parts are both Hermitian matrices. In 
a basis where the matrix lmG R (r, r';w) is diagonal, say for quantum number a, 
then 

ImG fi (a;w) = -tt ^ {e~ 0Km + e~ pKn ) (n\ c a \m) (m\ 4 \n) 6 (w - (£ m - E n - 

mn 

= -7TY / (^ Km +e- 0K ")\(n\c a \m)\ 2 d^-(E m -E n -^) (3. 

mn 

which proves that the matrix for the imaginary part is negative definite in addition 
to being Hermitian. The negative sign comes from the +irj in the original formula 
and is clearly a consequence of causality. 

Following Potthoff [26] we show that the retarded self-energy as a matrix has 
the same properties as G R . First, we need to prove that 

1 =XtiY (3.241) 



A±iB 

with X and Y both Hermitian and Y positive definite if A and B are both Her- 
mitian with B positive definite. This is true because 

' B-^—-—l—^B-^. (3.242) 



A±iB B~ 1 /' 2 AB- 1 / 2 ±iI 

Since £? _1 / 2 AB -1 / 2 is Hermitian as well, we can diagonalize it by a unitary trans- 
formation B~ x l 2 AB~ X I 2 = UcW where c is a diagonal matrix. Thus, 

= B-^U-^rfB- 1 ' 2 = B-^U^-U^B- 1 / 2 =X T iY (3.243) 
A ± %B c±il c 2 + 1 

with X and Y Hermitians since (l/tB-V2) t = B-^ 2 U. In addition, Y is positive 
definite since in the diagonal basis Y — > (c 2 + l) 1 . Now, define (G R ) 1 = 
(A - iBy 1 = X + iY and (G^)" 1 = (A - iB^y 1 = X + iY so that 

(G fi ) _1 =X + iY= (G*)' 1 - E fi = JT + iY - ReX R - iIm£ K (3.244) 

Then, given that X, Y and Xo, Yo ar e Hermitians, we have that Re T, R and Im S K 
are Hermitians. In addition, ImE fi is negative definite since Y$ is infinitesimal 
which implies that Y — Y can only be positive (or vanish in the non-interacting 
case). 



3.7 Three general theorems 



Risking to wear your patience out, we still have to go through three general the- 
orems used repeatedly in Many-Body theory. Wick's theorem forms the basis ot 
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the diagram technique in many-body theory. The linked-cluster theorems, or cu- 
mulant expansions, are much more general theorems that are also necessary to set 
up the machinery of diagrams. Finally, we state a variational principle for the free 
energy that allows us to give a physical meaning to Hartree-Fock theory as the 
best one-body Hamiltonian for any given problem. This variational principle is 
useful for ordinary system, but also becomes indispensable when there is a broken 
symmetry. 



3.7.1 Wick's theorem 



Wick's theorem allows us to compute arbitrary correlation functions of any Hamil- 
tonian that is quadratic in Fermion or Boson operators. That is clearly what we 
need to do perturbation theory, but let us look in a bit more details at how this 
comes about. We have already talked about Wick's theorem in the context of 
Hartree-Fock theory. Here we present the more general result valid at finite tem- 
perature. We will need to compute in the interaction picture 



Q(r) = - 



Tr 



-P h °T t (u(P,t)$(t)U(t,0)$ (0))" 



TV 



e -PHoT T U(f3,0) 



(3.245) 



Because U (r, 0) always contains an even number of fermions, it can be commuted 
with creation-annihilation operators without paying the price of minus signs so 
that 



S(r) 



Tr 



-^OT T (&(/9,0)$(T)5i t (0)) 
Tr[e-> 3H 0T T U(I3,0)] 



More specifically the evolution operator is, 



U(p,0)=T T [exp(-tfdT 1 V(T 1 )) 



(3.246) 



(3.247) 



Expanding this evolution operator to first order in the numerator of the Green's 
function one obtains 



-Tr 



-^T t (#(t)^ (0))] +^ 



dnTr 



-f3 H 



T °r T (v(Ti)^(r)^ (0))" 

(3.248) 

where in the case of a two-body interaction (Coulomb for example), V (tj.) contains 
four field operators. 

Wick's theorem allows us to evaluate expectation values such as those above. 
More generally, it allows us to compute expectation values of creation-annihilation 
operators such as, 

;a i (r i )a j (r j )4(r fe )4(r0)„ (3.249) 



as long as the density matrix e~^ H ° is that of a quadratic Hamiltonian. 

Note that since quadratic Hamiltonians conserve the number of particles, ex- 
pectation values vanish when the number of creation operators does not match 
the number of destruction operators. 

Lemma 9 If H = e±a\ai + e2a\a 2 then ^aia{a 2 a|^ = (aia}^J (a^a^j . 
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Proof: To understand what is going on, it is instructive to study first the problem 
where a single fermion state can be occupied. Then 



a\a\ 



Tr 



I a ia \ |0> 



Tr [e-f> H o] 
^«0| 01)010! (410)) 



(0||0)+ e -^((0|o 1 )(a t 1 |0) 



1 + e-P^ 



(3.250) 



(3.251) 



For two fermion states 1, 2, then the complete set used to evaluate the trace 
is 

|0>|0>, 0110)10), |0)<4|0), 4 10) 4 10> (3.252) 

so that 

t\ 1 l + e-P^ 1 

°*°V = l + e-^l + e-^ = TT^- (3 - 253) 

The easiest way to understand the last result is to recall that (l + a{) (l + a^j 

will generate the trace so that we can factor each subspace. The last result 
will remain true for an arbitrary number of fermion states, in other words 



Oio{ 



Furthermore, 



1 Um^l 1 + e 



1 + 



aia\a 2 al) - 



1 



n 



m^l,2 



1 



1 + 1 + e-P** n m ^i,2 1 + e ~ pe " 

1 1 



1 + e~P ei 1 + e"^ 2 



oio} 



0202 



(3.254) 

(3.255) 

(3.256) 
(3.257) 



Theorem 10 Any expectation value such as /oj (rj) aj (tj) a\ (t&) a\ cal- 



culated with a density matrix e P K ° that is quadratic in field operators can be 



puted as the sum of all possible products of the type (aj (tj) a\ (rk)) ( Oj (Ti) a] (ti) 

that can be formed by pairing creation an annihilation operators. For a given term 
on the right-hand side, there is a minus sign if the order of the operators is an odd 
permutation of the order of operators on the left-hand side. 

Proof: It is somewhat pretentious to call a proof the plausibility argument that we 
give below, but let us go ahead anyway. The trick to prove the theorem([23]) 
is to transform the operators to the basis where H is diagonal, to evaluate 
the expectation values, then to transform back to the original basis. Let 
Greek letters stand for the basis where H is diagonal. Using the formula 
for basis changes, we have, (with an implicit sum over Greek indices) 



,t . 



a l (u) aj (tj) a{ (r fe ) a] (tj) ) = 



(i\ a) (j\ (3) (o Q ( Ti ) a f3 ( Tj ) at ( Tfc ) a\ (t,)) o (l\ k) (5\ I) 
We already know from Eq. (3.141) that 



(3.258) 



(3.259) 



O q (Ti) 



e c « Ti o n 



4 (Ti) = ale 



(3.260) 
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so that 

[ai (Ti) a, (tj) 4 (Tfc) 4 ( T 0)„ (3-261) 





= (t| a) e- c ° Tl (j| 0) e _c <» T ' ^a aa(3 44) e ^ Tfc (7l fc ) eC<sT ' ( 5 I l ) ( 3 - 262 ) 
What we need to evaluate then are expectation values of the type 



a Q a ) g44) • (3.263) 



Evaluating the trace in the diagonal basis, we see that we will obtain a non- 
zero value only if indices of creation and annihilation operators match two 
by two or are all equal. Suppose (3 = 7, a = 8 and a^/J. Then, as in the 
lemma 

(«aa /3 44) = («a4) ( a /3 a /?) ( 3 ' 264 ) 

If instead, /3 = 5, a = 7 and then 

a Q a^44) = ~ (aaapa^a)^ = - (a a a) a ) Q /a^aU . (3.265) 



'o \ H 1 u \ M / 

The last case to consider is a = /?, /3 = 5, a = 7 

(fl a s«44) = °- (3.266) 

All these results, Eqs.(3.264)(3.265) and the last equation can be combined 
into one formula 

(a a apal l a i ^ = (a Q 4)o ( a ^ a ^) - S a ^5 0iS ) (3.267) 

= ( a «4) ( a a \) o - < a «4> ( a ^ a «) ( 3 - 268 ) 

which is easiest to remember as follows, 

(a a apa^al^ = ^4^4*4^ + (^aapa\a}}j (3.269) 

in other words, all possible pairs of creation and annihilation operators must 
be paired ("contracted") in all possible ways. There is a minus sign if an 
odd number of operator exchanges (transpositions) is necessary to bring the 
contracted operators next to each other on the right-hand side (In practice, 
just count one minus sign every time two operators are permuted). Substi- 
tuting Eq.(3.268) back into the expression for the original average expressed 
in the diagonal basis Eq. (3.262) we have 

(a i (T i )a j (T j )a] s (T k )a](T l )) o (3.270) 

= (tti (Ti) 4 Cn)> (dj (Tj) 4 (Tk) )„ - (^ (Ti) 4 (Tfc) )„ (oj (tj) 4 ( T 



By induction (not done here) one can show that this result generalizes to the 
expectation value of an arbitrary number of creation-annihilation operators. 

Definition 11 Contraction: In the context of Wick's theorem, we call each factor 
(ai (Ti)a\ (rfe)^ on the right-hand side, a "contraction" . 
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Since Wick's theorem is valid for an arbitrary time ordering, it is also valid for 
time-ordered products so that, for example 



H (Ti)aj (Tj)4 ( T k)4 (ti) 



J / 



(3.271) 



a 3 ( T i)4 ( T k) 



J / 



-<r T 



di (Ti) a{ (rfe) 



a j ( T j)4 ( T l) 



(3.272) 

The only simplification that occurs with time-ordered products is the following. 
Note that, given the definition of time-ordered product, we have 



a i ( T i) a l ( T fc) 



J / o 



4 (Tfc)«i (Ti) 



J / 



(3.273) 



Indeed, the left-hand side and right-hand side of the above equation are, respec- 
tively 



a-i (t;) 4 (T fe ) 

4 ( T k)at (Ti) 



= (ai(Ti)al(T k )) o 6(Ti-T k ) (3.274) 

- (4 ( Tfc ) a t (t;)) o ( Tk - Ti) (3.275) 

= -(al{T k )ai{Ti)) o 6(T k -Ti) (3.276) 

+ (a % ( Ti ) 4 (r fc )) o - T k ) (3.277) 



In other words, operators can be permuted at will inside a time-ordered product, in 
particular inside a contraction, as long as we take care of the minus-signs associated 
with permutations. This is true for time-ordered products of an arbitrary number 
of operators and for an arbitrary density matrix. 

On the other hand, if we apply Wick's theorem to a product that is not time 
ordered, then we have to remember that 



ai (Ti) 4 ( T fc)) + - (4 ( T k) a » () 



(3.278) 



as we can easily verify by looking at the special case T k = Ti or by going to a 
diagonal basis. We can anticommute operators at will to do the "contractions" 
but they cannot be permuted inside a contraction (ai (ji) a\ (r^)^ . 

In practice, we will apply Wick's theorem to time-ordered products. In nu- 
merical calculations it is sometimes necessary to apply it to objects that are not 
time-ordered. 

Example 12 To make the example of Wick's theorem Eq.(3.272) more plausible, 
we give a few examples, Suppose first that the time order to the left of Eq. (3.272) 
is such that the destruction operators are inverted. Then, 

(t t ^ (Ti) a, (tj) 4 (Tk) 4 (ti) ^ = - (aj (tj) ai (n) a\ (r fc ) a] (r z )^ (3.279) 

which means that since i and j have exchanged roles, in doing the contractions 
as above there is one more permutation to do, which gets rid of the extra minus 
sign and reproduces the right-hand side of Eq. (3. 272). More explicitly, to do the 
contractions as above, we have to change i for j on both the right- and the left-hand 
side of Eq. (3.270). Doing this and substituting above, we obtain 



aj (Tj) a, (Ti)a\ (r fc ) 4 (ti) 



a i ( T j) 4 (ti)) q ( a i ( T i) 4 ( T k))-( a i ( T i) 4 ( r fe)) Q («i ( T *) 4 (Tl) 

(3.280) 
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which we substitute in the previous equation to obtain exactly what the right-hand 
side of Eq. (3.272) would have predicted. To take another example, suppose that 
the time orders are such that 

(T T a, ( Ti ) a } (tj) 4 (Tfc) a\ (ti) ^ = - (a, ( Ti ) a\ (r fe ) aj (tj) a] (r;)^ . 

(3.281) 

Then, to do the contractions we proceed as above, being careful not to permute 
creation and annihilation operators within an expectation value 



v a » (t»)4 (Tk)aj {Tj)a\ (rj) / 

= - (Ti) 4 (Tj)) o (4 (Tfc) «i (Tj)) Q - (Oi (Ti) 4 (Tfe)) Q ( a i (Tj) 4 R^) 2 ) 

T/ie right-hand side of Eq. (3.272) gives us 



a % (r 4 )4 (t;) 



a j (Tj)4 (Tfc) 



a* (t 4 )4 (Tfc) 



a j (tj)4 On) 



« 4 (Ti) 4 Gn)) (4 (Tfc) % (Tj)) Q - (ai (Ti) 4 (Tfc)) () (aj (Tj) 4 (Ti) 7 q 

tmift i/ie minus sign in the first term because we had to exchange the order in one 
of the time-ordered products. 



J / o 
(3.283) 



3.7.2 Linked cluster theorems 



Suppose we want to evaluate the Green's function by expanding the time-ordered 
product in the evolution operator Eq. (3.247). The expansion has to be done both 
in the numerator and in the denominator of the general expression for the average 
Eq. (3.245). This is a very general problem that forces us to introduce the notion 
of connected graphs. A generalization of this problem also occurs if we want to 
compute the free-energy from 



In Z = In 



(rr[e-^ ^C8,0)]) =ln(z (c/(/3,0))J (3.284) 



In ( T T 



exp 



\nZ (3.285) 



In probability theory this is like computing the cumulant expansion of the char- 
acteristic function. Welcome to linked cluster theorems. 

These problems are special cases of much more general problems in the theory 
of random variables which do not even refer to specific Feynman diagrams or to 
quantum mechanics. The theorems, and their corollary that we prove below, are 
amongst the most important theorems used in many-body Physics or Statistical 
Mechanics in general. 



Linked cluster theorem for normalized averages 

Consider the calculation of 



(e-A*)) 



(3.286) 
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where the expectation () is computed over a multivariate probability distribution 
function for the variables collectively represented by x. The function / (x) is 
arbitrary, as is the function A (x). Expanding the exponential, we may write 



(e-/WA(x)) _E~ ^r((-/(x)) n A(x)) 



(e-/(*)> 



E~o^<(-/(x)) r 



(3.287) 



When computing a term of a given order n, such as ^ ((— / (x))™ A (x)), we may 
always write 



-((-/(x))"4(x)) = ^^^ ra+(; " 



n 



n\ eiml 



(-/(xV/A(x)) <(-/(x)n 



(3.288) 



where the subscript c on the average means that none of the terms in ^(— / (x)/ A (x) 
can be factored into lower order correlation functions, such as for example ^(— / (x))^ (A 

or ^(— / (x)/ -1 ^ ((— / (x)) A (x)) etc... The combinatorial factor corresponds to 

the number of ways the (— / (x)) n can be grouped into a group of I terms and a 
group of n — £ terms, the S n:7n+ g Kronecker delta function ensuring that indeed 
m = n — I. Using the last equation in the previous one, the sum over n is now 
trivially performed with the help of 5 n ,m+e and one is left with 



(3.289) 



(e-/WA(x)) _ E^oE^o^r((-/(x))M(x)) c ((-/(x)n 

< e - /(x) > ~ e:^((-/(x)d c 



The numerator can now be factored so as to cancel the denominator which proves 
the theorem 



Theorem 13 Linked cluster theorem for normalized averages: 



(e-'WA(x)) _^ 



E£ozr((-/(*))^M) =< e - /(x) ^(x)> 



(3.290) 



This result can be applied to our calculation of the Green's function since 
within the time-ordered product, the exponential may be expanded just as an ordi- 
nary exponential, and the quantity which plays the role of (— / (x)) , namely drV (t) 
can be moved within the T T product without costing any additional minus sign. 



Linked cluster theorem for characteristic functions or free energy 

We now wish to show the following general theorem for a multivariate probability 
distribution. 

Theorem 14 Linked cluster theorem (cumulant expansion). 



In (e-^ x )) = £^° =1 i ((-/ (x))™) c = (e- / ^) c - 1 



(3.291) 



The proof is inspired by Enz[24]. When / (x) = ik • x, the quantity (e~ lk ' x ) 
is called the characteristic function of the probability distribution. It is the gener- 
ating function for the moments. The quantities on the right-hand side, which as 
above are connected averages, are usually called cumulants in ordinary probability 
theory and ln(e~ lk ' x ) is the generating function for the cumulant averages. 
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Proof: To prove the theorem, we introduce first an auxiliary variable A 

^( 6 - A/(x) ) = ( e - A/(x) [-/(*)]) 
We can apply to the right-hand side the theorem we just proved 
(e-W h/ (x)]) = (e-^ [-/ (x)]) c (e-VW] 

so that 

1 9 /„-A/(x)\ _ / # -A/(x) 

(e-VW)d\\ I \dX 

Integrating both sides from to 1, we obtain 

ln(e- A K*>)|S = (e-'< x >) -1 

QED 



(3.292) 

(3.293) 
(3.294) 

(3.295) 



Example 15 It is instructive to check the meaning of the above result explicitly 
to second order 

In (e"*/(*) \ « l n A _ X f (x) + i (A/ (x)) 2 \ « f- (A/ (*)> + I ((A/ (x)) 2 \ V I (A/ 



(3.296) 
(3.297) 

(3.298) 



The above results will help us in the calculation of the free energy since we 
find, as in the first equations of the section on linked cluster theorems, 



(e- Xfix) ) c - 1 ~ — (A/ (x)) c + \ ((Xf (x)) 2 ) c 
so that equating powers of A, we find as expected, 

;(/(x)) 2 ) c = ((/(x)) 2 )-(/(x)) 2 . 



-Tin 



e -f£drV(T) 



OJ 



11=1 







drV (t) 



F = -TlnZ = -T 


[<T T 


'e-SSdrV(r) 


) 1 

1 Oc 


-TlnZ . 



TlnZ 

:99) 
(3.300) 



Oc 

(3.299) 



the subscript stands for averages with the non-interacting density matrix. The 
above proof applies to our case because the time-ordered product of an exponential 
behaves exactly like an ordinary exponential when differentiated, as we know from 
the differential equation that leads to its definition. 



3.7.3 Variational principle and application to Hartree-Fock theory 



It is legitimate to ask if there is a one-body Hamiltonian, in other words an 
effective Hamiltonian with a time-independent potential, whose solution is as close 
as possible to the true solution. To address this question, we also need to define 
what we mean by "as close as possible". The answer to both of these queries 
is provided by the variational principle for thermodynamic systems. We discuss 
below how Hartree-Fock theory comes out naturally from the variational principle. 
Also, it is an unavoidable starting point when there is a broken symmetry, as we 
will discuss more fully in a later chapter. 
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Thermodynamic variational principle for classical systems 

One can base the thermodynamic variational principle for classical systems on the 
inequality 

(3.301) 



e x > 1 



which is valid for all x, whether x > 0, or x < 0. This inequality is a convexity 
inequality which appears obvious when the two functions are plotted. We give two 
proofs. 

Proof 1: e x is a convex function, i.e. d 2 e x /d 2 x > for all values of x. At x = 
the functions e x and l + ias well as their first derivatives are equal. Since a 
straight line tangent to a convex curve at a point cannot intersect it anywhere 
else, the theorem is proven. QED 

Algebraically, the proof goes as follows. 

Proof 2: The equality occurs when x = 0. For x < —1, e x > while 1 + x < 0, 
hence the inequality is satisfied. For the remaining two intervals, notice that 
e x > 1 + x is equivalent to 



£ > 0. (3.302) 

n=2 



For x > 0, all terms in the sum are positive so the inequality is trivially 
satisfied. In the only remaining interval, —1 < x < 0, the odd powers of x in 
the infinite-sum version of the inequality are less than zero but the magnitude 
of each odd power of x is less than the magnitude of the preceding positive 
power of x, so the inequality (3.302) survives. QED 

Moving back to our initial purpose, let Hq be a trial Hamiltonian. Then take 
the trial density matrix corresponding to averages ()q. We will 
use the above inequality Eq. (3.301) to prove that 



-T\nZ < -TlnZ^ + {H - H 



(3.303) 



This inequality is a variational principle because H is arbitrary, meaning that 
we are free to parametrize it and then to minimize with respect to the set of 
all parameters to find the best one-particle Hamiltonian in our Physically chosen 
space of Hamiltonians. 

Proof Our general result for the free energy in terms of connected terms, Eq.(3.300), 
is obviously applicable to classical systems. The simplification that occurs 
there is that since all operators commute, we do not need to worry about 
the time-ordered product, thus with 

V = H-H (3.304) 

we have 

-TlnZ 5 . (3.305) 



F = -TlnZ = -T 



-0V\ 



e- pv >_ -1 

/ 0,c 



Using our basic inequality Eq.(3.301) for e ^ v we immediately obtain the 
desired result 



F <-T(-pvY +F T) (3.306) 

\ / 0,c 



which is just another way of rewriting Eq. (3.303). 
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It is useful to note that in the language of density matrices, 
the variational principle Eq. (3.303) reads, 



-T\nZ < Tr [p (H - pN)] + TTr [p In p ] 



(3.307) 



which looks as if we had the function (E — pN) — TS to minimize, quite a satis- 
factory state of affairs. 



Thermodynamic variational principle for quantum systems 

For quantum systems, the general result Eq. (3.303) applies but it is more difficult 
to prove because there is in general no basis that diagonalizes simultaneously 
each and every term in the expansion of T T exp — drV (r) . If V was not time 

dependent, as in the classical case^then matters would be different since V n would 
be diagonal in the same basis as V and one could apply our inequality Eq. (3.301) 
in this diagonal basis and prove the theorem. The proof of the variational principle 
in the quantum case is thus more complicated because of the non-commutation of 
operators. As far as I know, the proof is due to Feynman [25]. We do not repeat 
it here. 



Application of the variational principle to Hartree-Fock theory 

Writing down the most general one-body Hamiltonian with orthonormal eigen- 
functions left as variational parameters, the above variational principle leads to 
the usual Hartree-Fock eigenvalue equation. Such a general one-body Hamiltonian 
would look like 

^ = zZj d ^ to ("^) to c+ « c « ( 3 - 308 ) 

with (f> a (x) as variational wave-functions. In the minimization problem, one must 
add Lagrange multipliers to enforce the constraint that the wave-functions are not 
only orthogonal but also normalized. 

In a translationally invariant system, the one-body wave functions will be plane 
waves usually, so only the eigenenergies need to be found. This will be done in 
the following chapter. 

It does happen however that symmetry is spontaneously broken. For example, 
in an anti-ferromagnet the periodicity is halved so that the Hartree-Fock equa- 
tions will correspond to solving a 2 x 2 matrix, even when Fourier transforms are 
used. The matrix becomes larger and larger as we allow more and more general 
non-translationally invariant states. In the extreme case, the wave functions are 
different on every site! This is certainly the case in ordinary Chemistry with small 
molecules or atoms! 



3.8 Quantum impurities 



An important example that we will encounter in this School is that of a cor- 
related impurity in a bath of non-interacting electrons. Dynamical Mean-Field 
theory is based on this example. We will only set up the problem without solving 
it. The Numerical Renormalization Group approach (NRG) and Density Matrix 
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Renormalization Group are examples of approaches that can be used to solve this 
problem. 

We begin with the Anderson impurity problem. Including the chemical poten- 
tial we obtain, 

Ki = Hf + H c + Hf c - fiN (3.309) 
H f = ^2(e-^flf ia + u[fy^ (3.310) 

a 

H c = ^2^2(£k- fj)cl a c ka (3.311) 

a k 

H fc = 2XVaccL/ia + fr.c. (3.312) 

a k 

Note that the sum over k in the hybridization part of the Hamiltonian Hf c basi- 
cally tells us that it is the local overlap of the conduction band with the impurity 
that produces the coupling. 

Suppose we want to know the properties of the impurities, such as the local 
density of states. It can be obtained from the Green function 

Gff(r) = -(T T f ia (r)fl). (3.313) 



We will proceed with the equations of motion method. We first write the 
equations of motion for Ck Q and fi a 

d 

— c kQ = [H I ,c ka ] (3.314) 

= - (e k - M) c ka + Vikfia (3.315) 
g^fia = [HjJia] (3.316) 
= -(e-tifia-Uflji-afia + V&Cka (3.317) 

Proceeding like our in our earlier derivation of the equations of motion Eq.(3.153) 
we have 

§pSff(r) = S(r)({f ia (r),fl}) 

- (T T (- (e - A*) fia (r) - Ufl a (r) f^ a (r) f ia (r) + V^cJ) fl ) (3-318) 
= -S (r) - (e - M ) Gf S (r) + U (T T fj_ a (r) f^ a (r) f ia (r) f} a ) + V^G cf (k, i, r) 



where we defined 

G cf (k, i, r) = - (T T c kQ (r) fj a ) . (3.319) 



To eliminate this quantity, we write its equations of motion 

j^ c/ ( M ,r) = S(r)({c ka (r)jl}) 

- (t t (- (£k - fi) c ktt (t) + Vik/fa (r)) /t, ) 
= -*(T)-(ek-/i)ec/(k,*,r) + Vike//(T) (3.320) 

It can be solved by going to Matsubara frequencies 

G c f (k, i, ik n ) = ^V ik g f f (ik n ) . (3.321) 

ik n -(e k -fj,) 
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Substituting in the equation for Qff (ik n ) we obtain 

ik n - (e - (jl) + V£- t -Vn 



Q ff (ik n ) 



''ik n - (£k - M) 

- l ~ U [ ( T rfl a (r) fi-a (T) fia (t) £ 



(3.322) 



The equation to be solved has exactly the same structure as we would find for a 
single impurity except that now the "non-interacting" Green function is 



g° ff (ik n ) 1 =ik n -(e- ri + v* k 



1 



ik n - (£k - M) 



(3.323) 



This is in fact exactly the non-interacting Green function that we would find with 
U = 0. One can propagate from the impurity site back to the impurity site by 
going through the bath. One often defines the hybridization function A// (ik n ) 

by 

A ff (ik n ) = Vl . kn _l £k _ fi) V ik . (3.324) 
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